Numerical Simulations of GFRP-Reinforced Columns Having Polypropylene and Polyvinyl Alcohol Fibers

.e present investigation aims to propose a numerical model for assessing the complex damaging response of glass fiberreinforced polymer(GFRP-) reinforced concrete columns having hybrid fibers and confined with GFRP spirals (GFHF columns) under concentric and eccentric compression. Fiber-reinforced concrete (FRC) consists of polyvinyl alcohol fibers (PVA) and polypropylene fibers (PF). A total of six GFHF circular columns were constructed having a circular cross section of 250mm and a height of 1200mm. A commercial package ABAQUSwas used for the finite element analysis (FEA) of the GFHF columns by using a modified concrete damage plastic (CDP) model for hybrid fiber-reinforced concrete (HFRC). .e damaging response of GFRP bars was defined using a linear elastic model..e results depicted that the failure of GFHF columns occurred either in the upper or in the lower half portion with the rupture of GFRP longitudinal bars and GFRP spirals. .e decrease in the pitch of GFRP spirals led to an improvement in the axial strength (AS) of GFHF columns..e eccentric loading caused a significant reduction in the AS of columns. .e comparative study solidly substantiates the validity and applicability of the newly developed FEA models for capturing the AS of GFHF columns by considering the axial involvement of longitudinal GFRP bars and the confinement effect of transverse GFRP spirals. So, the suggested numerical model having a complex system of equations for HFRC can be used for the accurate analysis of HFRC members.


Introduction
e tendency of steel subjected to corrosion results in a substantial decrease in the service life of concrete structures made of steel reinforcement by lowering their strength and ductility. To avoid corrosion of steel bars, it is necessary to apply certain protective practices such as epoxy coatings and sacrificial cathodes which leads to additional cost. Glass fiber-reinforced polymer (GFRP) material became the most potential substitute having characteristics like lightweight, low density, high resistance, low thermal conductivity, electromagnetic susceptibility, and, most specially, high resistance to corrosion [1][2][3][4]. ese anisotropic composite materials are generally more effective in aggressive and corrosive conditions, to minimize the running costs while increasing the service life of reinforced concrete [2,5]. e decrease in the brittleness of plain concrete is the focus of advanced research. erefore, polyvinyl alcohol (PVA) fibers and polypropylene fibers (PF) were added into the concrete in the current work to improve its ductility by reducing the brittleness. Besides, previous studies portrayed that the concrete columns reinforced with steel bars and fibers portrayed increased strength and ductility [6,7]. ere is a need to investigate the performance of reinforced columns with GFRP bars and hybrid fibers by performing experimentation and numerical simulations which is the main goal of the present work.
Due to its superior performance, the use of GFRP bars in reinforced concrete (RC) structures over the last few decades has attracted the attention of many researchers.
For GFRP-RC slender columns, a greater lateral deflection was observed when compared with identical steel-RC columns [8].
e failure of GFRP lateral reinforcement typically occurs as a result of the bond-slip mechanism rather than rupture [9]. GFRP-RC columns have shown no balance points for the moment due to their linear elastic behavior [10]. e steel and GFRP-RC compression members demonstrated a similar behavior, except that the AS was found to be 7 percent lower for GFRP-RC columns [11]. e laterally confined GFRP-RC columns at a pitch of 76 mm showed the AS of 84 percent of their counterparts in steel-RC columns [12]. e GFRP bars perform well in concrete under compression due to the lower elastic modulus of GFRP bars compared with that of steel bars, and thus it is beneficial to use GFRP bars in concrete columns [13]. e axial and bending strengths were decreased, but column ductility was improved by utilizing the steel reinforcement with an equivalent amount of GFRP bars under various loading conditions [14]. e ductility of GFRP-RC columns has been improved substantially by decreasing the vertical spacing of stirrups [15][16][17]. e numerical simulations have become the interest of advanced research works. Different numerical models have been applied in different fields of sciences [18][19][20][21]. e finite element models (FEMs) are considered as the most effective methods for identifying the dynamic damaging behavior of fiber-reinforced polymer (FRP) composites and their contact mechanism with the concrete at minimal cost and time [22,23]. Several studies that carry out the finite element modeling of GFRP-RC columns could be found in the literature [24][25][26][27][28][29]. Using the concrete damaged plastic (CDP) model, the confinement of concrete compression members can be precisely described by the finite element analysis [30]. A FEM was introduced using ABAQUS which forecasted the axial response of the steel-RC and GFRP-RC columns accurately [2]. Simulation models of the finite elements accurately predicted the buckling behavior of pultruded GFRP columns [31]. Other models and simulation techniques have been used to predict the performance of output variables [32][33][34][35]. ANSYS analyzed various slender ratios for GFRP-RC short columns, and a close correlation was noted among tests and FEM estimates for failure mechanisms and postpeak collapse behavior [36,37]. In ABAQUS, the bond-slip behavior between the concrete and the FRP bars was accurately validated to obtain a strong correlation of tests with the simulated results in the postyield phases [25,38]. FEM can model accurately the failure modes, peak loads, corresponding deflections, and load-deflection behavior of GFRP-RC columns [2,24,39]. e structural efficiency of GFRP-RC columns with hybrid fibers (GFHF columns) has not been investigated by any researcher. Also, the finite element analysis (FEA) of these structural elements is scarce in the earlier findings. None of the previous works proposed a novel FEM for predicting the complex damage response of GFHF columns by considering a modified plasticity model for the HFRC. So, the objective of the current study is to quantitatively measure the structural behavior of corrosion resistant GFHF columns using the experimental data and numerical simulations of FEM with different pitches of GFRP spiral under different loading conditions. e simulation model of the structural performance of GFHF columns was presented due to the hybrid fibers reinforced concrete (HFRC), lateral GFRP spiral pitch, and diverse eccentric loadings. Using ABAQUS software, the FEM of all GFHF columns was done. e HFRC behavior was visualized using a modified CDP model, and the GFRP bar behavior was described by assuming that it was a linear elastic material.
is research will play a pivotal role in understanding the behavior of GFHF columns, experimentally as well as numerically, and putting them into practice in the construction sector.

Materials
In this work, good-grade OPC following the guidelines of ASTM C150 [40] was used for the production of HFRC. e fineness modulus of fine aggregate was 2.3, while the coarse aggregate of maximum size 10 mm was used. PF and PVA fibers lengths were 12 ± 1 mm and 10 ± 1 mm, and their diameters were 24 μm and 20 μm. A slump value of 90 mm was observed for the fresh concrete as per ASTM C143 [41]. Diverse HFRC ingredients are presented in Table 1. On the day of column construction, six cylinders with dimensions of 150 mm × 300 mm were produced and tested on the same day when specimens were examined. e average compression strength of the cylinders was 38 MPa when measured according to ASTM C39 [42] with a deviation of 2.24 MPa.

GFRP Bars.
Six longitudinal GFRP bars with a diameter of 14 mm and GFRP spirals with a diameter of 10 mm were used to strengthen the GFHF columns. e GFRP longitudinal bars and spirals were made utilizing E-glass fibers impregnated with additives, fillers, and thermosetting vinyl ester resin having 80 percent fiber quantity by volume. e physical and mechanical properties of GFRP bars are as shown in Table 2.

Column Construction.
A total of six circular HFRC columns with longitudinal GFRP bars and transverse GFRP spirals were built in this study. Under concentric and eccentric loading conditions, these columns were tested to failure to analyze the effect of hybrid fibers, the pitch of spirals, and the load eccentricity on the structural performance of the specimens. All the columns had a diameter of 250 mm and a height of 1200 mm. ese column measurements were used keeping in mind the size and compression efficiency of the testing unit. ese measurements also render the columns large enough to be treated as a fullsize test specimen.
All GFHF columns comprised six 14 mm longitudinal GFRP bars providing a 1.88% longitudinal reinforcement ratio. Usually, such a low reinforcement ratio is provided in the areas with no seismicity zones. e transverse 2 Complexity confinement was supported employing GFRP spirals with a diameter of 10 mm at a pitch of 38 mm and 75 mm, giving the volumetric ratios of 2.88% and 1.46%, respectively. e pitch of the GFRP spirals was prescribed to prevent GFRP bars from buckling elastically [43]. For all the columns, a concrete cover of 20 mm was provided. Figure 1 demonstrates the geometrical structure of a specimen with a spiral pitch of 38 mm, and the geometrical and test descriptions of all manufactured specimens are described in Table 3. e circular PVC pipes have been used as formwork, with an inner diameter of 250 mm and a wall thickness of 8 mm.
e reinforcement cages were positioned in the formwork, and HFRC was added and vibrated nonstop using an electrical vibrator during the construction activity. All the specimens were treated at room temperature and protected with sheets of polyethylene to prevent moisture loss.

Testing and Instrumentation.
In the testing machine, the testing of the samples was performed at a loading rate of 0.002 mm/sec. e specimens were prepared using steel collars (100 mm long and 10 mm thick) and plastered on the top and bottom surfaces before testing such that end crushing of the columns could be prevented and the load could be spread equally across the cross section. By applying a compression load on the top flattened surface, the concentric columns were loaded uniformly, whereas the eccentric columns were loaded by applying a line load at the appropriate eccentricity using a steel rod positioned between the loading plates [44] as shown in Figure 2. Readings of the axial load, axial deflection, and axial strains of the specimens were recorded using a data logger connected to the test compression unit. e failure patterns and failure modes of specimens were recorded using a video camera.

General Methodology.
is section specifically defines the FEA of GFHF columns by detailed simulation of HFRC and GFRP reinforcement behavior. For the FEA, a wellknown commercial software ABAQUS 6.14 was used. e HFRC simulations were produced by modeling it as a portion of 3D solid stress. e reinforcing bars were interpreted as 3D deformable wire elements. Modified concrete damage plasticity (CDP) model was used to describe the HFRC damage behavior that considers the influence of PF and PVA fibers. e bottom end of samples was restricted, and the top end was permitted for all degrees of freedom. e relationship between the reinforcement and HFRC was modeled by the use of "embedded area" constraint, which combines the compatible degrees of freedom of the wire elements with the requisite degrees of freedom of the concrete 3D stress elements [2]. Using the displacement control method, a uniform load of 15 mm on the upper surface of concentric columns was enforced while a line load of 15 mm was provided at the top of the eccentric specimens at the intended eccentricity keeping the boundary condition identical to that applied in the experimental study. For uniform load distribution, a steel plate with a thickness of 50 mm was tied on the top and bottom surfaces of the specimens. e interaction between the steel plate and the specimen surface was established using the ABAQUS "tie" constraints. e model specimens with the geometric and applied boundary conditions are shown in Figure 3.

Simulation of Plasticity, Tension, and Compression of HFRC.
e accurate simulation of HFRC behavior is very significant for predicting the behavior of GFHF columns. e HFRC elastic behavior can be predicted with two parameters: Poisson's ratio assumed to be 0.2 [45] and the elastic modulus computed as 4734 MPa is the compression strength of HFRC. e nonlinear and irreversible detrimental behavior of concrete can be precisely modeled using the CDP model, which considers the isotropic compression and tensile elasticity and plasticity of isotropic damage [24,47]. Hence, this model with some modifications has been used in the present research to describe the damaging behavior of complex HFRC.
In the CDP model, HFRC plasticity behavior can be described by five parameters: the eccentricity (e), the shape factor for yielding of the surface (K hf c ), the viscosity parameter (μ), the dilation angle of HFRC (ψ), and the ratio of biaxial to uniaxial compression yield stresses (σ hf bo /σ hf co ). All these variables for HFRC were changed excluding the parameter μ which was adjusted to obtain close predictions with the samples. e factor "e" has a default value of 0.1 for reinforced concrete in both plain and fiber [48]. us, this value was used in the present research for the parameter "e". For plain concrete, the value of the parameter (K c ) is always between 0.64 and 0.80 [49]. A strong estimate for the strength of plain concrete is obtained under low hydrostatic stresses when a value of 0.67 is used for K c . Additionally, an appropriate prediction is obtained under high hydrostatic  stresses when a value of 0.70 is used for this variable [48].
Hence, when dealing with HFRC, the relationship for K hf c can be reported by presenting some enhancing constants for compression and tensile meridians, i.e., k c and k t [45,48]: where k c and k t represent the constants for the compression and tensile meridians of the concrete. e model for the constant k c HFRC has not been well-defined yet due to a deficiency of experimental testing in this area. However, the constant k c was defined for PVA fiber-reinforced concrete that can be reported by using equation (2) and the constant k t for the HFRC can be reported by using equation (3) [45]: where V pf is the percent portion of PF by volume, V pvaf is the percent portion of PVA by volume, l pf is the total length of PF, l pvaf is the total length of PVA, d pf is the diameter of PF, d pvaf is the diameter of PVA fibers, λ pvaf is the fiber reinforcement index for PVA, and λ pf is the fiber reinforcement index for PF. Due to the prevention of cracks development and propagation in the concrete by adding the hybrid fibers, the enhancement in the tensile meridian will be higher than that of the compression meridian, i.e.,   4 Complexity k t /k c ≥ 1. erefore, inequality K hf c ≥ K c will always clutch for the HFRC. Making the amendments for HFRC, the parameter σ hf bo /σ hf co can be theoretically measured by using equation (6) to be placed in the CDP model [48]: e numerical value for the parameter σ hf bo /σ hf co was obtained as 1.46 after placing the values of enhancing constants (k c and k t ) in equation (6). When we use the values of enhancing constants as a unit representing the plain concrete, the value of 1.16 is secured for the parameter σ hf bo /σ hf co . e flow rule is measured by the parameter ψ that is very useful for the accurate simulations of the concrete in ABAQUS. e value of ψ will be less for HFRC as related to that of the plain concrete. e recommended value of ψ locates between 31°and 42°for the plain concrete [50]. A theoretical model for the dilation angle of HFRC (ψ hf ) can be reported in terms of λ pvaf and λ pf [48]: According to the elastic-plastic theory, the total strain (ε) of concrete produced due to the application of load can be reported in two parts: one is the elastic strain (ε el ) and second is the plastic strain (ε pl ) of concrete: e nonlinearity of the concrete is associated with plastic, compression, and tensile damage behaviors of concrete. e strength and stiffness reduction of the concrete increase when the strain produced in the concrete increases because the damaging phenomenon is associated with the elastic and plastic parts of strain, i.e., ε el and ε pl [48]. e damage mechanism of HFRC can be fully defined in numerical simulations by two parameters: one is the compression damage parameter (d c ) and the second is tension damage parameter (d t ). ese two parameters are used in the CDP model for the definition of damages produced by compression and tensile loads. Figure 4 reports the complete general stress-strain behavior of concrete subjected to compression loads. e compression stress of concrete (σ c ) can be expressed by the following equation: In this model, E o indicates the modulus of elasticity of HFRC that can be measured using the equation provided by ACI 318-95 [46], ε c indicates the compression strain, and ε pl c is the plastic part of the compression strain of HFRC. e parameter d c was theoretically measured by using equation (10) as recommended by Wang and Chen [51]: In this model, the parameter m c is the controlling parameter for the speed of damaging evolution of concrete under compression having a value of 0.1 [52] and ε in c,norm is the standardized inelastic strain of concrete exposed to compression that can be denoted as ε in c /ε in cu , where ε in cu is the strain equivalent to the standardized inelastic strain of concrete exposed to compression with a value of 0.033 [52]. us, the improved controlling parameter for HFRC under compression loads (m hf c ) can be reported by the following equation [48]: where a m1 and b m1 depict the constants that are associated with the characteristic parameter of different fibers with the For determining the maximum compression stress (σ hf co ) and the corresponding compression strain (ε hf co ) of the HFRC, the following models were employed in the present work [48]: ε hf co � ε co 1 + 0.705λ pvaf + 0.364λ pf .
e postcracking and propagation of cracks are reduced, and the ductility of concrete is improved by the addition of fibers. e postcracking performance of HFRC can be welldefined by the strain-hardening and a strain-softening portion of the stress-strain curve. e behavior of concrete before achieving the peak load can be well-defined by the strain-hardening, and the strain-softening defines the postcollapse behavior [53]. Figure 5 reports the complete stress-strain performance of HFRC under tensile loading from which the models for the compression stress of HFRC can be presented in the following form: where ε t depicts the tensile strain of HFRC and ε pl t is the plastic parts of the tensile strain of HFRC. e constant d t was well-defined by using equation (17) as recommended by Wang and Chen [51]: In this model, the parameter m t is the controlling parameter for the speed of damaging evolution of concrete under tension having a value of 0.05 [52] and ε ck t,norm is the standardized inelastic strain of HFRC exposed to tension that can be reported as ε ck t /ε ck tu , where, ε ck tu is the strain corresponding to the standardized inelastic strain of HFRC exposed to tension with a value of 0.0033 [52]. e improved controlling constant for HFRC in tension (m hf t ) can be reported by the following equation [48]: where a m2 and b m2 depict the constants that are linked to the characteristic parameter of PVA and PF with the values of 0.628 and 0.156, respectively [48]. erefore, the damage parameter in tension d hf t and tensile stress (σ hf t ) for the HFRC can be reported by the following equations, respectively: e maximum tensile stress (σ hf to ) and the corresponding tensile strain (ε hf to ) of the HFRC can be theoretically measured using the following models [48]: σ hf to � σ to 1 + 0.379λ pvaf + 0.020λ pf , ε hf to � ε to 1 + 0.498λ pvaf + 0.697λ pf . 6 Complexity perfect bond between the HFRC and GFRP bars was defined using the 'embedded region' constraint present in ABAQUS. Different properties of the reinforcement were taken according to Table 2. Poisson's ratio of GFRP material was considered as 0.25 [28], and a linear elastic model up to failure for employed for the definition of GFRP bars as recommended by the literature [2,25,28]. e linear elastic behavior of GFRP bars is reported in Figure 6. e GFRP bars were considered having their compressive strength as half of the tensile strength [54].

Calibration of the Finite Element Model.
For the finite element simulations, it is essential to validate the predictions using experimental results. In the present study, a control model (AA38-0ECC) was selected whose load-defection performance was employed for the calibration and validation purpose of the control FEM. e calibration process of the control FEM was performed for various geometrical and material parameters of GFHF columns. Different boundary conditions were also employed to select the most suitable constraints for the accurate predictions. e calibration for the plastic behavior of HFRC included the viscosity parameter (μ) only. All the other parameters of HFRC were theoretically measured using the modified equations provided in the previous section of this study. Various mesh sizes and the different types of elements were also assessed to secure the most accurate predictions compared with the testing measurements. When the control FEM was calibrated for different boundary conditions, the most accurate predictions were obtained for fixing all the degrees of freedom at the bottom and the top end was kept free for all the translational and rotational degrees of freedom. e different values used to examine the influence of μ on the load-deflection performance of control specimen (AA38-0ECC) were 0.0096, 0.0080, 0.0064, 0.0048, 0.0032, and 0.0016. For securing the close agreement with the testing results, the value of this parameter should be small [55]. e influence of μ on the load-deflection performance of the control specimen was reported in Figure 7(a). e value of 0.0032 for μ portrayed a close correlation with the testing results.
For avoiding the discontinuity of strain and stress spreading and to obtain accurate numerical results, the sensitivity analysis of the load-deflection performance of the control specimen due to the mesh size was carried out.
e size of the elements should be kept small enough to give accuracy and large enough to avoid the numerical convergence issues. Figure 7(b) reports the sensitivity analysis of load-deflection performance due to an investigation of mesh sizes of 70 mm, 60 mm, 50 mm, 40 mm, 30 mm, 20 mm, and 10 mm. e close agreement was found using a 20 mm size of elements. e studied elements of the HFRC were triangular, hexahedral elements, and tetrahedral elements available in ABAQUS. e best performance of the control model was obtained when the C3D8R elements for HFRC and T3D2R elements for the GFRP bars and spirals were employed as evidenced from the literature [2,24,48].

Column Strength of Specimens.
e column AA38-0ECC portrayed the ultimate axial compression strength with a value of 2678 kN. e axial compression strength portrayed by the column with a double spacing of spirals (AA75-0ECC) was 2403 kN. is value was 10.2% lower than the axial compression strength of the specimen with a half spacing of the spirals. It depicts that the decrease in the axial compression strength is larger for the GFHF columns with the enhancement in the spacing of GFRP spirals. Consequently, the GFHF columns with a lesser pitch of GFRP spirals or with larger transverse confinement are more effective in terms of axial compression strength. e axial compression strength is meaningfully exaggerated by the application of testing eccentricity. e average percent reduction in the axial compression strength was 24.8% and 58.2% for the GFHF columns with the testing eccentricity of 30 mm and 60 mm, respectively.
e axial compression strength of column AA38-30ECC was 2.7% larger than AA75-30ECC, but the column with 60 mm eccentricity portrayed a conflicting behavior. e axial compression strength was larger for the column with a higher spacing of GFRP spirals, i.e., the axial compression strength of AA75-60ECC was 10% larger than the axial compression strength of AA38-60ECC.
is depicts that the GFHF specimens under comparatively higher testing eccentricities are more effective in terms of axial compression strength when the spacing of GFRP spirals is higher.

Failure Modes.
When the axial compression load was applied on the top surface of the columns, their behavior was linear elastic up to 90% of the peak load. At this stage, there were no cracks developed on the surface of the HFRC, and the lateral confinement mechanism of GFRP spirals was not activated. When the load was increased form the peak load, the cracking of the HFRC cover was initiated by sparking the small vertical hairline cracks on the compression portion of the specimens producing a small sound. e vertical hairline cracks spread vertically, and their width improved with the application of a more enhanced load. After securing the maximum axial compression strength of the specimens, the rate of improvement in the axial strain was lifted by forming the long vertical cracks along with the height of the specimens. At this stage, the development in the HFRC cover happened, and the restriction influence provided by the lateral GFRP spirals was initiated. After getting to about 70% of the ultimate axial compression load in the postpeak failure stage, the transverse GFRP spirals achieved their maximum axial compression strength producing a sound of rupture. Most of the specimens failed in the upper or lower half regions. During the damaging phenomenon, GFRP bars were fractured in GFHF specimens. For the columns with a lesser spacing of spirals (AA38-0ECC, AA38-60ECC, and AA38-30ECC), the failure happened due to the rupturing of GFRP spirals and the crushing of HFRC core. e performance of eccentric columns was similar. e axial compression strength of the columns improved constantly until the transmission of hairline vertical cracks on the compressive portion of the columns happened. Mainly, the failure of eccentric specimens occurred due to the crushing of HFRC on the compression portion of the columns. e tension side of these columns depicted the flexural-tension failure mechanism that also caused the decrease of the axial compression strength of specimens. No serious spalling of the HFRC cover was observed after the failure of columns because of the bridging influence due to the incorporation of PVA and PF fibers in the concrete which avoid the cracking and decrease the crack spreading.

Effect of Pitch of Spirals.
As shown in earlier works, the AS of the columns was enhanced with the reduction in the pitch of lateral confinement [44,56]. e same pattern of AS reduction was found for concentric and eccentric specimens. As the pitch of the spirals decreased from 75 mm (a volumetric ratio of 1.46 percent) to 38 mm (a volumetric ratio of 2.88 percent), the increase was 4.52 percent in AS of the concentric column.
is improvement in the axial compression strength is due to the good confinement and restrained GFRP bars and effective lateral detention of the HFRC material by the tightly spaced spirals to take in more energy [56]. As the pitch of the spirals increased, the axial deflection contributing to the ultimate AS enhanced. When the spiral pitch was raised from 38 mm to 75 mm, the axial deflection enhancement equivalent to the ultimate AS was 0.58%. Figure 8 shows the pitch effect of spirals. e loading eccentricity of the GFHF specimen of 30 mm showed no substantial difference in the load-deflection response and an increase in the spiral pitch from 38 mm to 75 mm, whereas the GFHF specimen exposed eccentricity of 60 mm showed higher AS and high associated deflection at the larger pitch of spirals demonstrating that the GFHF columns are more effective under loading eccentricity with a greater pitch of GFRP spirals.

Effect of Loading Eccentricity.
Due to the application of loading eccentricity, a substantial reduction in the AS of all columns occurred. Figure 9 demonstrates the influence of loading eccentricity on the AS of the specimens tested. e AS of the GFHF columns with a spiral pitch of 38 mm was decreased from 2678 kN to 1933 kN and 1005 kN indicating a decrease of 27.82 percent and 62.46 percent in the AS when the eccentricity was increased from 0 mm to 30 mm and 0 mm to 60 mm. Furthermore, the AS of the GFHF columns with a spiral pitch of 75 mm decreased from 2403 kN to 1879 kN and 1106 kN showing a decrease of 21.80 percent and 53.97 percent in AS while the eccentricity increased from 0 mm to 30 mm and 0 mm to 60 mm, respectively. erefore, due to geometric imperfections and design defects, the minor eccentricities produced should be minimized because they cause a  Table 4 presents the observed data and the FEM estimates for the ultimate AS and subsequent deflection of all the GFHF specimens tested. e FEM of the 30 mm eccentric GFHF column (AA38-30ECC) showed a minimum difference of 1.22 percent concerning the AS. For the specimens AA75-30ECC and AA38-60ECC, respectively, the maximum and minimum disparities for axial deflection compared concerning the ultimate axial force were noticed. Such high percentages of errors may be due to the initial geometric inconsistency of the specimens during the manufacturing process which was not measured in the current study. e average error for the GFHF test specimen was 2.44 percent and 3.71 percent, respectively, for the axial force and associated axial deflection. e presented FEM mostly misjudged the ultimate AS for columns of the GFHF. Figure 10 illustrates a comparative study of the experimental and FEA analyses for the ultimate AS.

Peak Strength and Corresponding Deflection.
As for the eccentrically loaded columns, the proposed FEM correctly described the sensitivity of load-deflection. e average difference between the results obtained and the FEM responses for the ultimate strength and related deflection of the eccentric columns was 2.18% and 4.56%. ese minor differences represent that the modeling of HFRC's tensile behavior was done exactly in the current study, and the currently proposed FEM is sufficiently accurate to simulate GFHF column behavior.

Complete Load-Deflection Curves.
e maximum loaddeflection curves of the GFHF columns estimated from experiments and simulated from FEM are shown in Figure 11. ese curves demonstrate that the proposed FEM estimated accurately the behavior of GFHF column load-deflection in the elastic range with some variance in the postpeak collapse point. e FEM of the AA75-30ECC and AA75-60ECC specimens reported stiffer responses in the elastic limit, whereas the FEM of all the specimens exceeded the postpeak performance. e variations in the postpeak collapse stage shown by the FEM results can be due to the concept of HFRC degradation and damaging modes and the constraints supposed between the reinforcing bars and HFRC.
rough working on certain    parameters in detail, the findings can be further improved. However, the proposed FEM successfully estimated the overall response of the load-deflection curve of GFHF columns. Figure 12 presents the damage quantification of cracks based on experimental study and FEM simulations. e crack trends in finite element models can be illustrated precisely by the   positive principal plastic strain since the path of the cracks in concrete is perpendicular to this sort of strain [2,30,57,58]. us, in the present study, the simulation of GFHF column failure modes and cracking behaviors was done using such strains. It was concluded that the proposed FEM modeled all the specimens' failure modes and cracking trends in a well-conforming way with the experimental work. After obtaining the ultimate AS, concrete cover spalling was undertaken together with reinforcement yielding. e failure mainly occurred in the upper portion of the specimens which was correctly predicted by the FEM. In close correlation with the test results, the proposed FEM also forecast the failure mechanisms and cracking behavior of the eccentrically loaded columns. e presumption of a good bond between HFRC and GFRP is strong compared with that of HFRC steel reinforcement. GFRP bars have lower elasticity modules like that of HFRC than steel bars. erefore, strong alignment between the GFRP reinforcing bars and HFRC can be established.

Conclusions
e present research was carried out to investigate the structural performance of GFHF columns by conducting tests and numerical simulations. e numerical simulations were conducted using a modified damaged plastic model to predict the complex damaging behavior of HFRC. From this work, the following key points could be drawn: (1) e experimental tests and numerical simulations based on the complex equations depicted that the identical modes of failure for all the GFHF columns were observed. e failure mechanism either in the upper half section or in the lower half section with a ruptured sound of longitudinal GFRP bars and spirals of GFRP after reaching their ultimate strength.
(2) e loading eccentricity induced the AS of GFHF columns to decrease significantly. Owing to loading eccentricities of 30 mm and 60 mm, the average reduction in the AS of GFHF columns was 25 percent and 58 percent, respectively. (3) e reduction in the pitch of spirals led to an increase in the AS of GFHF column. As the pitch of the GFRP spiral dropped from 75 mm to 38 mm, the AS of concentrically loaded GFHF columns increased by 10 percent. (4) e finite element model had been proposed using a modified CDP model for the concrete having hybrid fibers. e relative difference between the FEM estimates and the experimental results of GFHF column was 2.4 percent for AS and 3.7 percent for the deflections at ultimate AS. e proposed FEM solidly exhibited its applicability to effectively describe the full load-deflection curves and the failure modes of specimens. us, the comparative studies firmly confirm the effectiveness of the proposed numerical models for forecasting the performance of GFHF columns using the suggested complex system of equations.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors have no conflicts of interest.