kian-et-al-seismic-testing-and-modeling-of-full-scale-substandard-rc-columns-retrofitted-with-sprayed-gfrm-with-and
kian-et-al-seismic-testing-and-modeling-of-full-scale-substandard-rc-columns-retrofitted-with-sprayed-gfrm-with-and
Nima Kian, S.M.ASCE1; Ugur Demir2; Ali Osman Ates3; Oguz C. Celik, A.M.ASCE4;
and Alper Ilki, Aff.M.ASCE5
Abstract: This study presents the experimental and analytical hysteretic behaviors of eight full-scale RC square and rectangular columns.
The columns were designed to have different shear spans that represent: (1) a column that complies with the Turkish Government Ministry of
Reconstruction and Resettlement’s (1975) seismic design code, Turkish Seismic Design Code (TSDC); (2) a substandard column; and (3) two
sprayed glass fiber–reinforced mortar (GFRM)-retrofitted counterparts of the substandard column with and without basalt mesh. The substan-
dard columns were designed to be subjected to relatively high shear ratios (i.e., the ratio of the shear force that corresponds to the moment
capacity to shear strength of the cross section) up to 0.85 and with a high axial load-to-capacity ratio of 0.75. All columns were tested under
constant axial load and reversed cyclic lateral displacement excursions. The results revealed that the columns that complied with the TSDC
showed satisfactory behavior for seismic performance, and the performance of the substandard columns was extremely poor. However, the
hysteretic performance of the substandard columns that were subjected to high axial stress and shear significantly improved after the proposed
retrofitting. Finally, a numerical model was developed in OpenSees to reproduce the hysteresis curves of the specimens. The slip of the
longitudinal bars at the column–foundation interface, strain penetration into the foundation, and buckling of the longitudinal bars in com-
pression were accounted for in the modeling. The results are in good agreement with the experimental hysteresis curves. The performance
levels of the columns are further specified, and the predictions of the current seismic codes were analyzed: (1) the European Committee for
Standardization’s 2005 code, Eurocode 8: Design of structures for earthquake resistance; Parts 1–3: Strengthening and repair of buildings
(EC8-3); and (2) the Turkish Government Ministry of Interior Disaster and Emergency Management Authority’s 2018 code, Turkish Building
Earthquake Code (TBEC). The TBEC provided more accurate estimates of plastic rotation capacities for substandard specimens. In contrast,
EC8-3 overestimated the plastic rotation capacity when shear stresses were relatively high due to lower shear span-to-depth ratios (a/d).
DOI: 10.1061/JCCOF2.CCENG-5227. © 2025 American Society of Civil Engineers.
Author keywords: Basalt mesh; Column; Ductility; High axial load; OpenSees; Reinforced concrete; Retrofitting; Seismic performance;
Sprayed glass fiber–reinforced matrix; Substandard.
(Ansal et al. 2009), an urgent need for retrofitting is inevitable. tively (where the shear span is a, and the effective depth of the cross
Some studies discussed major earthquakes and their potential ef- section is d ). Another originality of this study is the continuation of
fects on substandard buildings and proposed rapid risk assessment the tests until the lateral load resistance or axial load capacity is
techniques (Celep et al. 2011; Tapan et al. 2013; Aydogdu and Ilki lost. In the literature, most tests on RC columns have been termi-
2023; Gurbuz et al. 2023; Aydogdu et al. 2023). Some studies have nated soon after a 15%–20% loss in lateral load capacity. In addi-
enhanced the seismic performance of substandard members using tion, although CNR-DT215 (CNR [Italian National Research
different retrofitting methods in the last few decades (Dai et al. Council] 2018) and ACI 549 (ACI [American Concrete Institute]
2012; Fakharifar et al. 2015; Deng et al. 2018; Lu et al. 2019; 2020) have addressed the design of FRCMs under gravitational
Raza et al. 2019), where different applications that used fiber- loads, there is no information on the seismic retrofitting of existing
reinforced polymer (FRP) composites have been more preferable structural members with a composite material such as FRCM in the
in recent years (Ilki et al. 1998; Mirmiran et al. 1998; Ilki and Kum- regulations and standards, to the best of the authors’ knowledge.
basar 2001; Bousias et al. 2004; Ilki et al. 2004; Farrokh Ghatte Therefore, in this study, the efficiency of an FRCM is investigated
et al. 2019; Kian et al. 2021; Narlitepe et al. 2021; Tore et al. under pure axial stress and bending moments. In addition, in the
2021; Kian et al. 2024). Of note, the FRP method has some draw- section “Analytical Modeling”, to simulate the experimental behav-
backs, such as the emission of harmful greenhouse gasses during ior, a numerical model is developed in the OpenSees framework
the gluing process, the poor behavior of epoxy resins at tempera- (Mazzoni et al. 2006) to predict the hysteretic load–displacement
tures above the glass transition temperature, ineffectiveness on curves of all tested specimens that consider the buckling and slip
wet surfaces or at low temperatures, the difficulty of postearthquake of the nonribbed conventional steel plain longitudinal bars (con-
damage assessment of members that are wrapped with FRP, and the ventional steel bars without deformation), and the effects of strain
high cost of epoxy resin (Triantafillou et al. 2006; Bournas et al. penetration into the foundation. Comparisons between the experi-
2009; Awani et al. 2017). To address the problems that stem mental and numerical results show satisfactory agreement for the
from epoxy resins that are utilized as matrix material, a new branch load–displacement curves, which demonstrates the acceptable per-
in strengthening methods has been burgeoned by the substitution of formance of the proposed modeling approach. Finally, the accuracy
epoxy resins with cement-based mortar. This new composite mate- of seismic assessment documents [e.g., the European Committee
rial consists of a fabric with open mesh geometry and cement-based for Standardization’s 2005 code, Eurocode 8: Design of structures
matrixes. These composite materials are known as fiber-reinforced for earthquake resistance; Parts 1–3: Strengthening and repair of
cementitious matrix (FRCM), which consists of a cementitious buildings (EC8-3) (CEN 2005), the 2018 Turkish Building Earth-
matrix that is reinforced with fibers. This relatively low-cost mate- quake Code (TBEC) (AFAD 2018), and ASCE 41-23 (ASCE
rial is compatible with concrete (Colajanni et al. 2014; Awani 2023)] when estimating the plastic rotation limits for different seis-
et al. 2017; Ortlepp and Ortlepp 2017; Ates et al. 2019; mic performance levels are evaluated.
Gonzalez-Libreros et al. 2019; Faleschini et al. 2020) and masonry
substrates (Mezrea et al. 2017, 2021; Del Zoppo et al. 2019; Napoli
and Realfonzo 2022) and could be implemented at a wider range of Experimental Program
temperatures as well as on wet surfaces. In addition, after the re-
moval of the plaster, the concrete surface does not need to be pre- In the experimental work, eight full-scale RC columns with cross-
pared before strengthening, which eliminates the need for primer sectional dimensions of 300 × 300 mm and 300 × 600 mm were
usage as filler for the surface. Few studies were conducted on the cast. In each group, one of the columns represented the columns
evaluation of the seismic behavior of the FRCM confined columns that were constructed and complied with the TSDC (Turkish Gov-
(Bournas et al. 2009; Luleci et al. 2021; Kian et al. 2022; Ates et al. ernment Ministry of Reconstruction and Resettlement 1975) (i.e.,
2023b). Raoof and Bournas (2017) showed that a cementitious the code-compliant column). The three other columns were de-
composite material could surpass FRP systems at high tempera- signed to resemble the substandard columns that were built in Tür-
tures or in a fire because of the permeability and nonflammability kiye before the 2000s and did not comply with any seismic design
provided by a cement-based matrix. codes. One of the three substandard specimens was maintained
This study aims to investigate the cyclic behavior of the substan- as-is and tested without any intervention. The other two were retro-
dard RC square and rectangular columns that are subjected to high fitted with external jacketing in two different configurations: (1)
axial loads (0.75fcbh) as well as high shear demand (up to 85% of sprayed GFRM with a glass fiber content of 5% of the mixture
the shear capacity of substandard columns at flexural strength) be- weight; and (2) three layers of basalt mesh combined with sprayed
fore and after retrofitting with a sprayed glass fiber–reinforced ma- GFRM with a glass content of 3.5% of the mixture weight. The cy-
trix (GFRM) external jacket with and without basalt mesh. This clic test results of the code-compliant column served as control col-
method saves a great deal of time due to the easy application of umns to give intuition that whether performance improvement of
the matrix material to the existing concrete surface using a spraying retrofitted columns could capture that of the code-compliant one.
gun and requires less labor compared with traditional retrofitting The specimens were designated with two letters, where the first
methods (e.g., enlarging the cross section with concrete jacketing, letters refer to code-compliant (C), substandard (S), substandard
the application of FRCM using a trowel, and adding steel retrofitted with sprayed GFRM only (G), and retrofitted with basalt
Specimen Aspect ratio Explanation Axial load ratio ACI 318-19 TBEC a/d ACI 318-19 TBEC
C1 1 TSDC-compliant 0.40 0.68 1.1 3.81 0.38 0.43
S1 1 Substandard 0.75 0.20 0.65 3.81 0.62 0.80
G1 1 GFRM-retrofitted 0.75 0.20 0.65 3.81 0.72 0.75
GB1 1 GFRM + 3layers basalt mesh 0.75 0.20 0.65 3.81 0.71 0.73
C2 2 TSDC-compliant 0.40 0.99 0.85 2.21 0.34 0.33
S2 2 Substandard 0.75 0.63 0.55 2.21 0.80 0.71
G2 2 GFRM-retrofitted 0.75 0.63 0.55 2.21 0.84 0.83
Downloaded from ascelibrary.org by "Indian Institute of Technology, Hyderabad" on 04/28/25. Copyright ASCE. For personal use only; all rights reserved.
GB2 2 GFRM + 3 layers basalt mesh 0.75 0.63 0.55 2.21 0.85 0.84
mesh–reinforced GFRM (GB). The numbers refer to the cross- in the failure. Symmetrically, 10 nonribbed plain bars with a
sectional aspect ratio (e.g., either 1 or 2). Table 1 gives the speci- 16 mm diameter that corresponded to the longitudinal reinforce-
men matrix and details of the tested columns. In addition, ment ratios of 2.2% and 1.1% (calculated according to As/bh
the ratios of the transverse reinforcement of the tested columns to where the total area of longitudinal reinforcements in the cross sec-
the minimum required transverse reinforcement (ρsh/ρsh,min) that tion is As) were used in the square and rectangular specimens, re-
correspond to ACI 318-19 (ACI 2022) and TBEC (AFAD 2018) spectively. Plain bars with 10 mm diameters and spacings of 250
are listed in this table with the ratios of Ve/Vr calculated according and 80 mm were used as the transverse reinforcement that corre-
to these documents, where the shear force that corresponds to the sponded to 0.25% and 0.78% volumetric transverse reinforcement
moment capacity of the column is Ve and the shear capacity of ratios ( ρsh ) for the substandard and code-compliant specimens, re-
the column with the contribution of the concrete and transverse re- spectively. The volumetric transverse reinforcement ratio was cal-
inforcement is Vr. As given in Table 1, the axial load ratio is the culated according to ρsh = Ash/bks where the area of the transverse
ratio of the applied axial load-to-axial capacity of the RC column. reinforcement in the direction of the loading is Ash, the distance be-
All columns were cast together with overdesigned concrete tween the centroids of the extreme stirrups leg perpendicular to the
foundation blocks with plan dimensions of 1,200 × 1,200 mm and loading direction is bk, and the spacing of the stirrups is s. In addi-
1,800 × 14,00 mm for the square and rectangular column cross sec- tion, in the tested substandard columns, the ratio of the stirrup spac-
tions, respectively. The foundations had a height of 700 mm and ing to the maximum allowable stirrup spacing based on ACI
were intentionally over-designed to ensure sufficient strength for 318-19 (ACI 2022) and TBEC (AFAD 2018) were 1.90 and
testing the behavior of the columns. This was done because sub- 2.60, respectively. The geometry, reinforcement details, and
standard columns typically fail during earthquakes before the foun- mounted strain gauge (SG) positions for the square and rectangular
dations, especially when soil properties are not a significant factor columns are shown in Figs. 1 and 2, respectively.
Fig. 1. Reinforcing details and dimensions for the square columns (mm). Note: YFLA5, FLA5 = the commercial names of the strain gauges type that
were used for longitudinal and transverse bars, respectively with a gauge length of 5 mm and water resistance.
Fig. 2. Reinforcing details and dimensions for the rectangular columns (mm).
Table 2. Concrete mix proportions lower than 14 MPa. The compressive strength of the concrete
that was used for the construction of the footings was 35.6 MPa.
Crushed aggregate Sand
The average compressive strength was calculated from the uniaxial
(5–12 mm) (0–5 mm) Cement Water
Columns (kg/m³) (kg/m³) (kg/m³) (kg/m³)
compressive test results for three standard cylinder specimens (150
× 300 mm) in accordance with ASTM C39/C39M (ASTM 2021a).
Substandard 769 1,153 180 198 The first and last RC column specimens were tested at 62 and
Code-compliant 727 1,090 257 231
106 days after the specimens were cast, which corresponds to an
average concrete age of 84 days. Three tensile tests were carried
out on the longitudinal and transverse nonribbed plain reinforcing
Concrete and Reinforcement Material Properties
bars in accordance with ASTM A615/A615M (ASTM 2022a). The
yielding and ultimate stresses of the longitudinal ( fyl and ful) and
A ready-mixed concrete with average 28 and 90-day compressive
transverse ( fyw and fuw) bars were 324 and 404 MPa, and 361 and
strengths of 9.8 and 10.5 MPa for substandard columns and 19.3
430 MPa, respectively. The yielding and ultimate strains of the lon-
and 19.5 MPa for code-compliant columns were used, respectively.
gitudinal (εyl and εul ) and transverse (εyw and εuw ) bars were
Table 2 gives the concrete mix proportions for the substandard and
0.0015 and 0.28, and 0.0017 and 0.27, respectively. The COV
code-compliant specimens. Table 3 lists the concrete properties.
for the yielding and ultimate stresses and strains of the longitudinal
The concrete elastic modulus is Ec, the strain that corresponds to
and transverse reinforcement was 8%.
the peak stress is εco , the ultimate strain is εc85 (strain that corre-
sponds to 85% of the stress in the descending branch), the average
is μ, and the standard deviation is σ. The modulus of elasticity of
the concrete was calculated according to ASTM C469/C469M Material Properties of Retrofitting System
(ASTM 2022b), which defines the chord modulus as the slope be-
tween the stress levels that correspond to axial strains of 0.00005 The ingredients of the mortar in the GFRM were white standard
and 0.4fc. portland cement (900 kg/m3), fine silica sand (1,000 kg/m3), meta-
The aim of using two types of concrete was to assess the effec- kaolin (100 kg/m3), latex modifier (33 kg/m3), water (320 kg/m3),
tiveness of the retrofitting material when improving the seismic be- and polycarboxylate ether superplasticizer (2.4 kg/m3) with a mix
havior of substandard columns and to evaluate whether retrofitted ratio of 1:1.11:0.11:0.036:0.35:0.0026 by weight. The 28-day aver-
substandard columns could achieve or exceed the performance of age compressive strength, σ, and COV of five standard 50 mm cubic
code-compliant columns. Of note, the TSDC (Turkish Government specimens in accordance with ASTM C109/C109M (ASTM 2021b)
Ministry of Reconstruction and Resettlement 1975) did not allow were 47.5, 1.3 MPa, and 2.8% for the GFRM with a 3.5% (by
for the use of concrete with a characteristic compressive strength weight) short glass fiber content, and 41.4, 0.5 MPa, and 1.4% for
layers of basalt mesh were 9.7, 0.0091, 19.9 MPa, and 27.4 GPa, (a) (b) (c)
respectively. However, for the GFRM with a 5% glass fiber con-
tent, these were 7.42, 0.0077, 19.4 MPa, and 26.9 GPa (Ates
et al. 2023a). Of note, the cost of the retrofitting material with a
GFRM is approximately 25% less than that of the retrofitting tech-
nique where a combination of GFRM and basalt mesh is utilized.
The Short Nippon electric glass fibers (Shiga, Japan) that were
used in this study contain 17% zirconia, which assures resistance
against alkali in the cement. The elastic modulus, tensile strength,
and failure strain of the glass fibers are 75 GPa, 1,500 MPa, and
0.02, respectively, as provided by the manufacturer. These excep-
tional tensile characteristics contribute to the high tensile strength
of the GFRM, which provided adequate confinement to the col-
umns in this study. In addition, compared with carbon fibers,
glass fibers are cost-effective with alkali resistance. The technical
properties of the basalt mesh reinforcement for the weft and wrap
orientations with the same number of fibers in both directions pro-
vided by the (Spinteks) manufacturer are given in Table 4. The
(d)
height of the column is covered with a total width of basalt mesh
(1,000 mm). Fig. 3. Showing: (a) rounded corners of columns; (b) GFRM spraying;
Using a concrete grinding machine, the four corners of the col- (c) wrapping of first layer of basalt mesh; and (d) retrofitted column
umns were rounded off to a radius of 25 mm before retrofitting to cross section, basalt mesh, and glass fibers (mm). Note: R25 = the ra-
prevent premature stress concentrations on the corners of the retro- dius of rounded off corner.
fitting jacket. Then, using sandpaper, the rounded corners were
smoothed [Fig. 3(a)]. Fig. 3(b) shows the GFRM retrofitting pro-
cess, and Fig. 3(c) shows the basalt mesh–reinforced GFRM appli-
cation. Fig. 3(d) shows the basalt mesh grid and the glass fibers that cross-sectional analyses that were employed in the available mod-
were used in the mortar with cross-sectional representation. Of els in the literature (Triantafillou et al. 2006; Ates et al. 2019) to
note, bidirectional PVC-coated and silane fiber-sized dry basalt achieve the targeted structural seismic performance. To ensure
mesh reinforcement with mesh spacings of 25 mm in two directions the uniform thickness of the GFRM jacket, two wooden sticks,
were used in conjunction with the GFRM with a 3.5% glass fiber each 30 mm thick, were fixed to the top and bottom of the column,
content [Fig. 3(d)]. The mesh spacing was large to permit the which served as guides to maintain a consistent jacket thickness.
proper dispersion of the GFRM in the spaces in the basalt mesh re- The thickness of the GFRM was controlled using a depth gauge.
inforcement. Retrofitting was applied at heights of 1,000 mm from To avoid damage that was induced by cracks above the externally-
the bottom for the columns, 54 days after casting, using the sprayed jacketed plastic hinge zone during the reversed cyclic lateral load-
GFRM with or without the basalt mesh. Glass fibers with average ing, which was reported by other researchers (Tong et al. 2020;
filament diameters of 18 ± 2 μm and 32 mm lengths were chopped
Ates 2022; Ates et al. 2023b), the total height of the column
automatically by the spraying gun and mixed into the mortar that
(1,000 mm) was retrofitted instead of only the potential plastic
came from a tube. Then, the GFRM was sprayed (at a distance of
hinge region. For the specimens that were retrofitted with basalt
approximately 30–50 cm from the columns) onto the surface of
mesh–reinforced GFRM (which is referred to as basalt mesh–rein-
the members with a thickness of 30 mm [Figs. 3(b and c)]. The
forced FRCM), the process began by spraying the first GFRM layer
thickness of the GFRM jacket was determined using numerical
(approximately 7 mm thick) onto the concrete substrate. Then, a
basalt mesh was embedded in the GFRM, which ensured an even
Table 4. Basalt mesh reinforcement technical properties
spacing of 7–8 mm between each layer. This process was repeated
Property Value layer by layer until the final basalt mesh was embedded, and a final
3
Filament density (kg/m ) 2,620–2,650 GFRM layer was sprayed onto the column surface. Finally, an
Specific surface weight (g/m2) 303 overlap of 200 mm was provided in the transverse direction. The
Thickness (mm) 0.8–0.9 GFRM was evenly distributed over the concrete substrate using a
Width (mm) 1,000 roller with smooth, consistent pressure. Wooden sticks were used
Tensile strength (kN/m) 95 to maintain a 30 mm gap between the bottom of the jacket and
Young’s modulus (GPa) >90 the top of the foundation, which prevented early damage from
Note: Technical properties were provided by the manufacturer (Spinteks, the jacket coming into contact with the foundation under high
Denizli, Türkiye). axial loads and lateral displacement.
(a) (b)
(a) (b)
Fig. 5. Instrumentation of: (a) square; and (b) rectangular columns (mm).
deformation. In addition, the efficiency of the confinement due to transverse reinforcement ratios for the code-compliant and substan-
the smaller cross-sectional aspect ratio of the square column was dard specimens.
considered another reason for the higher drift capacity. For S1, con- By comparing the ultimate drift capacity of the substandard
crete cover crushing and spalling occurred at 1% and 1.5% drift ra- and code-compliant specimens, the necessity of retrofitting to im-
tios, respectively. During the first cycle of the 1.5% drift ratio, the prove the seismic performance of substandard columns was noted.
lateral load capacity loss was 24.2%, and an axial load capacity For the retrofitted Columns G1 and GB1, a vertical crack in the
loss of 47.6% was recorded during the second cycle of the 1.5% external jacket was observed on the edge of the external jacket
drift ratio [Fig. 6(c)]. For S2, concrete cover crushing and spalling for the 2% drift ratio; however, for G2 and GB2, it was observed
occurred at the 0.75% drift ratio. During the first cycle of the on the edge of the jacket for the 0.75% and 1% drift ratios, respec-
0.75% drift ratio, the lateral load capacity loss was 20.2%. Here, tively. The cracks were initiated from the edge of the jacket due to
the axial load decreased by 14%. Therefore, the failure of the sub- stress concentration in these regions. (Ates et al. 2019). Then, the
standard specimens was brittle due to spalling of the concrete crack width was widened, and the crack length was extended in
cover, buckling of the longitudinal reinforcements, and crushing the plastic hinge zone from the column–foundation interface up
of the core concrete sequentially. Of note, core concrete crushing to 300 mm for G1 and GB1 (Fig. 7) and 600 mm for G2 and
was detected after the cover concrete was removed following the GB2 (Fig. 8). This is due to the different plastic hinge zone
completion of testing [Figs. 7 and 8]. In addition, S2 experienced lengths of 300 and 600 mm for square and rectangular columns.
more brittle failure due to more shear-dominated behavior that The failure of G1 and G2 was captured by the increase in crack
was triggered by a lower a/d ratio compared with S1. Comparing width and rupture of the GFRM jacket. The failure modes of
the results between S1 and C1 revealed that cover crushing and GB1 and GB2 were accompanied by basalt mesh rupture at 3%
spalling in S1 occurred three times earlier than C1 in terms of and 1.5% drift ratios after widening the formed crack on the jacket
drift ratio, while failure of S1 occurred 2.66 times earlier than C1 at 2% and 0.75% drift ratios, respectively [Figs. 6(g and h), 7, and
in terms of drift ratio. The same values were obtained for the C2 8]. The basalt mesh rupture could be justified because the confin-
and S2 specimens for crushing and spalling. All these behavioral ing stress exceeded the tensile strength of the jacket. Then, by ex-
differences stem from different applied axial load-to-capacity and tending and widening the crack, the tensile stress transfers
(a) (b)
(c) (d)
(e) (f )
(g) (h)
Fig. 6. Hysteresis curves of specimens: (a) C1; (b) C2; (c) S1; (d) S2; (e) G1; (f) G2; (g) GB1; and (h) GB2.
through the cracks to the basalt mesh, which results in basalt rup- was promising to overtake and reach the target lateral load capac-
ture. Compared with S1 and S2, a more gradual failure was wit- ity of the code-compliant specimen (C1).
nessed for the retrofitted specimens, similar to C1 and C2. The As shown in Figs. 7 and 8, substandard columns experienced
lateral load capacity of S1 was 26% less than C1. When compared buckling in the longitudinal rebars, which led to early concrete
with S1, the lateral load capacities of G1 and GB1 were increased cover spalling and a brittle failure mode with strength degradation.
by 35% and 32%, respectively. The lateral load capacity of code- This failure occurred due to the high axial load level, which re-
compliant Specimen C1 was successfully achieved by retrofitting sulted in a loss of axial stiffness and led to the loss of the lateral
the square columns. The lateral load capacity of S2 was 27% less stability of the longitudinal reinforcing bar. This instability was ex-
than C2. Compared with S2, the lateral load capacities of the G2 acerbated by the weak lateral confining pressure and the high un-
and GB2 specimens were 15% and 16% higher, respectively. supported length of the longitudinal bar [(s/db) where the
Therefore, in square specimens, the applied retrofitting method diameter of the longitudinal bar is db]. However, in the
in
in
in
GFRM-retrofitted columns, although the bars reached their com- reported by Ates et al. (2023b). The results of this study indicate
pressive yield strain and buckled under severe axial load, they that, in addition to the s/db ratio, the onset of bar buckling depends
could sustain axial load at larger drift ratios. This was possible be- on the confinement that is provided by the external jacket (e.g., G1,
cause the external confinement restricted the cover spalling and ex- GB1, G2, and GB2) and internal transverse reinforcement (e.g., C1
erted inward lateral pressure on the concrete, which prevented and C2). A summary of the test results for the square and rectangu-
sudden bar buckling by restricting the lateral deformation, as lar specimens is given in Table 5.
ternal jacket fractured in GB1 and GB2, the basalt meshes that bore
the tensile stresses provided confinement to the core concrete and
prevented crack extension and widening in contrast to G1 and
G2. In addition, the ultimate tensile strain of the GFRM with
3.5% glass content and three layers of basalt (GB specimens)
mesh was 18% higher than that of the GFRM with 5% glass content
(G specimens) (Ates et al. 2023a). The ductility increase was 29%
less for the G2 with respect to G1. For GB2, the ductility increase
was 23.4% less compared with the GB1 specimens. This was be-
cause of a lower a/d ratio in the rectangular columns (42%),
Fig. 9. Yield and ultimate displacements, dissipated energy, and secant which decreased the confining efficiency due to predominant
stiffness definition. Note: Pmin = maximum lateral force in pulling di- shear deformations. In addition, the confining efficiency that
rection; Pmax = maximum lateral force in pushing direction; and Pi = stemmed from the retrofitting on the rectangular columns was
the lateral force corresponding to target drift at i-th cycle. less than the square columns, as reported by Demir et al. (2019).
(a) (b)
(c) (d)
(e) (f )
(g) (h)
Fig. 10. Strain profiles: (a) C1; (b) C2; (c) S1; (d) S2; (e) G1; (f) G2; (g) GB1; and (h) GB2. The horizontal solid line refers to the column–foundation
interface; C = compression; and T = tension.
there is a lower strain distribution along the longitudinal reinforce- Specimens C1 and C2 reached compression yield strain values
ment, indicating that bar slippage inside the column above h was at 3% and 2% drift ratios, respectively [Figs. 10(a and b)]. Yielding
insignificant. This agreed well with the damage pictures (Figs. 7 and postyielding compression strains for S1 and S2 occurred at the
and 8), where above the plastic hinge zone, there was no splitting sections 150 and 300 mm above the footing at 1% and 0.5% drift
or cracking that originated from the bar slippage. ratios, respectively, which was consistent with the observed
and h)]. These results indicate that a GFRM jacket could drastically lated using 1.2VeLV/AG, where the shear modulus is G and was cal-
restrict the compressive damage to longitudinal bars. culated using Ec/2(1 + νc) (Pujol 2002). Here, the distance from the
It was assumed that the vertical displacement that was measured actuator to the middle of the first gauge length is Lv, the cross-
by the first LVDT from the surface of the foundation with a gauge sectional area is A, and Poisson’s ratio is νc. The elastic shear dis-
length of 30 mm was composed of slip in the longitudinal bars placements for the S1 and S2 specimens were 0.18 and 0.24 mm,
(Rodrigues et al. 2013; Seifi et al. 2018). The vertical displacement respectively.
that was measured by two subsequent LVDTs with gauge lengths The results revealed that a significant proportion of the top dis-
of 120 and 150 mm were assumed to be composed of the flexural placement (24%–51%) originated from the column–foundation
(a) (b)
(c) (d)
(e) (f )
(g) (h)
Fig. 11. Contribution of lateral displacement components on top displacement; (a) C1; (b) C2; (c) S1; (d) S2; (e) G1; (f) G2; (g) GB1; and (h) GB2.
Stiffness Degradation
The average peak-to-peak secant stiffness (Ki) versus drift ratio of
the second cycle in the push and pull directions for the square and (b)
rectangular specimens are shown in Fig. 12. The secant stiffness
Fig. 13. Residual displacement of (a) square; and (b) rectangular spec-
is defined as the ratio of the lateral load to the corresponding dis-
imens per drift ratio.
placement at each cycle (Fig. 9). Rodrigues et al. (2013) showed
(a)
(b) (c)
Fig. 14. (a) Modeling details that consider the slip proportion in the top displacement; (b) uniaxial Material Pinching 4 model; and (c) defined stress–
strain curve for concrete and reinforcement. Note: Sec. 7 = section 7 at which 7th integration point was assigned; and N, M = normal applied axial
load and bending moment.
rotations in the column–foundation interface that were induced by stress–strain curve of concrete until the stress level of 0.2fc
strain penetration. Then, using the uniaxialMaterial Pinching4
model, a moment–rotation spring was defined and associated Gfc 0.8fc
ε20 = − + εco (4)
with the zero-length element. The parameters for the moment–rota- 0.6fc LIP Ec
tion spring model of the specimens were determined using several
steps. First, theoretical tensile force values for the longitudinal bars where LIP = length of the softening integration point. Of note, the
were derived from the moment–curvature analysis for various mo- fracture energy in tension is defined as the required energy to be ab-
ment values. Next, to determine the longitudinal rebar slip that cor- sorbed to create a unit area of crack (Oshtolagh et al. 2023). How-
responded to each tensile force, the so-called good bond stress (τb)– ever, it is accepted for localization in compression (Coleman and
slip curve given in the Fib Model Code (fib 2010) was selected Spacone 2001). In Eq. (4), the compression fracture energy (Gfc )
using a trial-and-error process. Then, this curve was applied to is defined as the area under the postpeak stress–strain curve of con-
19 joint link elements that are associated with a single bar in crete compressive strength, which is from the peak stress (fc ) until
SAP 2000 (CSI 2010). A nonlinear static displacement-controlled the end of the softening branch (0.2fc). Therefore, by calibrating the
load case was defined to the bar’s end node, which was in the col- ε20 (strain that corresponds to 0.2fc in the descending branch), a
umn–foundation interface to implement the tensile analysis by pull- constant energy release was ensured. Therefore, the numerical
ing it out to obtain the appropriate force–slip curve. The other end model’s postpeak behavior improved, as shown in Fig. 15, and a
of the bar was fixed using constraints in all degrees of freedom due good agreement was achieved with the experimental result.
to sufficient anchorage that was provided by enough development As shown in Fig. 16(a), the inclusion of buckling in bars refined
length. Then, by dividing the obtained bar slip by the neutral axis the postpeak response of the substandard columns. In addition,
depth, the rotation of the column due to the longitudinal bar slip Figs. 16(b and c) show the example recorded stress–strain curve
was calculated. Finally, with the rotation and corresponding mo- of longitudinal reinforcement for the S1 and G1 specimens. Ac-
ment values at each step, the moment–rotation curve for the rota- cording to the experiments, the reinforcing bar could reach a
tional spring was obtained. yield plateau in the retrofitted and code-compliant columns; how-
ever, they buckled under severe axial stresses in the substandard
columns. Figs. 16(b and c) show that they agreed well with the ob-
Employed Regularization Method and Numerical served behavior in the tests. As shown in Fig. 16(b), the bar buck-
Refinements led severely after reaching its peak stress. However, as shown in
Fig. 16(c), the bar did not buckle until larger drift ratios and strains
Strain softening due to localization in force-based elements during were experienced.
the postpeak response occurred in S1. Therefore, the postpeak Of note, the modeling of substandard RC columns that have ma-
load–displacement of the numerical model had a steeper descend- terial properties out of the engineering range (e.g., a compressive
ing branch than the experimental result. Coleman and Spacone concrete strength of less than 12 MPa) is extraordinarily cumber-
(2001) proposed two methods of regularization to obtain the objec- some. Therefore, capturing all occurrences that were witnessed dur-
tive response in the postpeak zone: (1) constant fracture energy cri- ing the test was challenging. Considering buckling of the
terion on the element force–displacement level; and (2) geometric longitudinal reinforcements, solving the strain and curvature local-
scaling on the local moment–curvature level. Of note, the number ization at the first integration point as well as accounting for the
and location of the integration points affect the overall section cur- longitudinal rebar slip contribution to the top displacement are
vature demand and element response. In addition, they provided an some superiorities (e.g., ability to capture key damage mechanisms
equation to control the linear postpeak softening branch of the or failure modes, strain penetration, and reinforcement buckling in
Specimen θd θ DL θ SD θ NC PL θ SH θ KH θ GO PL
S1 0.012 0.0038 0.008 0.011 NC 0 0.003 0.004 GO
C1 0.033 0.0033 0.017 0.022 NC 0 0.009 0.012 GO
G1 0.024 0.0039 0.008 0.011 NC 0 0.004 0.006 GO
GB1 0.028 0.0039 0.008 0.011 NC 0 0.004 0.006 GO
S2 0.004 0.0032 0.007 0.009 DL 0 0.002 0.003 GO
C2 0.021 0.0032 0.011 0.015 NC 0 0.008 0.010 GO
G2 0.016 0.0033 0.007 0.009 NC 0 0.003 0.005 GO
Downloaded from ascelibrary.org by "Indian Institute of Technology, Hyderabad" on 04/28/25. Copyright ASCE. For personal use only; all rights reserved.
oval profile to promote uniform stress distribution and prevent lo- displacement that stemmed from rotational slip (bond degrada-
calized failures. tion) at the column–foundation interface due to strain penetra-
tion (24%–51%).
3. The achieved displacement ductilities for the retrofitted speci-
mens were in the range of μD = 4.2–5.4 for the G and GB
Conclusions specimens. The ductility of G1 and GB1 increased up to 80%
compared with the substandard Specimen S1, which ap-
Eight full-scale RC columns with varying cross-sectional aspect
proached the 94% increase in ductility that was exhibited by
and a/d ratios were tested under a constant axial load and reversed
C1 with respect to S1. Similarly, the ductility of the
cyclic lateral displacements. A novel retrofitting method was im- GFRM-retrofitted rectangular columns increased up to 62%
plemented using a GFRM with a 5% glass content as well as a compared with their substandard counterpart S2, with a lower
GFRM with a 3.5% glass content that was reinforced with three a/d and higher shear demand. Of note, GB1 and GB2 showed
layers of basalt mesh. This retrofitting significantly enhanced the slightly higher ductility increments of 30% and 27% compared
performance of substandard columns under high axial load ratios with G1 and G2, respectively.
(0.75) and shear demand (Ve/Vr of 0.71 and 0.80), which brought 4. The numerical model was validated accurately against experi-
their performance closer to that of code-compliant columns. The mental results that reproduced the hysteretic lateral load–
experimental results were discussed comparatively, which focused displacement curves and incorporated the main failure mecha-
on key behavioral metrics, such as strength, stiffness, energy dissi- nisms that were observed in the tests, such as longitudinal bar
pation, and observed failure mechanisms. An accurate numerical buckling and slip due to strain penetration at the column–foot-
model that incorporates the nonlinear effects, which includes strain ing interface.
penetration and buckling of longitudinal rebars in OpenSees, was 5 The rotational plastic demand of the retrofitted specimens in-
developed, and global hysteretic responses (e.g., ultimate displace- creased compared with the substandard ones. In addition, the
ment and lateral load capacity) were validated against the experi- success of the codes when evaluating the rotational plastic de-
mental data. The success of EC8 (CEN 2005) and TBEC (AFAD mand of the specimens was scrutinized across the three col-
2018) in estimating the performance limits of columns was as- umns. The results indicated that the plastic rotation demand of
sessed. Based on the experimental and modeling results, the fol- substandard columns could be properly predicted by TBEC
lowing major conclusions are drawn. (AFAD 2018), although EC8-3 (CEN 2005) overestimated the
1. The S1 and S2 columns showed a brittle failure mode that was plastic rotation demand by 60% compared with the test results
characterized by cover spalling, concrete crushing, and longitu- of S2. This discrepancy was attributed to the high axial
dinal bar buckling, which highlights the necessity of retrofitting load-to-capacity ratio of 0.75 in the substandard columns,
to mitigate the lack of confinement. Retrofitted square columns which requires more careful consideration in EC8-3 (CEN
(G1 and GB1) with a higher a/d ratio of 3.81 and lower shear 2005). In addition, ASCE 41-23 considers no plastic rotation ca-
demand (Ve/Vr of 0.72) exhibited a 35% increase in lateral pacity for S1 and S2, despite the fact that these columns demon-
load bearing capacity compared with S1 (Ve/Vr of 0.62). The lat- strated a limited but observable plastic rotation capacity.
eral load bearing capacity of C1 (Ve/Vr of 0.38) was 32% higher
than that of S1. For rectangular retrofitted columns (G2 and
GB2) with a lower a/d of 2.21 and higher shear demand (Ve/
Vr of 0.85), the lateral load capacity increased by 16% compared Data Availability Statement
with S2 (Ve/Vr of 0.80); however, C2 (Ve/Vr of 0.34) had a lat-
eral load capacity 36% higher than that of S2. The greater effi- All data that support the findings of this study are available from the
ciency of retrofitting in square columns was attributed to the corresponding author upon reasonable request.
20% lower shear demand due to the higher a/d ratio and the ad-
verse relationship between the confinement efficiency and
cross-sectional aspect ratio. Acknowledgments
2. Despite the relatively high shear demand in substandard col-
umns, the top displacement was primarily governed by flexural The authors are most thankful to the Fibrobeton company and its
displacement, along with significant contributions from board member, Muhammed Marasli, for their continuous support
Notation
μ = average;
The following symbols are used in this paper: μD = displacement ductility;
A = cross-sectional area; νc = concrete Poisson’s ratio;
As = total area of longitudinal rebars; ρsh = volumetric transverse reinforcement ratio;
Ash = area of the transverse reinforcement alongside the ρsh,min = minimum required volumetric transverse
lateral load; reinforcement ratio;
a = shear span; σ = standard deviation;
b = column dimension (width); τb = bond stress; and
bk = distance between the centroid of the extreme stirrups φaverage = average curvature.
leg;
d = effective depth of cross section;
d0 = distance from the first LVDT to bottom of the column;
db = diameter of longitudinal bar; References
dl = measured left displacement;
dr = measured right displacement; ACI (American Concrete Institute). 2013. Guide for testing reinforced con-
Ec = elastic modulus of concrete; crete structural elements under slowly applied simulated seismic loads.
Es = the slope of compression stress-strain curve of ACI 374.1. Farmington Hills, MI: ACI.
reinforcement in descending branch after buckling; ACI (American Concrete Institute). 2019. Building code requirements for
Ets = the slope of tensile stress-strain curve of concrete in structural concrete. ACI 318. Farmington Hills, MI: ACI.
ACI (American Concrete Institute). 2020. Guide to design and construction
descending branch;
of externally bonded fabric-reinforced cementitious matrix and steel-
fc = concrete compressive strength;
reinforced grout systems for repair and strengthening concrete struc-
ful = ultimate stress of longitudinal rebar; tures. ACI 549. Farmington Hills, MI: ACI.
fuw = ultimate stress of transverse rebar; ACI (American Concrete Institute). 2022. Seismic evaluation and retrofit of
fyl = yield stress of longitudinal rebar; existing concrete buildings—Code and commentary. ACI 369.
fyw = yield stress of transverse rebar; Farmington Hills, MI: ACI.
G = shear modulus; AFAD (Ministry of Interior Disaster and Emergency Management
Gfc = fracture energy in compression; Authority). 2018. Turkish building earthquake code. Ankara,
h = column dimension (height); Türkiye: AFAD.
Ki = peak-to-peak ith cycle secant stiffness; Ali, M. A., and E. El-Salakawy. 2016. “Seismic performance of
LG = height of the measurement segment; GFRP-reinforced concrete rectangular columns.” J. Compos. Constr.
Lh = horizontal distance between two LVDTs; 20 (3): 04015074. https://doi.org/10.1061/(ASCE)CC.1943-5614
LIP = length of softening integration point; .0000637.
Ansal, A., A. Akinci, G. Cultrera, M. Erdik, V. Pessina, G. Tönük, and G.
Lu = unsupported length of longitudinal bar in height
Ameri. 2009. “Loss estimation in Istanbul based on deterministic earth-
(stirrup spacing);
quake scenarios of the Marmara Sea region (Turkey).” Soil Dyn.
Lv = distance from actuator to middle of the first gauge Earthquake Eng. 29 (4): 699–709. https://doi.org/10.1016/j.soildyn
length; .2008.07.006.
l = clear height of the column to horizontal load (shear ASCE. 2023. Seismic evaluation and retrofit of existing buildings.
span); ASCE41-23. Reston, VA: ASCE.
lp = plastic hinge length; ASTM. 2021a. Standard test method for compressive strength of cylindri-
Mw = moment magnitude; cal concrete specimens. ASTM C39/C39M-21. West Conshohocken,
s = stirrup spacing; PA: ASTM.
Ve = shear force that corresponds to the bending moment; ASTM. 2021b. Standard test method for compressive strength of hydraulic
Vr = shear capacity of cross section; cement mortars (using 2-in. or [50-mm] cube specimens). ASTM
Δslip = contribution of slip displacement component; C109-21. West Conshohocken, PA: ASTM.
Δt = total displacement; ASTM. 2022a. Standard specification for deformed and plain carbon-steel
bars for concrete reinforcement. ASTM A615/A615M-22. West
Δflexure = The contribution of flexural displacement;
Conshohocken, PA: ASTM.
δi = target displacement of ith cycle; ASTM. 2022b. Standard test method for static modulus of elasticity and
δres = residual lateral displacement at zero lateral force after Poisson’s ratio of concrete in compression. ASTM C469/C469M.
unloading; West Conshohocken, PA: ASTM.
δ+u = ultimate displacement in pushing; Ates, A. O. 2022. “Improvement of seismic performance of reinforced con-
δ−u = ultimate displacement in pulling; crete columns using glass fiber reinforced sprayed mortar with/without
δ+y = yield displacement in pushing; textile reinforcement.” Ph.D. thesis, Dept. of Civil Engineering,
δ−y = yield displacement in pulling; Istanbul Technical Univ.
square columns strengthened with self-compacting concrete-filled State-of-the-art review.” Sustainability 11 (11): 3208. https://doi.org/10
CFRP-steel tubes.” J. Bridge Eng. 24 (2): 04018119. https://doi.org/ .3390/su11113208.
10.1061/(ASCE)BE.1943-5592.0001345. Rodrigues, H., A. Arêde, H. Varum, and A. G. Costa. 2013. “Experimental
Luleci, M. C., B. Sari, U. Demir, M. Marasli, and A. Ilki. 2021. “Seismic be- evaluation of rectangular reinforced concrete column behaviour under
havior of substandard extended rectangular RC columns jacketed with biaxial cyclic loading.” Earthquake Eng. Struct. Dyn. 42 (2): 239–
sprayed GFRM.” In COMPDYN Proc., 8th Int. Conf. on Computational 259. https://doi.org/10.1002/EQE.2205.
Methods in Structural Dynamics and Earthquake Engineering, 111– Saatcioglu, M., and S. R. Razvi. 1992. “Strength and ductility of confined
119, Athens, Greece: ECCOMAS Proceedia. concrete.” J. Struct. Eng. 118 (6): 1590–1607. https://doi.org/10.1061/
Mazzoni, S., F. McKenna, M. H. Scott, and G. L. Fenves. 2006. “Opensees (ASCE)0733-9445(1992)118:6(1590).
command language manual.” Pac. Earthquake Eng. Res. Center Seifi, A., A. Hosseini, M. S. Marefat, and M. Khanmohammadi. 2018.
264 (1): 137–158. “Seismic retrofitting of old-type RC columns with different lap splices
Mezrea, P. E., M. Ispir, I. A. Balci, I. E. Bal, and A. Ilki. 2021. “Diagonal by NSM GFRP and steel bars.” Struct. Des. Tall Spec. Build. 27 (2):
tensile tests on historical brick masonry wallets strengthened with fabric e1413. https://doi.org/10.1002/TAL.1413.
Sezen, H., M. Asce, and J. P. Moehle. 2004. “Shear strength model
reinforced cementitious mortar.” Structures 33: 935–946. https://doi
for lightly reinforced concrete columns.” J. Struct. Eng. 130 (11):
.org/10.1016/J.ISTRUC.2021.04.076.
1692–1703. https://doi.org/10.1061/(ASCE)0733-9445(2004)130:11
Mezrea, P. E., I. A. Yilmaz, M. Ispir, E. Binbir, I. E. Bal, and A. Ilki. 2017.
(1692).
“External jacketing of unreinforced historical masonry piers with open-
Tapan, M., M. Comert, C. Demir, Y. Sayan, K. Orakcal, and A. Ilki. 2013.
grid basalt-reinforced mortar.” J. Compos. Constr. 21 (3): 04016110.
“Failures of structures during the October 23, 2011 Tabanlı (Van) and No-
https://doi.org/10.1061/(ASCE)CC.1943-5614.0000770.
vember 9, 2011 Edremit (Van) earthquakes in Turkey.” Eng. Fail. Anal.
Mirmiran, A., M. Shahawy, M. Samaan, H. E. Echary, J. C. Mastrapa, and
34: 606–628. https://doi.org/10.1016/J.ENGFAILANAL.2013.02.013.
O. Pico. 1998. “Effect of column parameters on FRP-confined con-
Tong, T., H. Lei, S. Yuan, and Z. Liu. 2020. “Experimental investigation
crete.” J. Compos. Constr. 2 (4): 175–185. https://doi.org/10.1061/
and seismic vulnerability assessment of low flexural strength rectan-
(ASCE)1090-0268(1998)2:4(175). gular bridge piers retrofitted with ultrahigh-performance concrete
Napoli, A., and R. Realfonzo. 2022. “Compressive behavior of masonry jackets.” Eng. Struct. 206: 110132. https://doi.org/10.1016/J
columns confined with FRCM systems: Research overview and analyt- .ENGSTRUCT.2019.110132.
ical proposals.” J. Compos. Constr. 26 (3): 04022019. https://doi.org/10 Tore, E., C. Demir, M. Comert, and A. Ilki. 2021. “Seismic collapse perfor-
.1061/(ASCE)CC.1943-5614.0001200. mance of a full-scale concrete building with lightly reinforced col-
Narlitepe, F., N. Kian, U. Demir, C. Demir, and A. Ilki. 2021. “Seismic umns.” J. Struct. Eng. 147 (12): 04021207. https://doi.org/10.1061/
strengthening of sub-standard rc columns through a novel hybrid thin (ASCE)ST.1943-541X.0003178.
jacketing method.” In Proc., 18th fib Symp. Concrete Structures: New Triantafillou, T. C., C. G. Papanicolaou, P. Zissimopoulos, and T.
Trends for Eco-Efficiency and Performance, 1223–1232, Lausanne, Laourdekis. 2006. “Concrete confinement with textile-reinforced mor-
Switzerland: International Federation for Structural Concrete (fib). tar jackets.” ACI Struct. J. 103 (1): 28–37. https://doi.org/10.14359/
Ortlepp, R., and S. Ortlepp. 2017. “Textile reinforced concrete for strength- 15083.
ening of RC columns: A contribution to resource conservation through TSDC (Turkish Government Ministry of Reconstruction and
the preservation of structures.” Constr. Build. Mater. 132: 150–160. Resettlement). 1975. Specifications for structures to be built in disaster
https://doi.org/10.1016/J.CONBUILDMAT.2016.11.133. areas. Turkish Seismic Design Code. Ankara, Turkey: Turkish
Oshtolagh, M. R., M. Farzam, N. Kian, and H. Sadaghian. 2023. “Behavior Government Ministry of Reconstruction and Resettlement,
of recycled steel fiber-reinforced concrete beams in torsion- experimental Earthquake Research Institute.
and numerical approaches.” Comput. Concr. 32 (2): 173–184. https://doi Verderame, G. M., G. Fabbrocino, and G. Manfredi. 2008. “Seismic re-
.org/10.12989/cac.2023.32.2.173. sponse of r.c. columns with smooth reinforcement. part II: Cyclic
Ozcan, O., B. Binici, and G. Ozcebe. 2010. “Seismic strengthening of rect- tests.” Eng. Struct. 30 (9): 2289–2300. https://doi.org/10.1016/J
angular reinforced concrete columns using fiber reinforced polymers.” .ENGSTRUCT.2008.01.024.
Eng. Struct. 32 (4): 964–973. https://doi.org/10.1016/j.engstruct.2009 Zhu, L., K. J. Elwood, and T. Haukaas. 2007. “Classification and seismic
.12.021. safety evaluation of existing reinforced concrete columns.” J. Struct.
Paulay, T., and M. J. N. Priestley. 1992. Seismic design of reinforced con- Eng. 133 (9): 1316–1330. https://doi.org/10.1061/(ASCE)0733
crete and masonry buildings, 768. New York: Wiley. -9445(2007)133:9(1316).