ACI Structural Journal Sept. - 2023 V. 120 No. 5
ACI Structural Journal Sept. - 2023 V. 120 No. 5
ACI Structural Journal Sept. - 2023 V. 120 No. 5
5
SEPTEMBER 2023
ACI
STRUCTURAL J O U R N A L
@seismicisolation
@seismicisolation
CONTENTS
STATEMENT OF OWNERSHIP
Editorial Board
Michael Kreger, Editor-in-Chief
ACI Structural Journal
University of Alabama
Catherine French September 2023, V. 120, No. 5
University of Minnesota
Mary Beth Hueste a journal of the american concrete institute
Texas A&M University an international technical society
David Sanders
Iowa State University 3 Ultra-High-Toughness Concrete Retrofitted Boundary Column Shear
Gustavo Parra-Montesinos
University of Wisconsin–Madison
Walls: Tests and Capacity Prediction, by Hongmei Zhang, Giorgio
Monti, Yuanfeng Duan, Zhiyuan Chen, and Chen Gu
Board of Direction
19 Hysteretic Model of Coupler Box Assembly for Seismic Retrofitting of
President
Antonio Nanni
Severely Damaged Reinforced Concrete Buildings, by Naveen Kumar
Kothapalli, R. Siva Chidambaram, and Pankaj Agarwal
Vice Presidents
Maria Juenger 37 Seismic Performance of Five-Spiral Concrete Columns with Large-
Michael J. Paul Diameter Longitudinal Reinforcement, by Yu-Chen Ou and Brandon Li
Directors
Oscar R. Antommattei
49 Rehabilitation of Exterior Beam-Column Joint by Geopolymer Mortar
Peter Barlow under Quasi-Static Loading, by Arshad Hussain Choudhury and Aminul
Michael C. Brown Islam Laskar
Arturo Gaytan Covarrubias
Anthony R. DeCarlo Jr. 63 Strength and Behavior of Glass Fiber-Reinforced Polymer-Rein-
John W. Gajda forced Concrete Box Girders without Web Reinforcement under Pure
Carol Hayek Torsion, by Ibrahim T. Mostafa, Salaheldin Mousa, Hamdy M. Mohamed,
Kamal H. Khayat
Robert C. Lewis and Brahim Benmokrane
Anton K. Schindler
Matthew R. Sherman 77 Experimental Study on Drilled Shaft Reinforcing Bar Anchorage in
Lawrence L. Sutter Footings, by Yousun Yi, Hyunsu Kim, Jongkwon Choi, Juan Murcia-Delso,
and Oguzhan Bayrak
Past President Board Members
Jeffrey W. Coleman
Cary S. Kopczynski 93 Design of Glass Fiber-Reinforced Polymer-Reinforced Concrete
Charles K. Nmai Columns per ACI CODE-440.11-22, by Zahid Hussain and Antonio Nanni
Executive Vice President 109 Cover Spalling in Reinforced Concrete Beams Subjected to Pure
Frederick H. Grubbe Torsion, by Allan Kuan, Evan C. Bentz, and Michael P. Collins
@seismicisolation
@seismicisolation
The rapid growth of the construction industry in Asia and the conse- increases the initial dimensions (Woods et al. 2016); besides,
quent updating of design specifications put forward higher perfor- it presents the problem of transferring the bending action
mance requirements for structural components, which results in a from the existing to the new elements (Woods et al. 2016).
large number of existing shear walls that are not compliant with the Fiber-reinforced concrete (FRC) exhibited excellent
current seismic standards. A prospective retrofitting method, which
strength, toughness, and damage tolerance (Li et al. 2001;
is based on replacing the existing boundary concrete or attaching
Boshoff et al. 2009) due to crack prevention and strain-hard-
external boundary columns to nonconforming shear walls, is
experimentally studied. Four shear-wall specimens were designed ening effects. Fischer and Li (2002a) emphasized that the
according to the current Chinese design code: one using plain most important impact of high-performance FRC compos-
concrete boundary columns and three using ultra-high-toughness ites on the structural response is to maintain the integrity of
boundary columns (UHTBCs), adopting three different strength- and provide lateral stability to steel reinforcements, thereby
ening strategies relevant to the boundary size and the connection ensuring adequate cyclic inelastic deformation and energy
form. Cyclic performance, damage patterns due to UHTBCs, and dissipation. Experiments involving columns (Mirmiran
connection form are discussed based on the experimental results, and Shahawy 1997; Fischer and Li 2002a, 2003), beams
from which it was ascertained that shear walls with UHTBCs show (Li and Xu 2009; Prota et al. 2004), and beam-column
improved seismic performance, compatible with the requirements joints (Parra-Montesinos 2005; Prota et al. 2004) have
of the current seismic design code, even for the reduced-boundary
proved that structural elements with FRC exhibit improved
UHTBCs and non-connection specimens. The predictive equation
seismic resistance. The bridging effect helps the structural
for the sectional moment capacity of shear walls with UHTBCs
was discussed as a practical tool for retrofitting applications. This elements maintain the lateral bearing capacity (Naaman and
study highlights the most important features of a rapid retrofitting Najm 1991) and alleviate damage evolution (Fischer and
measure to improve the resilience of existing nonconforming shear- Li 2002b). Thanks to the high tensile-strength-to-weight
wall structures, while also proving to be an effective measure for ratio and durability, FRC is ideal for repairing and retro-
newly constructed structures. fitting infrastructures and historical buildings (Mustafaraj
et al. 2020). Other authors (Xu et al. 2012; Li et al. 2009)
Keywords: boundary columns; cyclic loading test; retrofitting; seismic proposed ultra-high-toughness cementitious composites
performance; shear wall; ultra-high-toughness.
(UHTCCs) containing high-strength short fibers, which
showed excellent performance in bending elements and
INTRODUCTION
proved to be a viable alternative to plain concrete, thanks to
Background
their strain-hardening behavior and simple applicability to
Reinforced concrete (RC) shear walls are widely used in
existing structural elements.
mid- and high-rise structures as seismic-resistant elements.
Applying FRC composites to improve the flexural
Under the combined action of flexure, shear, and torsion, the
capacity of shear walls has recently attracted research
damage tends to concentrate in the outer boundary zones of
interest (Parra-Montesinos 2005; Parra-Montesinos et al.
RC shear walls (Mutō 1969; Zhang et al. 2014). Cracked
2006; Dazio et al. 2008; Olsen and Billington 2011; Li et al.
concrete weakens its compression capacity, which in turn
2014), where the positive effects on the seismic performance
expedites concrete crushing, leading to the failure of the
were demonstrated. FRC composites are also considered to
shear walls (Fischer and Li 2002a; Zhang et al. 2014). This
apply in shear-wall retrofitting (Kesner and Billington 2005;
shear and flexure interaction combined with local damage
Liang et al. 2013; Dang et al. 2014; Lu et al. 2018; Zhang
makes retrofitting of such shear walls a challenging endeavor.
et al. 2020).
To improve the flexural and shear performance of existing
and newly constructed shear walls, different retrofit strate-
RESEARCH SIGNIFICANCE
gies have been conceived—for example, the insertion of RC
In most of the existing shear walls constructed before 2001,
boundary columns (Al-Gemeel and Zhuge 2019), attaching
due to material deterioration and old specification-based
a steel profile (Chrysanidis and Tegos 2020), or even
enlarging the boundary elements. Attaching a steel profile ACI Structural Journal, V. 120, No. 5, September 2023.
to an existing concrete structure involves serious prob- MS No. S-2021-064.R5, doi: 10.14359/51738829, received July 15, 2022, and
reviewed under Institute publication policies. Copyright © 2023, American Concrete
lems of effective collaboration and bond durability (Zhou Institute. All rights reserved, including the making of copies unless permission is
obtained from the copyright proprietors. Pertinent discussion including author’s
et al. 2010). Inserting additional jacketed concrete columns closure, if any, will be published ten months from this journal’s date if the discussion
is received within four months of the paper’s print publication.
@seismicisolation
@seismicisolation
designs, the bending capacity may not conform to current lengths of SW-1 and UCSW-2 were 200 mm, and those of
specifications. Besides, the application of strengthening UCSW-3 and UCSW-4 were 150 mm. The experiments were
measures is often impeded by limited available space. In conducted in the State Key Laboratory of Disaster Reduction
this study, an alternative seismic strengthening method in Civil Engineering at Tongji University, Shanghai, China.
using ultra-high-toughness boundary columns (UHTBCs) The specimen configurations and details are shown in
on noncompliant shear walls with weak moment capacity Fig. 1. The boundary columns of Specimens UCSW-2
is presented. With this scope, four shear walls—one of plain to UCSW-4 were cast with UHTCC. The RC core wall,
concrete representing the current seismic design code (as a top beam, and foundation were cast with plain concrete.
reference specimen), two UHTBC shear walls with equal Except for the boundary column type, UCSW-2 was iden-
and reduced boundary elements, respectively (as strength- tical to SW-1, which was tested to determine the effect of
ened specimens by inserting additional elements or replacing the UHTBCs on the seismic performance. UCSW-3 was
the boundary concrete), and one UHTBC shear wall without designed with reduced boundary column widths and longitu-
shear connections (representing a case of rapid retrofit- dinal bars to investigate the possibility of reduced boundary
ting)—were cyclically tested to validate the corresponding columns. For Specimens UCSW-2 and UCSW-3, there were
strengthening strategies. key slots and shear reinforcing bars for connection to the
This study provides in-depth considerations on this rapid inner core wall. Half-grouted sleeve connections were used
strengthening measure that proved to be effective in retro- to connect the longitudinal bars of the UHTBC elements
fitting slender shear walls, through the effective use of a and the RC foundation (JG/T 398-2012 [2012]). No connec-
high-performance material, with the significant advantage tion was set for Specimen UCSW-4, which was designed
that the overall wall configuration is only slightly affected. to explore the possibility of no-connection retrofitting for
Furthermore, it also improves the bending capacity without ancient and protected buildings.
damaging the inner portions. The UHTBCs can be prefab-
ricated and assembled on site, and even replaced at a later Material properties
stage if needed. The developed predictive equation of the The UHTCC used in the boundary columns of the tested
moment capacity also provides a reference for design. shear walls was composed of: PO 52.5 cement, fly ash, silica
fume, and fine sand; water-reducing additive; water; and
EXPERIMENTAL TESTS high-strength polyvinyl alcohol (PVA) fibers (fiber content:
Specimen design 2%). The literature demonstrates that, as multiple cracking
The experimental shear-wall specimens were designed phenomena occur, UHTCC exhibits tensile strain-hardening
referring to an existing high-rise building built in the 1990s in behavior and attains a maximum strain capacity larger than
Shanghai in a second-category seismic site. In this area, the 0.03, even with typical fiber volume fractions of 2% or less
earthquake intensity level is Grade 7, which corresponds to a (Xu and Cai 2010; Li et al. 2009). The main properties of the
design peak ground acceleration (PGA) = 0.35g (g is gravity high-strength PVA fibers include: length = 12 mm; diameter
acceleration) having an exceedance probability of 10% in = 39 μm; tensile strength ftc = 1620 MPa; tensile modulus of
50 years. A half-scale shear wall (SW-1) was designed and elasticity Ef = 42.8 GPa; ultimate strain = 0.07; and density
constructed with plain concrete as a reference wall, repre- = 1.3 g/cm3.
senting the shear wall satisfying the current seismic design Figure 2 shows the specimens and test views of dog-bone-
code (GB 50011-2010 [2010]). Three half-scale shear- shaped UHTCC specimens for tensile strength and cubic
wall specimens (UCSW-2, UCSW-3, and UCSW-4) were specimens of both plain concrete and UHTCC for compres-
designed with the same geometry and reinforcement as the sive strength. Figures 3(a) and (b) show the stress-strain
reference shear wall but strengthened with UHTBCs. For all curves of the tested material specimens. The main proper-
specimens, the 1:2 scale factor applies to all geometry quan- ties, including the elastic modulus and compressive strength
tities, while the materials are designed to be the same as in of the plain concrete and UHTCC, are presented in Table 2.
the original structure. Therefore, the load-bearing capacity Longitudinal reinforcement to the boundary columns was
is 1:4 of the original structure. Detailed descriptions of the formed from HRB400, and the main material properties of
specimens are presented in Table 1. The boundary column @seismicisolation
@seismicisolation
Fig. 2—Dog-bone tensile and cubic compressive tests: (a) dog-bone-shaped specimen and tensile test; (b) UHTCC specimen;
and (c) plain concrete specimen.
specimen was seriously damaged, and it decreased margin- bars’ dowel action, shear cannot be transferred between the
ally faster than that of UCSW-2 near the “no collapse” stage. boundary columns and core wall, which eventually end up
This indicates that the seismic performance of UCSW-3 carrying the horizontal force independently. To obtain an
satisfied the requirements of the current seismic design code. improved energy dissipation effect, connecting shear keys
Moreover, the specimen with a weak connection are essential (Park 2006).
(UCSW-4) exhibited high toughness and collapse-resistance
capacity in the experiment. Although the lateral force of Crack development
UCSW-4 was smaller than those of UCSW-2 and UCSW-3 The toughness of FRC composites is essentially due to
before the drift ratio reached 1/100, it was not smaller than the cracking prevention by the bridging effect of the fibers.
that of SW-1. The lateral bearing force of UCSW-4 at a For this reason, the crack growth during the experiment was
drift ratio of 1/50 was 34.3% and 5.6% larger than those of investigated. Figure 9 shows the crack distribution profile of
SW-1 and UCSW-3, respectively. The results indicate that each specimen at the yield, peak, and ultimate states and at
the seismic performance was significantly improved through drift ratios of 1/200 (~10 mm) and 1/100 (~20 mm), around
the connection-strengthening strategy involving UHTBCs, which is approximately “medium damage” and “serious
particularly in the large-deformation stage. damage.” The fractal index D (Mandelbrot 1982), as in
The hysteresis curves of Specimens UCSW-2, -3, and -4 Eq. (1), was used to evaluate the crack distribution as follows
show a pinching effect caused by slippage occurring at the
interface between the boundary columns and core wall. As ln N(r)
_
demonstrated by other authors (for example, Park [2006] D = − lr→0
im (1)
ln r
and Sørensen et al. [2017a,b]), shear keys can be prevented
from slipping by adopting concrete key slots and shear where r represents the size of a square grid cell; and N
reinforcing bars, which, however, show different damage represents the number of grids with cracks inside. The
modes while differently affecting the shearing capacity. vertical cracks appearing at the interface between the
In this study, Specimens UCSW-2 and -3 have shear key boundary column and the core wall are caused by a weak
slots and shearing reinforcing bars, as shown in Fig. 8(a). interface connection, which did not significantly impact
The shear reinforcing bars are placed at 160 mm spacing. the final failure and can be improved in actual use; thus,
Specimen UCSW-4 is constructed directly by precast and they were disregarded. In addition to the fractal index, the
on-site cast concrete, as shown in Fig. 8(b). In principle, crushing ratio, defined as the crushing area divided by the
under shear force, the concrete shear strength at the inter- total shear-wall surface area, was extracted from the images
face is firstly overcome, then the concrete key slots and for each specimen and was used to evaluate the concrete-
shear reinforcing bars continue to bear the shear force. In crushing performance at ultimate (Table 5).
the tests, the interfaces of Specimens UCSW-2 and UCSW-3 Figure 9 reveals smaller horizontal cracks in the UCSW-2
slightly cracked along the key, indicating that key slots and specimen than in SW-1. The fractal indexes of the three
shear reinforcing bars play an important role in the ultimate UHTBC shear walls were significantly larger than that of
state and that they are activated after concrete cracks. The SW-1, and the fractal-index growth rate of SW-1 was signifi-
pinching effect is more significant for UCSW-4 because cantly higher than those of the UHTBC shear walls before
there are no connecting bars to resist shear. After the shear the lateral drift ratio reached approximately 1/133. Subse-
strength is overcome and the concrete cracks, the aggregate quently, the fractal-index growth rate was lower than those
interlock contribution is weak and, without the reinforcing of the UHTBC shear walls. The crushing areas of UCSW-2
and
@seismicisolation
@seismicisolationUCSW-3 were significantly reduced compared with
@seismicisolation
@seismicisolation
where H is the wall height; and h is the section depth. where λ0 = 2 is the aspect ratio; ft = 4.5 MPa is the concrete
The curvatures at yield and at the peak lateral force are tensile strength; bw and h0 are the shear-wall section width
shown in Fig. 12, as compared to those calculated according and length; N = 600 kN is the axial force; Aw and A are the
to the model (Bohl and Adebar 2011) and those obtained shear web and entire section area, respectively; and Ash =
from the test results. The experimental yield and peak point 56 mm2, s = 125 mm, and fyh = 501 MPa are the trans-
curvatures are obtained from strain gauges applied on the verse reinforcement area, its vertical spacing, and its yield
outer steel bars along the wall height. The corresponding strength, respectively.
lateral displacements are approximately 14 to 17 mm and According to the aforementioned codes, the average
32 mm (the lateral top drift is 1/62.5), respectively. The yield shear capacity of the four specimens tested in this study is
curvature is calculated when the outer strain reaches the estimated as 895 kN and 445 kN, respectively. It is worth
yield value (0.002), while the peak point curvature is found noticing that these two values are significantly different, thus
when compressive concrete attains its strength according proving the large discrepancies in the shear capacity equa-
to the model. The peak point curvature depends on the tions of different codes (De Domenico et al. 2023). Nonethe-
boundary length, so the values obtained for UCSW-3 and less, because the shear capacity computed with both equa-
-4 are slightly lower than SW-1 and UCSW-2. The plastic tions is higher than the bending capacity exhibited in the
hinge length is obtained by Eq. (4) as 360 mm for all four tests, which is at most 400 kN, this check confirms that the
specimens, as shown in Fig. 12. shear walls were correctly designed to prevent shear failure.
In all four specimens, the yield curvature distribution is in As far as the flexural capacity is concerned, both the
good accordance with the model results for SW-1, UCSW-2, Chinese code (GB 50010-2010) and ACI 318-19 are based
and UCSW-3, which indicates that the model is reasonable on the cross-section equilibrium and stress-strain relation-
for shear walls with reliable shear keys. However, for spec- ship at nominal strength. In both, it is assumed that: 1) the
imen UCSW-4, the yield curvature is lower than the model section strain is linearly distributed; 2) the longitudinal rein-
results. The weak interface connection between the boundary forcing bars in the tension boundary column have yielded;
@seismicisolation
@seismicisolation
be satisfied, and the vertical force equilibrium can be stated Because the post-cracking effect of FRC shows a signif-
as icant impact on the seismic performance, two models are
discussed and compared with the currently used model. As
N = αβ1fccbwl0 + 0.5β2fcbwlc + Fs – Fs′ – fsw′Asw′ (7) stipulated in the codes, the post-cracking tensile strength of
concrete is not considered in the sectional bearing capacity
where α is the ratio of the peak stress to the concrete (ACI 318-19; GB 50011-2010), and the fiber-reinforcement
strength, equal to 1.0 for concrete strength lower than effect is not considered in Model 1. In Model 2, the boundary
50 MPa; β1 and β2 are the stress block shape factors of columns’ steel bars are assumed to reach the ultimate
compressive concrete, equal to 0.8 for concrete strength strength according to the actual steel bar strain. In Model 3,
lower than 50 MPa, according to GB 50011-2010; the first the post-cracking FRC tensile strength is also accounted for
term αβ1fccbwl0 represents the compression resultant of based on the assumption of Model 2. The descriptions of the
the assumed uniformly distributed concrete stress in the models are listed in Table 6.
boundary column; and the second term 0.5β2fcbwlc represents In Model 2, the section moment is
the compression resultant of the linearly distributed concrete
@seismicisolation
@seismicisolation
*
Post-cracking tensile strength is not considered.
M = fsu′As′(hw0 – l0) + ΣfsAsls + Σfsw′Asw′lsw′ the reduction factor of the fiber-reinforced concrete of post-
(10) cracking tensile strength is approximately 0.5 and the error
_1
( _lc _
)
l0 _1
+ 2 β2fcbw lc 3 + 2 + 2 N(hw0
− l0 ) is between 3.8 and 10%. When the reduction factor is 0.75,
the error is between 0.7 and 7.3%.
In Model 3, the equilibrium of vertical force gives As the calculation comparison based on different failure
patterns indicates, the tensile strength of plain concrete
N = αβ1fccbwl0 + 0.5β2fcbwlc + ΣfsyAs – Σfsu′As′ shows a minor contribution to the bearing capacity, which
– fsw′Asw′ – γfctbwl0 (11) meets with the current assumption. When the boundary
column concrete is replaced by the UHTCC, the strength of
The section moment is then steel bars in the boundary area and the post-cracking tensile
strength of UHTCC contribute significantly to the lateral
M = fsu′A
s′h w0
_l0
(
−2 + ∑fsA )
sls+ ∑fsw′A
sw′lsw′
(12)
bearing capacity. With the fiber confinement, the edge longi-
tudinal steel bars are more stable against buckling. Conse-
+ γfct bw l0 ( hw0
_1
( _
)
lc _
l0 _1
− l0 ) + 2 β2fcbw lc 3 + 2 + 2 N( hw0
− l0 ) quently, the boundary column steel bars are more likely to
reach the ultimate strength.
where fs and fsu′ are the yield and ultimate strengths of longi- DISCUSSION
tudinal steel bars in the boundary area on the tension and Strengthening effect of UHTBC
compression sides; As′ and As are the area of the boundary-area The experimental results indicated that replacing plain
steel bars on the tension and compression side, respectively; concrete in the boundary columns with the UHTCC can
ls′ is the distance from the rotation center; γ is the reduction significantly improve the seismic performance of shear
factor of the tensile strength of FRC in the confined boundary walls.
zone; and Σfsw′Asw′lsw′ represents the moment caused by The UHTBC used in this study demonstrates its effective-
distribution bars, which can be calculated according to the ness for strengthening the boundary elements of shear walls
plane sections assumption. by mitigating the local cracks, reducing their configuration,
Table 7 shows the calculated maximum section moment and preventing rapid crack development and local concrete
through Models 1 to 3. The sectional bearing moment and crushing. The shear walls with UHTBC elements exhibited
the corresponding lateral bearing force are also listed for improved lateral resistances and deformation capacities. The
comparison. In Model 3, the reduction factor γ of the tensile peak lateral force and the corresponding lateral displace-
strength of FRC in the confined boundary region is set as 0.5 ment of UCSW-2 were 22% and 19% larger, respectively,
and 0.75, separately. than those of SW-1. The UHTBC specimens maintained a
As Table 7 shows, the maximum lateral bearing capacity good resistance capacity even when the lateral bearing force
obtained from Model 1 is acceptable for plain concrete shear dropped to 70% of the peak value. The secant stiffness of
walls, but the error is 23 to 32% for the FRC-RC shear walls. UCSW-2 was significantly higher than that of SW-1 (by
When considering the ultimate strength of the outer-edge 51%, 13%, and 95% at drift ratios of 1/500, 1/100, and 1/50,
steel bars, the error is lower yet still 10 to 17%. In Model 3, respectively). The UHTBC elements increased the energy
considering the concrete tensile strength, for the two kinds dissipation capacity, particularly in the large-deformation
of confined column shear walls, the error is significantly stage. The energy dissipations of the single hysteresis loops
reduced. As the literature suggested (Zhang et al. @seismicisolation
2022),
@seismicisolation
@seismicisolation
@seismicisolation
Cyclic tests are conducted on interlinked reinforcing bar resisting large bending moments, which may lead to bond-
coupler-box assemblies, adopted to retrofit buckled reinforcing slip failure. Welding can be an alternative technique for
bars at the plastic hinge locations of columns in multi-storied connecting broken reinforcing bars with new bars. However,
reinforced concrete building frames. The efficacy of the proposed it also requires removing a large amount of intact concrete
retrofitting technique is evaluated by comparing the hysteresis
on both sides of the damaged section to expose the required
behavior, computed parameters of performance index, and failure
weld length and an all-around clear space for functioning.
mechanism of the reconstructed frame with the original frame. An
energy-based strength deterioration hysteresis model is developed The length of lap weld for connecting reinforcing bars as per
on the basis of cyclic test results for analytically computing the strength, and codal requirements, IS 9417:2018 (2018), is
post-yield behavior of retrofitted reinforced concrete (RC) frame approximately 600 mm (23.62 in.), which is higher than the
with the proposed coupler-box assembly. The experimental test concrete removal area required for insertion of the proposed
results manifest that the coupler-box assembly can be a promising coupler sleeve, 500 mm (19.69 in.)—that is, two times the
futuristic approach for seismic retrofitting of severely damaged proposed coupler sleeve length. Further, welding causes
reinforced concrete buildings, where buckling of longitudinal rein- metallurgical changes with the recrystallization of micro-
forcing bars at the plastic hinge location of columns is inevitable, structures in particles and is not preferred in the potential
and the process of restoration is challenging under existing gravity hinge region of primary structural members. The study on
loads. The suggested retrofitting mechanism restrains the section
mechanical properties of welded deformed reinforcing steel
from any movement against rotation and helps in shifting the yield
bars conducted by Ahmed (2015) shows that the strength
location of reinforcing bars. The main advantage of adopting the
coupler-box is that there is no observed slip of reinforcing bar from and elongation of welded reinforcing bars decrease up to
the sleeve, and the entire retrofitted section remains intact even 40% and 60%, respectively.
after a lateral storey drift of 6%, which is larger than the collapse The use of reinforcing bar couplers or coupling sleeves
prevention drift level of 4% as per Federal Emergency Manage- in joining the ends of reinforcing bar segments is a preva-
ment Agency guidelines. lent practice in the modernized construction of buildings or
bridges, which safeguards the reinforcing materials at their
Keywords: coupler-box assembly; deterioration strength hysteresis model; lap-splice regions. The invention of grouted splice sleeves
energy dissipation; hysteresis behavior; reinforcing bar coupler sleeve;
seismic retrofitting.
dates back to the late 1960s by Yee (1970, 1986, 2009).
Several studies in the past on coupler/coupling sleeves are
INTRODUCTION based on their shape, diameter, length, material character-
Buckling/yielding of longitudinal reinforcing bars at the istic, bonding properties of the grout material, anchorage
plastic hinge locations of multi-storied building frames under capacity, and so on. Hybrid couplers with a combination
severe earthquake excitations is one of the most common of mechanical and grouting techniques were tested by
failure modes. As a result, the column member shortens Hope (1987). Threaded end reinforcing bars into grout-
and loses its capability to carry even its prescribed gravity filled coupler sleeves were suggested by Lancelot (1995).
load, which sometimes leads to the progressive collapse of Mechanical couplers in the form of tubular members with
buildings. Seismic retrofitting of buckled reinforcement at and without internal threading, along with a distinctive
the hinge location of columns in such multi-story buildings number of external bolts, were proposed by Holdsworth
proves to be challenging as there is no mechanism to realign (1997, 2006). Annular-shaped seals at the ends of steel
the buckled bars under high compressive loads. Under these sleeves for locking reinforcing bar elements along with
circumstances, removal of the buckled portion of rein- grout matter were recommended by Mochizuki and Nihei
forcing bar is inevitable, and new reinforcement is intro- (1998). High-strength concrete grout in addition to large-
duced either by welding, splicing, or through mechanical headed reinforcing bars through pipe couplers was suggested
couplers. The conventional method of lap splicing at hinge
ACI Structural Journal, V. 120, No. 5, September 2023.
locations is neither recommended nor practical as it requires MS No. S-2021-120.R7, doi: 10.14359/51738831, received March 4, 2023, and
large demolition of intact concrete in an existing structural reviewed under Institute publication policies. Copyright © 2023, American Concrete
Institute. All rights reserved, including the making of copies unless permission is
member on either side of the cut reinforcement. Moreover, obtained from the copyright proprietors. Pertinent discussion including author’s
closure, if any, will be published ten months from this journal’s date if the discussion
splicing reinforcing bars at hinge locations is incompetent in is received within four months of the paper’s print publication.
@seismicisolation
@seismicisolation
of reinforcing bars. This leads to a conventional failure or coupler-stirrup column, which is absent in the coupler-box
rupture in reinforcing bars of the coupler-stirrup column column.
with an asymmetric behavior. The irregular behavior in Fig. Interlinking the proposed coupler sleeves essentially
6(a) is due to the rupture of reinforcing bars at the right edge unites the strength of the standalone coupler under flexure,
of the column face—that is, for classification, rupture occurs which minimizes the chance of flexural yielding and the
in reinforcing bars that experience tensile stresses under restored section is subjected to shear force where the
the negative half-cycle, which causes a sudden reduction strength of the coupler sleeve is resistive. Thus, the entire
in lateral strength post 3.2% lateral drift. This aspect is not box assembly functions as a rigid link section which shrinks
noticed during the negative half-cycle of testing. Moreover, the plastic hinge length and compels the adjacent rein-
the pinching effect in the hysteretic loops is noticed in the forcing bar segments to yield away from the coupled region.
Energy dissipation capacity of the box assembly column is
@seismicisolation
@seismicisolation
( )
post their elastic state are computed and plotted against
βi = i (5)
corresponding lateral drifts, as shown in Fig. 18. The varia-
E
t −
∑
E
j
j=1
tion in stiffness degradation pattern is similar in both model
frames, with a marginally reduced rate in the RET frame. where Ei represents the hysteretic energy dissipated in the
The variation in strength degradation parameter of frame i-th cycle; ∑Ej indicates the summation of energy dissipated
models shows a deviation in the RET frame post 5% drift by in all previous cycles until i; and Et is the energy dissipation
virtue of its declination in lateral strength capacity. Thus, it capacity expressed as a function of twice the elastic strain
is affirmed that the RET frame regains strength and stiffness energy at yielding—that is, Fyδy
up to its original level by adopting the recommended retro-
fitting technique. Et = γFyδy (6)
Fref,i
− = (1 − βc,i) Fref,i−1
− (12)
where Fref,i
+/− is the deteriorated reference strength after expe-
dition i; Fref,i−1
+/− is the deteriorated reference strength before
expedition i; and βc,i is computed using the relevant γc factor
in the aforementioned equations.
The computed parametric values, along with the cyclic
degradation modes of basic and post-capping strength for
tested RC portal frames, are presented in Fig. 22, where
the value of degradation parameter γ is calibrated based
on the experimentally obtained hysteresis behavior of the
individual frame. The value of γ is high for the CON frame
specimen owing to its larger cumulative energy dissipation
capacity. The model ductility (δc/δy) is relatively high by 7%
for the CON frame, whereas the peak strength of the model
backbone is comparatively high by 7% for the RET frame.
The post-capping stiffness coefficient (αc) for the CON frame
model is half the value of the RET frame, which signifies the
Fig. 20—Backbone curve for hysteretic models (Ibarra and gradual post-peak behavior. However, the strain-hardening
Krawinkler 2005).
Fig. 21—Peak-oriented models with computed deterioration modes (Ibarra and Krawinkler 2005).
@seismicisolation
@seismicisolation
@seismicisolation
@seismicisolation
Check out the Chapters Special Section from the November 2020 Concrete International:
www.concrete.org/publications/concreteinternational.aspx
Student Chapters
Join or form an ACI Student Chapter to maximize your influence, knowledge
sharing, and camaraderie! ACI has 240+ student chapters located throughout
the world, each providing opportunities for students to:
• Connect with their peers and participate in concrete-related activities such as: student
competitions, ACI Conventions, ACI Certification Programs, ACI Educational Seminars, local
chapter meetings, social events, and community service projects;
• Network with members of local chapters, many of whom have been in the industry for
decades and can help to develop professional relationships and offer career advice;
• Win recognition for their universities through the University Award; and
• Learn about the many scholarships and fellowships offered by the ACI Foundation and by
ACI’s local chapters.
@seismicisolation
@seismicisolation
www.concrete.org/chapters
ACI STRUCTURAL JOURNAL TECHNICAL PAPER
Title No. 120-S80
Large-diameter longitudinal reinforcement (D43 [No. 14]) and consists of one central large spiral and four small spirals
multi-spiral reinforcement were examined for use in the construc- at the corners of the column cross section. The reinforced
tion of concrete columns to reduce labor demand. Two full-scale concrete design code of Taiwan14 is the only design code
columns were tested to investigate their seismic performance. The globally that includes the five-spiral reinforcement. The code
test column was reinforced with large-diameter longitudinal rein-
requires that the maximum distance between the inner faces
forcement (D43 [No. 14]) and the control column with normal-
of small and large spirals in the overlapping region shall not
diameter longitudinal reinforcement (D36 [No. 11]). Both columns
were laterally confined by five-spiral reinforcement. Test results be less than 0.3 times the inner diameter of the small spiral
showed that the seismic performance of the test column was similar and 60 mm (2.36 in.) to ensure proper interlocking between
to the control column. Both columns showed ductile behavior with small and large spirals. Moreover, the volumetric ratio of
good energy dissipation. Based on the test observations, a buckling each spiral of five-spiral reinforcement in the plastic hinge
model was developed to predict the buckling behavior of longitu- region shall satisfy the requirements for conventional spiral
dinal bars laterally supported by five-spiral reinforcement. Based reinforcement.7,14
on the proposed model, it was found that for the same amount and The large-diameter longitudinal reinforcement used
within the practical spacing range of five-spiral reinforcement, the in this research has a diameter of 43 mm (1.69 in.) (D43
calculated buckling slenderness ratio and hence the compressive [No. 14]). Large-diameter longitudinal reinforcement can
stress-strain behavior of D43 (No. 14) bars is similar to that of
reduce the number of mechanical splices and hence decrease
D36 (No. 11) bars.
labor demand associated with splicing operations. Previous
Keywords: buckling; columns; cyclic test; five-spiral reinforcement; studies have investigated the bond behavior of large-diam-
large-diameter bars; multi-spiral reinforcement; reinforced concrete; eter bars. Ichinose et al.15 studied the bar size effect on the
seismic. bond strength of bars with a diameter up to 52 mm (2.05 in.).
Research results showed bond strength tended to decrease
INTRODUCTION with increasing bar diameter. However, this size effect
Due to the shortage of workers, the construction industry decreased with increasing confinement to the bars provided
in Taiwan has been developing methods to decrease labor by concrete or transverse reinforcement. Steuck et al.16 inves-
demand. In this research, two ways to reduce labor demand tigated the bond behavior of bars grouted in corrugated steel
are employed for the construction of reinforced concrete ducts for a bar diameter up to 57 mm (2.24 in.). Research
columns, including multi-spiral transverse reinforcement results showed that the bar size effect on bond strength was
and large-diameter longitudinal reinforcement. insignificant. Note that good confinement was provided for
Multi-spiral reinforcement is developed to allow non- bars tested in Steuck’s study. The bond failure was governed
circular columns to use and take advantage of spiral rein- by shearing along the bar-grout interface, which reduced
forcement. Several forms of multi-spiral reinforcement the size effect according to Ichinose’s study. The bond tests
have been developed and studied, including two-spiral rein- by Ichinose et al. and Steuck et al. were carried out using
forcement for oblong columns,1-6 five-spiral reinforcement monotonic loading. Murcia-Delso et al.17 studied the bond
for square columns,7-9 six-spiral reinforcement for rect- behavior of bars with a diameter up to 57 mm (2.24 in.) using
angular columns,10 seven-spiral reinforcement for oblong cyclic loading. The bars were embedded and well confined
columns,11-13 and 11-spiral reinforcement for rectangular in a large circular concrete member reinforced with a dense
columns.12 It has been demonstrated that multi-spiral rein- transverse reinforcement. This simulated the confined situa-
forcement can be fabricated by a machine and assembled tion for a column longitudinal bar embedded into an enlarged
rapidly on a construction site, significantly reducing labor pile shaft. A development length equation was proposed for
demand.7 Moreover, multi-spiral reinforcement shows large-diameter bars based on the test results and analytical
a better confinement effect than rectilinear tie reinforce- study.
ment. Even with a smaller amount of transverse reinforce-
ment, columns with multi-spiral reinforcement show better
seismic performance than columns with rectilinear trans- ACI Structural Journal, V. 120, No. 5, September 2023.
MS No. S-2021-347.R2, doi: 10.14359/51738834, received April 5, 2022, and
verse reinforcement. In this research, five-spiral reinforce- reviewed under Institute publication policies. Copyright © 2023, American Concrete
ment is used. Figure 1 shows the five-spiral column tested Institute. All rights reserved, including the making of copies unless permission is
obtained from the copyright proprietors. Pertinent discussion including author’s
in this research. It can be seen that five-spiral reinforcement closure, if any, will be published ten months from this journal’s date if the discussion
is received within four months of the paper’s print publication.
@seismicisolation
@seismicisolation
Fig. 4—C36 damage pattern at: (a) 1.5% drift ratio; (b) 3% drift ratio; (c) buckling of longitudinal bars (7% drift ratio); and
(d) fracture of longitudinal bars (8% drift ratio).
At the 7% drift ratio, for C36, buckling of longitudinal strength dropped to 47% of the peak strength during the
bars occurred during the third cycle of loading, as shown 8% drift loading. The testing was terminated at the 9% drift
in Fig. 4(c). The maximum lateral strength of the column ratio when the lateral strength dropped to 38% of the peak
dropped to 83% of the peak strength during the 7% drift strength.
loading. At the 8% drift ratio, fracture of longitudinal bars For C43, buckling of longitudinal reinforcement also
occurred, as shown in Fig. 4(d). The maximum lateral occurred at the 7% drift but during the second cycle, slightly
@seismicisolation
@seismicisolation
Fig. 6—Force and displacement diagram for: (a) C36; and (b) C43. (Note: 1 kN = 0.22481 kip; 1 mm = 0.0394 in.)
earlier than C36. The maximum lateral strength decreased The relationship between the lateral force and drift ratio
to 84% of the peak strength during the 7% drift loading. (hysteretic behavior) of C36 and that of C43 are shown in
Figure 5(c) shows the damage condition of the column at Fig. 6. Note that the axial loading system tilted as the column
the 7% drift ratio. At this drift ratio, fracture of five-spiral deformed laterally during testing, causing an additional
reinforcement was observed. This was not found in C36. moment at the column base. This additional moment was
Fracture of longitudinal bars also occurred at the 8% drift removed from the force presented in Fig. 6. Both columns
ratio, causing the lateral strength to drop to 34% of the peak showed ductile behavior with drift capacities significantly
strength during the 8% drift loading. The testing was stopped higher than 3%, a seismic performance target used in the
at the 9% drift ratio when the lateral strength dropped to literature for columns.21 The difference in behavior between
23% of the peak strength. the two columns was not significant until the 7% drift ratio,
at which the longitudinal bars started to buckle. After the
@seismicisolation
@seismicisolation
Fig. 7—(a) Equivalent viscous damping ratio; (b) strain of longitudinal reinforcement; (c) strain of large spirals; and (d) strain
of small spirals.
7% drift ratio, the lateral strength of C43 degraded faster material properties,19 and Mtest/Mn are also listed in Table 2.
than that of C36. The number of longitudinal bar fractures C43 exhibited a slightly higher overstrength (Mtest/Mn) than
for C36 was the same as that for C43. However, because C36. In general, the performance difference between the two
the area of one longitudinal bar of C43 was larger than that columns was not significant.
of C36, the strength degradation for C43 related to each bar
fracture was higher than C36. Energy dissipation and reinforcement strain
The envelope of the hysteretic behavior was idealized The hysteretic energy dissipation in terms of the equiv-
by a bilinear relationship based on FEMA 365.22 The first alent viscous damping ratio for both columns is shown in
segment of the bilinear relationship intersected the envelope Fig. 7(a). C43 exhibited an energy dissipation capacity
response at 60% idealized yield force Vy, which is the force similar to C36. The use of D43 (No. 14) longitudinal bars
at the idealized yield point. The second segment ended at did not significantly alter the energy dissipation capacity of
the ultimate drift ratio ∆u, defined as the drift ratio when the the column. Fig. 7(b), (c), and (d) show the strain responses
envelope response drops to a point with a lateral force equal of the longitudinal reinforcement, the large spiral, and the
to 80% of the peak lateral strength. The intersection of the small spiral, respectively, from the plastic hinge region of
two linear segments, which defines the idealized yield point, the columns. The locations of the strain gauges used in these
was selected so that the area below the envelope response plots can be found in Fig. 2. The longitudinal reinforcement
and below the idealized bilinear relationship were similar. of C43 generally showed a smaller strain response than C36
The displacement ductility μ of each column is defined as for the same drift ratio. The longitudinal reinforcement of
the ratio of ∆u to the displacement of the idealized yield C43 yielded at the 1% drift ratio, later than C36, which
point ∆y. The Vy, ∆y, ∆u, and μ of the bilinear relationship occurred at the 0.75% drift ratio. However, at the 1.5% drift
for the two columns are listed in Table 2. C36 had a slightly ratio, the longitudinal reinforcement of both columns showed
larger displacement capacity (∆u and μ) than C43. The similar strain responses. The strain gauges of both columns
maximum column moment measured from testing Mtest, the failed later at this drift ratio. The strain responses of the large
nominal moment strength Mn calculated based on measured and small spirals of C43 were similar to those of C36. The
@seismicisolation
@seismicisolation
tial energy of the model consists of the energy stored in the = 0 ⇒ P = + ∑ 1 − cos (3)
∂δ L3 2π 2 i =1 n
buckled bar, the energy stored in the equivalent springs, and
a decrease in the potential energy due to work done@seismicisolation
by the
@seismicisolation
=
1 L
∫ EI (v ′′ ) dx + k ∑ 1 − cos
2 1 L
− 2 P ∫0 (v ′ ) dx
2 bars (shown in blue color) in each of the two small spirals
2 0 2 2 i =1 L in the compression side and the one bar (shown in red color)
(5) in the large spiral are prone to buckle. Bars in overlapping
regions are assumed not to buckle because they are further
where sa is additional buckling length beyond (n – 1) × s. away from the extreme compression fiber and are subjected
Fig. 9—(a) Locations of buckled bars; (b) deformed parts of small and large spirals; and (c) equilibrium to determine pushing
force from buckled bar. (Note: Full-color PDF can be found at www.concrete.org.)
@seismicisolation
@seismicisolation
to strong confinement from a large and a small spiral. More- sectional area of a spiral. For the large spiral of C36 and C43,
over, because the overlapping region is subjected to strong nb = 1 and 0, respectively. For the small spiral of C36 and
confinement and is typically designed with longitudinal bars C43, nb = 2. As stated previously, the spirals of the columns
to enhance the capability of the region to interlock the large tested yielded at the 2 to 3% drift ratios, earlier than the
and small spirals, the overlapping region is assumed to act as buckling of the longitudinal bars, which occurred at the 7%
a stiff column and hence provides a strong anchorage for the drift ratio. This supports the use of the post-yield modulus of
small spiral anchored to the region. As a result, buckling of elasticity for spirals.
the two longitudinal bars in each small spiral is assumed not The stiffness of the equivalent spring k for the large or
to affect the large spiral. Due to different constraints from the small spiral can be obtained by the following equation with
large or small spiral and different numbers of buckled bars F from Eq. (9)
within the large or small spiral, the stiffness of equivalent
θ
springs for the large spiral and that for the small spiral need
to be calculated and compared. The bar within a spiral with F
( )
Et At cos ( π − θ m ) + 1 m
2π (10)
k= =
a smaller stiffness buckles first. δ Dnb
When the bar in the large or small spiral buckles, the large
where k is equivalent spring stiffness of a spiral. For the large
or small spiral is pushed to deform outwardly, as shown in
spiral of C36, k = 563.34 MPa-mm (3216.75 psi-in.) based
Fig. 9(b). Due to the constraint by concrete, only a portion
on Eq. (10). The large spiral of C43 has no longitudinal bars
of the spiral is deformed. The deformed part of the large
for buckling consideration (nb = 0). For the small spiral of
spiral is assumed to be one-half of the large spiral and has
C36 and C43, k = 284.15 MPa-mm (1622.54 psi-in.) based
a central angle θm of 180 degrees. This assumption is the
on Eq. (10). Therefore, the governing value of k for both C36
same as that used by Dhakal and Su.26 The effect of the small
and C43 is 284.15 MPa-mm (1622.54 psi-in.). And, for both
spiral is ignored. The deformed part of the small spiral is
columns, buckling of longitudinal bars in small spirals is the
assumed to be one-half of the small spiral, excluding the
governing buckling behavior.
overlapping region, as shown in Fig. 9(b). The portion of the
small spiral within the overlapping region is assumed not
Comparison with test observation
to deform due to the strong confinement effect. The central
Table 3 lists the values of important parameters for the
angle of the deformed part is smaller than 180 degrees. With
analytical buckling behavior of the D36 (No. 11) and D43
these assumptions, the tensile strain developed in the large
(No. 14) bars of C36 and C43, respectively. As stated previ-
or small spiral due to the push from the buckled bar can be
ously, the governing value of k is from the small spiral.
calculated as
This means buckling of longitudinal bars within the small
spiral governs the buckling behavior. This is consistent with
θm ( D + δ) − θm D θmδ the experimental observation, as shown in Fig. 4 and 5 for
ε = = (8)
2πD 2πD C36 and C43, respectively. The longitudinal bars within the
large spiral of C36 did not show significant buckling until
where ε is strain of a spiral due to bar buckling; θm is angle
the end of testing. Note that for C43, no longitudinal bars
corresponding to the deformed part of a spiral; and D is
were designed within the large spiral close to the extreme
diameter of a spiral. For the large spiral of C36 and C43,
compression fiber of the section (Fig. 2). Experimental
θm = 180 degrees and D = 820 mm (32.28 in.). For the small
observation also showed that longitudinal bars within the
spiral of C36 and C43, θm = 143 degrees and D = 240 mm
overlapping regions did not expose and were still well
(9.45 in.).
confined by surrounding concrete and spirals until the end
From the free-body diagram shown in Fig. 9(c), the
of testing. This is consistent with the assumption used in
pushing force from the buckled bar F can be determined by
the proposed buckling model that the overlapping regions
equilibrium with ε from Eq. (8)
provide a strong anchorage for small spirals to restrain buck-
ling of longitudinal reinforcement.
θm δ
( ) ( )
nb F = Et At ε cos ( π − θ m ) + 1 = Et At cos ( π − θ m ) + 1
2πD
The calculated buckling lengths are 477.70 and 582.59 mm
(18.81 and 22.94 in.) for C36 and C43, respectively, close
(9) to measured buckling lengths—that is, 487.5 and 562.5 mm
(19.19 and 22.14 in.), respectively—from the experiment.
where nb is number of buckled bars; Et is post-yield Because the test data are still limited, more tests are needed
modulus of elasticity of spiral steel, taken as 1.5% the to further validate the applicability of the proposed model to
initial modulus of elasticity of spiral steel ; and At is cross-
27 the buckling behavior of longitudinal bars laterally supported
by five-spiral reinforcement. According to Dhakal and
@seismicisolation
@seismicisolation
@seismicisolation
@seismicisolation
@seismicisolation
@seismicisolation
Building the Future
ACI STRUCTURAL JOURNAL TECHNICAL PAPER
Title No. 120-S81
Most of the studies conducted on the rehabilitation of reinforced to establish a strength hierarchy between members framing
concrete (RC) beam-column joints are on pre-1970 structures. into the joint, thereby preventing brittle joint failure in the
Recently, it was reported that seismically designed beam-column column and promoting ductile beam failure (strong-column,
joints might also suffer damage under lateral loading. On the other weak-beam).8 Several techniques such as concrete jacketing,
hand, there is an increasing interest among researchers to study
steel jacketing, haunch solution, fiber-reinforced polymer
the effectiveness of geopolymer as an alternative repair material.
(FRP) wrapping, removal and replacement technique, and so
To date, no study has been conducted to examine the performance
of geopolymer for the rehabilitation of seismically detailed beam- on are available to rehabilitate damaged beam-column joint.
column joints following the removal and replacement method under Concrete jacketing has been found successful in restoring
cyclic loading. In the present investigation, two groups of exterior the cyclic performance of damaged beam-column joints.1
beam-column joints with different flexural strength ratios were However, concrete jacketing is a time-intensive process that
rehabilitated with geopolymer mortar. For comparison, another set increases the structural weight and size of the members.
of beam-column joints (one from each group) were rehabilitated Construction time decreases due to the use of prefabricated
with cement mortar following the same rehabilitation technique steel plates and steel haunch retrofit systems. Both tech-
and testing. Test results indicated that geopolymer rehabilitated niques change the brittle joint shear of control specimens
specimens exhibited 20 to 21% higher initial stiffness, 19 to 22% to ductile beam failure of rehabilitated specimens.9-11 Along
higher displacement ductility, 24 to 37% higher cumulative energy
with increased weight and size of repaired joint, corrosion
dissipation, 14 to 17% higher initial equivalent viscous damping
is a major concern. Moreover, drilling through slabs poses
ratio, 21 to 26% higher ultimate equivalent viscous damping ratio
at failure, and 10 to 14% lower damage index compared to spec- a practical implementation problem for haunch retrofit
imens rehabilitated with cement mortar. However, irrespective systems. FRP has shown exemplary performance in the
of repair material, removal and replacement technique was only rehabilitation of damaged beam-column joints.1 Neverthe-
able to partially restore the cyclic performance of rehabilitated less, high cost, debonding due to incompatibility between
specimens. the epoxy and concrete substrate, decreased performance in
high temperatures, and wet environments are some of the
Keywords: beam-column joint; cyclic loading; geopolymer; removal and demerits associated with FRP.
replacement method; seismically detailed.
FEMA 30812 suggests that the removal and replacement
technique should be carried out for rehabilitation of highly
INTRODUCTION
damaged specimens subjected to severe earthquakes. In
An abundant number of literatures are available on the
this technique, damaged concrete is removed and replaced
repair and strengthening of beam-column joints damaged
with high-strength, low-shrinkage repair material.1,13-16 Lee
under earthquake loading. Most of the studies are limited to
et al.13 adopted a removal and replacement method to repair
substandard pre-1970 structures with inadequate detailing.1
extensively damaged exterior beam-column joints under
With the advent of the seismic code of practice, beam-
severe earthquake loading. Concrete of the damaged joint
column joints are being constructed with special confining
was removed followed by subsequent replacement with
hoop reinforcement and adequate anchorage. Such seismi-
higher strength concrete (compressive strength 48 MPa
cally detailed beam-column joints exhibit better strength,
[6.96 ksi]). Experimental results showed that the removal
ductility, and higher energy dissipation capacity.2 Recently,
and replacement method was able to restore structural integ-
some studies revealed that reinforced concrete (RC) struc-
rity of the rehabilitated specimens up to control specimens.
tures designed as per seismic code of practice exhibited
Karayannis et al.14 removed the fragmented joint concrete
damage in the joint region when maximum force acting
of exterior beam-column subassemblies. Voids were filled
on the joint exceeded the code expectation level.3-5 Only a
with high-strength (83 MPa [12.04 ksi]), rapid-hardening,
limited number of studies are available on the rehabilita-
and low-shrinkage cement paste. When tested under cyclic
tion of seismically designed and detailed structures.3,6,7 The
loading, a significant increase in ultimate load, stiffness,
failure of such structures opened new research avenues in
the field of rehabilitation of ductile beam-column joints. ACI Structural Journal, V. 120, No. 5, September 2023.
Rehabilitation of damaged structures provides an econom- MS No. S-2021-494.R2, doi: 10.14359/51738835, received December 2, 2022, and
reviewed under Institute publication policies. Copyright © 2023, American Concrete
ically viable alternative to demolition and reconstruction. Institute. All rights reserved, including the making of copies unless permission is
obtained from the copyright proprietors. Pertinent discussion including author’s
The aim of seismic rehabilitation of beam-column joints is closure, if any, will be published ten months from this journal’s date if the discussion
is received within four months of the paper’s print publication.
@seismicisolation
@seismicisolation
with various molarities—8M, 10M, and 12M—of activator strength ratio, test specimens were grouped into two catego-
solution were prepared. The binder (GGBS + FA) to sand ries: Group 1 (flexural strength ratios = 1.2); and Group 2
proportion was maintained at 1:3. For rapid strength gain (flexural strength ratios = 1.4). The flexural strength ratio of
and high compressive strength, the ratio of alkali to binder a joint is defined as the ratio of summation of nominal design
was kept constant at 0.65.25 The codal provisions of IS strength of all the columns to the summation of nominal
1727-195726 were followed to prepare trial mixtures in stan- design strength of all the beams framing into the joint.32
dard 50 cm (7.75 in.) mortar cubes and to evaluate setting Such flexural strength ratios were achieved by varying beam
times. Test results of trial mixtures are provided in Table 2. and column reinforcement. The structural design of beam
According to Duan et al.20 and Yun and Choi,27 desired prop- and column components of specimens was done as per
erties for suitable repair materials are early strength, rapid provisions laid down in IS 456-2000.33 The joint panel of the
setting (low initial and final setting time), and high compres- specimens was designed and detailed as per codal provisions
sive strength. Among all the trial mixtures, T3 displayed the of IS 13920-2016.32 Special confinement reinforcement in
earliest strength gain, least initial and final setting time, and the form of rectangular hoop reinforcement was provided in
highest compressive strength. Therefore, T3 was selected as and around the joint region by taking into the provision of
the most suitable repair material to carry out the rehabilita- codal requirement of IS:13920-2016.32 The cross-sectional
tion work. area of the bar (Ash) of rectangular hoop reinforcement as
per IS:13920-201632 is calculated by the following equation
Cement mortar (CM)
[ fy ]
= maximum of 0 . 18 svh [A − 1]or ⌊0.05svh ⌋
Grade 43 OPC was used to prepare mortar mixture for fck
_ _ Ag fck
_
repair work.28 Cement to sand ratio was maintained at Ash
fy k
1:3 for mortar preparation, and water-cement ratio (w/c) (1)
was fixed at 0.35. For higher workability, a high-range
water-reducing admixture (HRWRA) at the rate of 1.5% by where sv is center-to-center spacing of stirrup; h is longer
weight of cement was used in the mixture. All these ingredi- length of rectangular stirrup measured up to outer face; Ag is
ents were mixed to prepare 50 cm (7.75 in.) mortar cubes to gross cross-sectional area of column; Ak is area of concrete
determine the compressive strength. Cubes were demolded confined in the core; fck is characteristic compressive strength
after 24 hours and cured under water. Compressive strength of concrete; and fy is 0.02% proof strength of transverse steel
was determined as per IS:1727-196726 and was found to be reinforcement bars.
25 MPa (3.63 ksi), 34 MPa (4.93 ksi), and 52 MPa (7.54 ksi) Adequate anchorage length was provided following
at the third, seventh, and 28th days, respectively. IS:13920-201632 to prevent pullout failure of beam rein-
forcement. In all specimens, both top and bottom bars of
Test specimens beams were bent at 90 degrees heading towards joint core.
In the present study, one-third scaled-down beam-column As per IS:13920-2016,32 the total anchorage length was
joint specimens were used. ACI 374.1-0529 allows scaling determined by Eq. (2)
down of specimens up to one-third of full-size specimens,
as the scaled-down specimens are large enough to represent Anchorage length =
load transfer mechanism along with all the real material (Ld + 10Ø – allowance for 90-degree bend) (2)
complexities and behavior of full-scaled specimens. The
size of columns and beam components were kept uniform where Ø is bar diameter; and Ld is development length given
at 150 x 120 mm (4.72 x 5.91 in.) for all specimens. Beam- by Eq. (3) as per IS:456-2000
column joints were designed on the assumption that the
point of contraflexure occurs at the beam midspan and _Ø σs
column midheight.30 Column length measured from one
Ld = 4 τ (3)
bd
contraflexure point to another was 1400 mm (55.12 in.). The
length of the beam measured from inner face of column to where Ø is bar diameter; σs is stress in longitudinal
the point of contraflexure in beam was 690 mm (27.17 in.). reinforcing bars; and τbd is design bond stress.
Beam-column subassemblies were designed following the Dimensions and reinforcement details of test specimens
strong-column, weak-beam concept.31 Based on the flexural are provided in Fig. 1.
@seismicisolation
@seismicisolation
Fig. 2—Schematic representation of cyclic test setup. (Note: 1 kN = 0.2248 kip; 1 mm = 0.0394 in.)
Test setup and experimental procedure experiment started by applying a displacement of magnitude
The experimental setup to test beam-column joints 1 mm (0.039 in.) to the beam-end. Displacement progressed
under cyclic loading is shown in Fig. 2. Quasi-static cyclic with an increment of 1 mm (0.039 in.) having 0.25 Hz
loading was applied to beam-column joints by means of a frequency until failure of the joint. The failure was marked
servo-controlled hydraulic actuator having a maximum by a drop in peak load by 25% as per ACI 374.1-05.29 Each
capacity of ±100 kN (22.48 kip) and peak displacement of displacement level is comprised of three cycles of push and
±90 mm (3.54 in.). The ends of the columns were mounted pull loading. Typical loading history is presented in Fig. 3.
horizontally on a roller steel support to simulate hinged- Similar loading protocol was also adopted by Mukherjee
boundary condition, and the beam was placed in the vertical and Joshi,35 Chidambaram and Agarwal,36 and Kheni et al.,37
plane, as shown in Fig. 2. The beam was hinged connected to among others. The displacement applied by the actuator
the actuator through swivel base assembly. A uniform axial was gauged with linear variable differential transformers
load of 10% gross capacity of the column was applied to (LVDTs) fitted in the actuator, and the corresponding load
the column with a hydraulic jack to replicate gravity loading was recorded by a coaxially mounted load cell. A data
condition.34 acquisition system collected all test data and stored it in the
In this study, the loading protocol was selected following computer connected to the actuator.
codal provision mentioned in ACI 374.1-05.29 The
@seismicisolation
@seismicisolation
Fig. 4—(a) Typical specimen after removal of damage concrete from joint region; and (b) typical beam-column joint after
rehabilitated by GM. @seismicisolation
@seismicisolation
Fig. 8—Load-displacement curve of test specimens. (Note: 1 kN = 0.2248 kip; 1 mm = 0.0394 in.)
Table 3—Peak load and joint shear stress of shear force; Lb and db are length and effective depth of beam,
specimens respectively; Lc is column length; and Dc is the overall
Specimens Peak load, kN τIS, MPa τjh, MPa τjh/τIS
column depth. The value of τjh is presented in Table 3. The
value of designed shear stress τIS is dependent on compres-
1C 9.10 5.37 2.31 0.43
sive strength of joint concrete. As in rehabilitated speci-
1GM 10.40 8.31 2.64 0.32 mens, joint core was strengthened with higher compressive
1CM 10.66 8.65 2.70 0.31 strength repair materials (fck = 48 MPa [6.96 ksi] for GM and
2C 10.27 5.37 2.60 0.49
fck = 52 MPa [7.54 ksi] for CM) compared to control spec-
imens (fck = 20 MPa [2.9 ksi]); the value of τIS for rehabili-
2GM 11.51 8.31 2.92 0.35
tated specimens increased significantly. On the other hand,
2CM 12.16 8.65 3.08 0.36 due to the same geometrical dimensions of all specimens, τjh
Note: 1 kN = 0.2248 kip; 1 MPa = 0.145 ksi.
is dependent primarily on P. The τjh of rehabilitated spec-
_ imens also increased with an increase in peak load owing
τIS = 1.2√fck
(4) to the use of high-strength GM and CM. A large increase
in τIS compared to τjh decreased the value of τjh/τIS so that
where fck is the characteristic strength of concrete. it is less than 1. Apparently, the ratio of τjh/τIS for all reha-
For the exterior beam-column joint, horizontal shear bilitated specimens indicates that specimens have adequate
force attained in the joint during loading is taken as follows shear strength capacity and enabled adjoining beam to reach
(Murty et al.47) its ultimate capacity. However, rehabilitated specimens
failed due to beam failure accompanied by joint failure, as
τjh = h ( − L )
_ P _ Lb _
Lb + 0.5Dc observed from Fig. 6. This may be because Eq. (5) does not
(5)
A core db c reflect the effect of cold joints formed due to rehabilitation
and the presence of yielded hoop reinforcement in the joint.
where P is the peak load of the test specimens; Ahcore is hori-
zontal cross-sectional area of the joint resisting horizontal
@seismicisolation
@seismicisolation
|F+i
_ |
− F−i
Ki = (6)
|δ+i− δ−i|
where F+i and F–i are maximum positive and negative load at
corresponding displacement δ+i and δ–I; respectively.
The stiffness calculated from Eq. (6) is plotted against
drift angle to obtain stiffness degradation curve, as shown in Fig. 10—Typical stress-strain graph of hoop reinforcement.
Fig. 9. Control specimens 1C and 2C displayed initial stiff- (Note: 1 MPa = 0.145 ksi.)
ness of 2.561 and 2.722 kN/mm (14.62 and 15.83 kip/in.),
and 2CM were able to resist diagonal tensile cracks forma-
respectively. GM and CM were only able to restore 79 to
tion at lower drift levels. Nevertheless, with increasing
80% and 65 to 67% of initial stiffness in rehabilitated speci-
drift, cracks at the cold joint interface expanded and caused
mens in comparison to control. This was due to the presence
development of new cracks in the joint. On the other hand,
of the cold joint interface between new repair mortars and
yielded hoop reinforcement also contributed to formation of
old concrete. When load was applied on rehabilitated speci-
crack in the joint. When yielded hoops were reloaded under
mens, the cold joint being a weak interface between concrete
reversed cyclic loading, yielded hoops yield at lower stress
and the repair material, a small amount of force opened
level and sustain lower ultimate strain. This is evident from
the interface. This led to earlier crack formation at lower
the stress-versus-strain graph of fresh hoops (tested before
displacement levels, as observed in the previous section
cyclic loading) and yielded hoops (taken out from control
where crack was first observed in rehabilitated specimens
specimens after cyclic loading), as shown in Fig. 10. Similar
at 1 mm (0.039 in.), irrespective of repair materials. On the
stress-strain curve of fresh and yielded reinforcing bar was
other hand, initial crack at 1C and 2C were observed at 2 and
also reported by Mukherjee and Jain.48 Yielding of yielded
3 mm (0.079 and 0.118 mm) displacement levels, respec-
hoops at a lower stress level produced loss of confinement
tively. Thus, earlier crack at lower displacement lowered
and deterioration of bond at lower beam-tip load compared
initial stiffness of all rehabilitated specimens. Higher initial
to control specimens. As a result, confinement of joint
stiffness of 1GM and 2GM than 1CM and 2CM was due to
hoops deteriorates significantly and becomes ineffective
higher bond strength of geopolymer with concrete substrate
against volumetric expansion of joint core. Such expansion
compared to OPC. As a result, a higher amount of force
44
of joint core cannot be resisted by the brittle nature of GM
is required to form first crack in 1GM and 2GM thereby
and CM, resulting in a large number of cracks in the joint.
leading to 20 to 21% increase in initial stiffness compared to
These factors contributed to higher crack development in
1CM and 2CM, respectively.
joint, thereby initiating higher stiffness degradation of 1GM,
Figure 9 reveals that 1GM, 2GM, 1CM, and 2CM exhib-
1CM, 2GM, and 2CM at higher drift.
ited lower initial stiffness, yet they displayed lower stiffness
degradation at lower drift level compared to controls. This
Pinching width ratio (PWR)
can be attributed to higher compressive strength of repair
The pinching width ratio (PWR) is a dimensionless param-
materials which increases the strength of diagonal strut,
eter developed by Mostofinejad and Akhlaghi49 to measure
as performance of diagonal strut depends on compressive
pinching associated with hysteresis loops under cyclic
strength of joint core material. Hence, 1GM, 2GM,@seismicisolation
1CM,
@seismicisolation
Fig. 13—Comparison of cumulative energy dissipation. (Note: 1 kNm = 0.0088 kip-in.; 1 mm = 0.0394 in.)
@seismicisolation
@seismicisolation
The torsional behavior of solid reinforced concrete (RC) members deterioration, particularly in North America, resulting in
reinforced with fiber-reinforced polymer (FRP) bars has been the costly repairs, rehabilitation, and a considerable reduction in
subject of several experimental studies. No experimental research, the service life span. Fiber-reinforced polymer (FRP) rein-
however, seems to have focused on RC box girders reinforced with forcing bars, on the other hand, can be used as an alternative
FRP bars under a pure torsional moment. This paper reports the
to steel reinforcement to combat corrosion, thereby lowering
results of an experimental investigation on the torsional strength
maintenance costs and extending service life. There are
and behavior of full-scale RC box girders reinforced with longi-
tudinal glass FRP (GFRP) bars. All specimens measured 380 mm numerous advantages to using FRP bars instead of steel
(15 in.) in height, 380 mm (15 in.) in width, 100 mm (4 in.) wall bars, such as higher tensile strength, weight one-quarter to
thickness, and 4000 mm (157.48 in.) in length. They were tested one-fifth lighter than steel, noncorroding nature, and longer
under pure torsional loading over a clear span of 2000 mm service life than steel (ACI 440.1R-15 [ACI Committee 440
(78.74 in.). The test specimens consisted of four RC box girders 2015]). The last two decades have seen a rapid increase in
with longitudinal GFRP bars and one RC box girder with longitu- the use of FRP reinforcing bars in many applications, such
dinal steel bars as a reference. All the specimens were constructed as bridges, piles, parking garages, marine structures, water
without web reinforcement to study the contribution of the longi- tanks, and tunnels (El-Salakawy et al. 2004; Eladawy et al.
tudinal reinforcement to torsional strength. The test variables 2019; Mohamed and Benmokrane 2014; Mousa et al. 2018;
included the longitudinal reinforcement ratio (ranging between
Mohamed et al. 2020).
1.10 and 2.74%) and the type of longitudinal reinforcement (GFRP
Due to the scarcity of experimental studies on torsion
or steel). The test results indicate that increasing the GFRP longi-
tudinal reinforcement ratio increased the torsional strength after members reinforced with FRP bars, the FRP-RC design
the initiation of the first diagonal crack, especially for specimens guidelines (ACI 440.1R-15; JSCE 1997) do not include
with a high reinforcement ratio. In addition, theoretical torsional any provisions pertaining to torsion. In addition, the torsion
moment-twist curves were developed and gave predictions consis- provisions in FRP-RC design codes (CSA S806-12 [2017];
tent with the experimental test results. Lastly, the ultimate torsional CSA S6 [2019]; AASHTO 2018) are mainly modifications
strength of the GFRP-RC box girders without web reinforcement of steel-RC codes. These modifications include the mechan-
was estimated with the CSA S806-12 (R2017) design equation with ical properties of FRP reinforcement (modulus of elasticity
a modification related to the GFRP tensile strain limit. and stirrup bending strength). The design of torsion in
FRP-RC codes (CSA S806-12; CSA S6:19; AASHTO 2018)
Keywords: cracking pattern and modes of failure; design codes; glass
fiber-reinforced polymer (GFRP) bars; reinforced concrete (RC) box is based on a thin-walled tube space-truss analogy. Once RC
girders; reinforcement ratio; theoretical modeling; torsional crack width members have cracked under torsion, the torsional strength
and toughness; torsional loading; torsional moment-twist response; is mainly provided by closed stirrups and longitudinal bars.
torsional strength. CSA S806-12 provides an equation to estimate the torsional
strength based on the contribution of the GFRP longitudinal
INTRODUCTION reinforcement; this equation is a function of the total area
Reinforced concrete (RC) box girders can be found in and tensile strength of the GFRP longitudinal bars. The
many civil engineering applications. The most notable are experimental program of the current study was intended to
cable-supported bridges, pedestrian bridges, curved bridges, evaluate the accuracy of the CSA S806-12 design equation.
and modern elevated structures for light rail transport (Rahal Valuable research work has been carried out in the last
and Collins 1995). There are economic and structural bene- decade to investigate the torsional behavior of the FRP-RC
fits to using box girders, such as low self-weight of struc- members with a solid cross section (Deifalla et al. 2014;
tures and higher torsional stiffness (Jeng et al. 2013). Torsion Hadhood et al. 2020; Mohamed and Benmokrane 2015;
in bridges could occur as a consequence of the geometric Shehab et al. 2009; Zhou et al. 2017). These studies
complexities of horizontally curved bridges and/or the large
eccentricity of gravity loads. Such structures are usually ACI Structural Journal, V. 120, No. 5, September 2023.
exposed to aggressive and/or harsh marine environments, MS No. S-2022-137.R3, doi: 10.14359/51738836, received November 28, 2022, and
reviewed under Institute publication policies. Copyright © 2023, American Concrete
leading to corrosion of the steel reinforcement. This type of Institute. All rights reserved, including the making of copies unless permission is
obtained from the copyright proprietors. Pertinent discussion including author’s
corrosion in concrete structures is a prime cause of structural closure, if any, will be published ten months from this journal’s date if the discussion
is received within four months of the paper’s print publication.
@seismicisolation
@seismicisolation
Note: ρL = AL/Ac, where AL is total area of longitudinal reinforcement and Ac is area of outer perimeter of concrete cross section. 1 MPa = 0.145 ksi.
shown in Fig. 1(b). Once the concrete was cast, the spec- box girders reinforced with longitudinal GFRP bars but
imens were covered with wet burlap and plastic sheets to without transverse reinforcement. A total of five large-scale
prevent moisture loss and were cured for 7 days. All speci- RC box girders—including four reinforced with longitudinal
mens were stored in the laboratory until the day of testing. GFRP bars and one reinforced with longitudinal steel bars
Figure 1 shows the casting, curing, and storage process of as a control specimen—were tested under pure torsional
the test specimens. loading. The test matrix was designed to study the effect of
longitudinal reinforcement on the torsional behavior of RC
Specimen design box girders reinforced with GFRP bars. The test parame-
The experimental program of this study was designed to ters included the type of longitudinal reinforcement (GFRP
provide experimental data on the torsional behavior@seismicisolation
of RC versus steel) and the longitudinal reinforcement ratio. The
@seismicisolation
Note: Tcr is torsional moment corresponding to first diagonal crack; Tu is maximum torsional moment; Kun(exp) is uncracked torsional stiffness (Tcr/φcr); and Kcr(exp) is cracked
torsional stiffness ([Tu – Tcr]/[φu – φcr]). 1 kN∙m = 0.7376 kip∙ft; 1 kN∙m2 = 2.42 kip∙ft2; 1 kN = 0.225 kip; 1 rad/m = 0.305 rad/ft.
Fig. 6—Effect of test parameters on torsional moment-twist behavior. (Note: 1 kN∙m = 0.7376 kip∙ft; 1 rad/m = 0.305 rad/ft.)
Effect of type of longitudinal reinforcement on
torque-twist response
Figure 6(b) shows the effect of longitudinal reinforce-
ment type on the torsional moment-twist behavior. Spec-
imen BGW-12 and its counterpart steel-reinforced specimen
(BSW-12) were designed to have the same longitudinal
reinforcement ratio (1.66%). The figure indicates that the
cracking strength and the pre-cracking torsional stiffness
were almost the same for both specimens. Table 3 indicates
that the torsional cracking strength of BGW-12 and BSW-12
was 36.55 and 36.88 kN∙m (26.97 and 27.22 kip∙ft), respec-
tively. After cracking occurred, the specimens’ torsional
performance was highly dependent on the axial stiffness
of the longitudinal reinforcing bars, which is a function of
the area A and modulus of elasticity E of the longitudinal
reinforcement. Figure 6(b) indicates that the steel-reinforced
specimen (BSW-12) experienced less of a drop in torsional
strength than its counterpart GFRP-reinforced specimen
(BGW-12) after the first diagonal crack appeared. In contrast,
BSW-12 achieved a slight increase in torsional strength of
Fig. 7—Ratio of post-cracking to cracking strength versus
7% after the first crack, followed by a reduction of 10% in
reinforcement ratio.
ultimate torsional strength. The reduction in the torsional
strength of the GFRP-RC box girder (BGW-12) was almost (200 GPa [29,000 ksi]) to that of the GFRP bars (62.5 GPa
3.29 times that of its steel counterpart (BSW-12). This value [9063 ksi]). Afterward, the two specimens showed degrada-
is similar to the ratio of the modulus of elasticity of the@seismicisolation
steel tion in stiffness and strength as the angle of twist increased.
@seismicisolation
_T
Kun(Theo)
= φ (1)
Fig. 10—Torsional moment-crack width. (Note: 1 kN∙m = where T is the applied torque to the girder; φ is the twist
0.7376 kip∙ft; 1 mm = 0.0394 in.) per unit length; and Kun(Theo) is the theoretical pre-cracking
the longitudinal reinforcement ratio from 1.10% in BGW-8 torsional stiffness (Kun(Theo) = RGC); R is the reduction factor
to 2.74% in BGW-20 increased the torsional toughness by (R ≈ 0.60); G is the shear modulus of concrete and is equal to
approximately 325%. The steel-reinforced specimen exhib- Ec/(2(1 + υ)); Ec is the modulus of elasticity of the concrete;
ited torsional toughness approximately 12% higher than its υ is the Poisson’s ratio; and C is the Saint-Venant torsional
counterpart with GFRP reinforcement. It can be concluded constant.
that the torsional toughness significantly increased as the The reduction factor (R) considers the stiffness loss
longitudinal reinforcement ratio increased. observed in the laboratory tests, which was almost 20 to
40% of the Kun(Theo) value. It was already reported by Leon-
THEORETICAL INVESTIGATION hardt and Schelling (1974) and adopted by Bernardo and
Theoretical prediction of torsional Lopes (2008). The torsional reinforcement can be neglected
moment-twist curve in calculations of the Kun(Theo) (Bernardo and Lopes 2008;
In RC members subjected to the pure torsional moment, Hsu 1973). For rectangular hollow sections, the Saint-
the torque is generally resisted by longitudinal rein- Venant torsional constant (C) is calculated as
forcement, transverse reinforcement, and concrete diag-
onal struts. Therefore, the torsional moment-twist curve _
4A2h
C = u (2)
comprises three distinct zones (Bernardo and Lopes 2008).
The first zone is characterized by the linear-elastic behavior where A is the area enclosed by the centerline of the wall
of materials and represents the torsional behavior until the of the hollow section; h is the wall thickness of the hollow
first diagonal crack occurs. The cracking strength in this section; and u is the perimeter of the centerline of the wall of
zone can be predicted theoretically with the skew-bending the hollow section.
theory, Bredt’s thin-tube theory, and the theory of elasticity. Cracking torsional strength—The cracking torsional
The slope of the curve in the first zone represents uncracked strength (Tcr) can be theoretically predicted with three
torsional stiffness, which can be predicted with Saint- different theories (the skew-bending theory, Bredt’s thin-
Venant stiffness. The second zone corresponds to the torsional tube theory, and the theory of elasticity), in addition to three
behavior in a cracked state with linear-elastic behavior of codes (CSA A23.3:19, CSA S6:19, and AASHTO 2018),
the materials, and the slope of the curve represents the post- which are briefly described as follows.
cracking torsional stiffness. The space-truss analogy with an Skew-bending theory—Hsu (1968) derived an equation
inclination angle of diagonal concrete struts of 45 degrees for the nominal torsional strength of solid concrete members
can theoretically predict the torsional behavior in the second based on the bending mechanism of torsional failure. This
zone. The third zone represents the nonlinear behavior of theory contends that the failure of a torsional member is
the materials with a softening effect because at least one of reached when the tensile stress induced by a 45-degree
the reinforcement types (longitudinal or transverse) reaches bending component of torque on the wider face reaches
the ultimate strength or because of the nonlinearity of the a reduced modulus of rupture of concrete. The torsional
diagonal concrete struts. In contrast, based on the experi- cracking strength is given by
mental test results from this study, the ultimate torque for
the specimens reinforced with relatively lower GFRP longi- ____
fc′ ( x )for h ≤ 4
3 _
4h _x
tudinal reinforcement ratios (1.10 and 1.66%) coincided
T
cr = 6y(
x
2
+ 10)
√ (3)
with the cracking torque. Therefore, the torsional moment-
twist curve for these specimens terminated in the first where x and y are the shorter and longer dimensions of
zone. In contrast, the specimens reinforced with relatively the cross section (in.), respectively; and fc′ is the concrete
@seismicisolation
@seismicisolation
@seismicisolation
@seismicisolation
This paper presents an experimental study on the anchorage To ensure full tensile yield capacity of the tie elements
behavior of drilled shaft reinforcement subjected to tension in comprising a strut-and-tie model, sufficient development
drilled shaft footings loaded under combined axial force and length (anchorage length) should be provided beyond the
uniaxial bending moment. Four large-scale tests were conducted point at which the tie meets the strut anchoring it. Current
on drilled shaft footing specimens employing an equivalent loading
STM provisions1,2 establish the critical section at which
condition introducing tension in the drilled shaft reinforcement.
the bar should be developed based on the geometry of the
Three different anchorage details were tested: straight bars,
hooked bars, and headed bars. The drilled shaft reinforcement extended nodal zone bounded by the edge of the strut in a
was capable of developing its full yield strength in tension in all singular node (for example, a node with a well-bounded
the tests, regardless of the anchorage detail. The tensile stresses bearing face). Specifically, the critical section is defined at
in drilled shaft bars were primarily developed in the region of the the point where the centroid of the tie element intersects the
embedment length closest to the interface between the drilled shaft extended nodal zone, as shown in Fig. 2(a). For smeared
and the footing, while negligible stress and slip were measured in nodes, it is not possible to identify the critical section for
the vicinity of the unloaded end of the bars for all three anchorage tension development based on the extended nodal zone
details. Finally, a critical section is also proposed in this study to owing to the absence of a well-defined nodal geometry. An
perform the anchorage check for the drilled shaft reinforcement alternative method was proposed by Yi et al.3 to determine
in drilled shaft footings designed with the strut-and-tie method.
the critical section of the column ties in Fig. 2(b), corre-
The definition of the critical section provides a safe estimate of the
sponding to a drilled shaft footing subjected to combined
available development length of the drilled shaft bars according to
the findings of the experimental program. axial force and a moderate uniaxial bending moment. This
method assumes a large compression field bounded by the
Keywords: bar anchorage; bond; development length; drilled shaft footing; idealized diagonal struts near the column tie element of
drilled shaft reinforcement; large-scale tests; strut-and-tie method (STM). the 3-D strut-and-tie model, and the stress field is used to
define the critical section for the column tie element at the
INTRODUCTION smeared node. The proposed compression field is considered
A drilled shaft footing is a component of a pier founda- to perform a similar role to the extended nodal zone in deter-
tion transferring forces from the pier to a group of drilled mining the critical section, as shown in Fig. 2(b). The safety
shafts. Drilled shaft footings are generally considered to of this critical section criterion was verified experimentally.3
behave as D-regions due to their reduced shear span-depth However, the same criterion cannot be used for the case
ratios and are therefore recommended to be designed using represented by the 3-D strut-and-tie model in Fig. 1(b) due
the strut-and-tie method (STM). Three-dimensional (3-D) to its different configuration of struts and ties. Therefore, it
strut-and-tie models are required to represent internal forces remained unclear how to perform the anchorage check for
in drilled shaft footings supported on a grid of drilled shafts. the column and drilled shaft tie elements in the 3-D strut-
The configuration of the 3-D strut-and-tie model depends on and-tie model for drilled shaft footings under combined
the loading combination applied through the pier. Figure 1 axial force and a large uniaxial bending moment.
presents two strut-and-tie models for a drilled shaft footing Most previous experimental studies4-10 on drilled shaft
subjected to uniaxial flexure and vertical compression footings were conducted for the simplest loading condition:
loading. The model in Fig. 1(a) corresponds to a column a pier subjected to uniaxial compression. To the authors’
subjected to combined axial force and a moderate axial knowledge, a limited number of experimental studies3,11
bending moment, resulting in tension at one face of the have been conducted on drilled shaft footings subjected to
column and nonuniform compression in drilled shafts, while uniaxial compression combined with the moderate flexural
the model in Fig. 1(b) corresponds to a case of combined moment, inducing tension at the column reinforcement
axial force and a large uniaxial bending moment, resulting and nonuniform compression in drilled shafts. However,
in tension not only at one face of the column but also at no experimental studies have been reported to investigate
two of four drilled shafts. Both column and drilled shaft tie
elements in the model are connected with smeared nodes, in
which the nodal geometry cannot be clearly established due ACI Structural Journal, V. 120, No. 5, September 2023.
MS No. S-2022-257.R1, doi: 10.14359/51738837, received January 21, 2023, and
to the absence of well-defined geometrical constraints (for reviewed under Institute publication policies. Copyright © 2023, American Concrete
Institute. All rights reserved, including the making of copies unless permission is
example, bearing area or strut width limitation). obtained from the copyright proprietors. Pertinent discussion including author’s
closure, if any, will be published ten months from this journal’s date if the discussion
is received within four months of the paper’s print publication.
@seismicisolation
@seismicisolation
Fig. 2—Critical section for anchorage verification of column tie elements in drilled shaft footings subjected to combined
axial force and moderate uniaxial bending moment.
the response of drilled shaft footings subjected to flexural drilled shaft footings subjected to combined axial force and
moments large enough to induce tension in the drilled shafts. a large uniaxial bending moment, resulting in tension in two
Therefore, due to this limited experimental database, of four drilled shafts. They proposed anchorage details for
studies12-20 on STM-based analytical models for drilled drilled shaft reinforcement under tension based on conser-
shaft footings that correspond only to uniaxial compres- vative assumptions and best construction practices given
sion have been conducted. A few studies21,22 proposed the lack of in-depth experimental or computational studies
STM-based design recommendations and examples for on this problem. The design example provided by Widianto
@seismicisolation
@seismicisolation
time of design. Therefore, this study defined the available drilled shaft of the specimens. The reinforcement details are
embedment length (lp,a), an extended embedment length of summarized in Table 2 and illustrated in Fig. 7.
drilled shaft reinforcement that extends up to the elevation
of the top mat reinforcement, to determine an adequate size Material properties
for the drilled shaft reinforcement. Given the geometry of Two batches of concrete were used to fabricate each
the footing specimens, the size of the drilled shaft bars was footing specimen. The mixture proportions were identical for
determined to provide a realistic ratio between the available all batches of concrete with a design compressive strength
embedment length (lp,a) and bar diameter (db,s) according to of 3.6 ksi (24.8 MPa). ASTM A706/A706M-1624 Grade 60
the footing database. The number of bars was then deter- reinforcing bars were used for the drilled shaft reinforce-
mined to have a similar drilled shaft reinforcement ratio ment. ASTM A706/A706M24 bars were chosen to be able to
(ρs) as that observed in the footing database. To reproduce weld them to a steel plate needed to apply the tensile load,
a similar available embedment length ratio (lp,a/db,s = 43.1) as will be discussed later. The rest of the reinforcement in
and drilled shaft reinforcement ratio (ρs = 1.05%) averaged the test specimens were ASTM A615/A615M-2025 Grade 60
from the footing database, five No. 6 reinforcing bars (ρs = reinforcing bars. The mechanical properties of the concrete
1.09%) were provided and extended up to the top mat rein- and reinforcement were determined from the average results
forcement (lp,a/db,s = 46.3) for each 16 in. (406 mm) diameter of materials testing on a minimum of three samples of each
@seismicisolation
@seismicisolation
loads corresponding to the yielding of all the drilled shaft shaft plate on the specimen measured unstable data after
reinforcing bars. As shown in Table 4, these load values are drilled shaft reinforcement yielding due to the conical cracks
within 10% of the average value of 347 kip (1.5 MN) of the formed around the drilled shaft plate. Therefore, the west-
four tests. The maximum stresses of the drilled shaft rein- side plate was selected for VII-TK-ST (straight drilled shaft
forcement were computed based on strain gauge data at the reinforcement), and the east-side plate was selected for the
maximum applied load of 400 kip (1.8 MN). For all spec- other specimens to represent the behavior of the drilled shaft
imens, the maximum stresses were commonly developed reinforcement embedded in the footing. The load in Fig. 10
at the measurement location closest to the top surface of corresponds to the load measured by the load cell in the same
the footing (2 in. [51 mm] below the top surface), and they side’s actuator. Figure 10 also presents the results in terms
exceeded 90% of the tensile strength of the bars in all cases. of the average stress level on the drilled shaft reinforcement
The theoretical capacity of the specimens according to the by dividing the applied load by the total area of reinforce-
3-D STM would correspond to the yielding of the drilled ment connected to one drilled shaft plate. As shown, all four
shaft reinforcement at a load of 300 kip (1.3 MN). The differ- tests presented a consistent bar stress-elongation response.
ence between the theoretical yield load and the experimental The average drilled shaft reinforcement stress exceeded the
values of the load for first yielding and complete yielding of yield stress at a load of approximately 300 kip (1.3 MN) and
the bars, indicated in Table 4, can be explained by the fact reached a tensile stress of approximately 90 ksi (621 MPa)
that the distribution of tensile forces among the bars was not at a load of 400 kip (1.8 MN).
perfectly uniform, as will be discussed later based on the The reinforcement slip measured at the bottom tip of the
analysis of bar stress data. drilled shaft reinforcement showed a trend to increase as
The overall test response is also analyzed in terms of the load increased; however, the values were very small (of the
load applied to one drilled shaft versus the average elon- order of 10–5 in. [10–4 mm]) for all anchorage types. Hence,
gation of its drilled shaft bars in Fig. 10. To this end, the the slip was negligible.
measurements obtained from the linear potentiometers
attached on the drilled shaft plates were averaged to calcu- Visual observations
late the average displacement of the plate with respect to During each test, all faces of the tested span of the footing
the footing, which represents the average elongation (or slip) were visually inspected. During testing of VII-TK-ST, a
of the drilled shaft reinforcement at the top surface of the series of horizontal cracks formed at midheight of the north
footing. Some linear potentiometers installed at one drilled face of the specimen as the load increased. Those cracks are
@seismicisolation
@seismicisolation
Fig. 16—Proposed critical sections for column and drilled shaft reinforcement in drilled shaft footings under combined axial
force and large uniaxial bending moment.
@seismicisolation
@seismicisolation
This paper is an attempt at a better understanding of design limitations to non-seismic zones and structures not requiring
provisions of ACI CODE-440.11-22, building code for the design fire resistance. The addition of provisions for compressive
of glass fiber-reinforced polymer (GFRP)-reinforced concrete members is a critical development for practitioners inter-
(RC) columns. Sway and a non-sway column examples origi- ested in nonmetallic reinforcement as it allows designing
nally designed with steel reinforcement were redesigned using
and construction of a building entirely with GFRP-RC.
GFRP longitudinal bars and ties as per provisions of ACI CODE-
Though ACI CODE-440.11-22 permits the design of
440.11-22 to analyze the effect of changing reinforcement type.
Columns were designed with both low-modulus (Ef = 6500 ksi), columns using GFRP bars, due to their lower reliability, the
and high-modulus (Ef = 8700 ksi) GFRP bars. A parametric study minimum compressive strength properties of GFRP bars
was carried out by varying the concrete compressive strength, are not specified in ASTM D7957.12 As stipulated in ACI
the cross-section aspect ratio, and the resultant load eccentricity. CODE-440.11-22, in pure compression, their presence can
GFRP-RC columns require larger cross-section dimensions and be treated as having the same stiffness and strength as those
more reinforcement area than steel-RC columns irrespective of of the surrounding concrete. However, in the presence of
the GFRP elastic modulus when subjected to the same demand. moment, GFRP reinforcement may effectively contribute
The concrete strength has a significant effect on the dimensions to the column capacity. Therefore, this study is carried out
of GFRP-RC columns, and rectangular sections were found to be to show the implications of current Code provisions on the
more efficient than square sections with the same gross concrete
design of GFRP-RC compressive members.
area in the presence of moment. GFRP-RC columns subject to high
eccentricity loads take advantage of GFRP tensile properties and,
thus, are more efficient. RESEARCH SIGNIFICANCE
The recently published ACI CODE-440.11-2211 allows
Keywords: building code; concrete columns; eccentricity; glass fiber-rein- the design of columns with GFRP reinforcement. Due to
forced polymer (GFRP) reinforcement. remaining knowledge gaps in the behavior of GFRP-RC
columns, some Code provisions were only analytically
INTRODUCTION developed and verified by incorporating differences in material
Fiber-reinforced polymer (FRP) bars, being a competitive properties with steel-RC. This study shows the implications of
option for reinforced concrete (RC) members in aggressive Code provisions and highlights the areas for further research.
environments, were not allowed in compression members in
the previous editions of the ACI 440 Guide.1 The primary METHODOLOGY
reason for this exclusion was a lack of information regarding In this study, a column part of a sway frame from the ACI
the behavior of FRP-RC members subjected to compressive Reinforced Concrete Design Handbook,13 a Companion to
loads. However, researchers have been actively investigating ACI 318-19,14 is selected and redesigned using GFRP rein-
the behavior of glass FRP (GFRP)-RC columns during the forcement. This column is part of an interior, continuous,
last decade and have found GFRP-RC columns to be permis- six-bay frame, and built integrally with a 7 in. (178 mm)
sible structural elements. In fact, several experimental deep slab, as shown in Fig. 1(a). The constituent mate-
studies investigated the effect of the compressive behavior rials selected for column design are listed in Table 1.
of longitudinal GFRP bars by testing RC columns2-5 with The concrete strength fc′ is 5000 psi (35 MPa), while the
an overall positive assessment of their feasibility. Jawaheri GFRP type is compliant with material specification based
Zadeh and Nanni6 provided information on flexural stiff- on ASTM D7957.12 For the non-sway case, a column from
ness in frame analysis for GFRP-RC that resulted in close a frame part of an industrial building was taken from the
correspondence to limits proposed by Bischoff.7 Similarly, textbook by Wight and Macgregor,15 as shown in Fig. 1(b).
Hadhood et al., among other researchers, proposed a 1%
8 9
This is a laterally braced column with a beam on one side, as
minimum reinforcement necessary to maintain section shown in Fig. 1(b). The concrete strength fc′ for this column
integrity to achieve a nominal capacity of columns. Khor- was 4000 psi (28 MPa) (as given in the textbook), while
ramian and Sadeghian10 performed structural tests validating
the performance of GFRP-RC columns with reinforcement ACI Structural Journal, V. 120, No. 5, September 2023.
MS No. S-2022-307.R3, doi: 10.14359/51738838, received April 10, 2023, and
ratios as high as 5.3%. Given, these significant advances in reviewed under Institute publication policies. Copyright © 2023, American Concrete
research over the past decade, the new ACI CODE-440.11-22 Institute. All rights reserved, including the making of copies unless permission is
obtained from the copyright proprietors. Pertinent discussion including author’s
Building Code permits the use of GFRP-RC columns with
11
closure, if any, will be published ten months from this journal’s date if the discussion
is received within four months of the paper’s print publication.
@seismicisolation
@seismicisolation
ΦSn ≥ U (1)
(0.01Ef)
ffu
moments.
Design strength = m
in (b) The stability index does not exceed 0.05. The stability
index for a given story, Q, shall be calculated as shown
Code Section 10.6.1 specifies a minimum reinforcement
of 1% of the gross concrete area (Ag) to provide resistance to ∑ Pu δo
_
Q = (3)
bending and possibly concrete creep. Similarly, maximum Vus
lc
reinforcement of 8% is specified to avoid congestion of
reinforcing bars and to ensure that concrete can be properly where ∑Pu is total factored vertical load; Vus is horizontal
consolidated. story shear; δo is first-order relative lateral displacement
The minimum number of bars is indicated by the Code between the top and bottom of that story; and lc is height of
Section 10.7.3, given as the column from the center to center of the joints.
GFRP-RC columns are more susceptible to the slender-
Minimum number of longitudinal bars = ness effects than steel-RC due to the lower stiffness of GFRP
⎧4 Rectangular or circular ties reinforcement compared to steel bars; therefore, more strict
limits are imposed when checking slenderness effects for
⎨
⎪
3
Triangular ties
GFRP-RC columns. Slenderness effects can be neglected in
⎩6
⎪
Enclosed by spirals both sway and non-sway frames if the following conditions
of Code Sections 6.2.5.1(a) or 6.2.5.1(b) are satisfied, given
For longitudinal reinforcement, the minimum clear herein as Eq. (4a), (4b), and (4c)
spacing between bars is specified in the Code Section 25.2.3 (a) For columns not braced against sidesway
as follows
k lu
_
1.5 in.(38 mm)
r ≤ 17 (4a)
( 4 / 3dagg
)
Minimum spacing between bars = max 1.5db
(b) For columns braced against sidesway
where Bdns is the ratio of maximum factored sustained axial (Mu)top, kip-ft 145 38
load to maximum factored axial load; and If is the moment (Mu)bott, kip-ft 197 94.4
of inertia of the GFRP bars about the centroid of the cross
Note: 1 kip = 4.44 kN; 1 kip-ft = 1.35 kN-m.
section.
Code Section 6.6.4.5 implies amplifying M2 for the effects laterally unbraced; therefore, slenderness effects were
of member curvature in a non-sway frame given as checked as per ACI CODE-440.11-22, Eq. (6.2.5.1a) (that
is, klu/r < 17). The unsupported length of the column, as
Mc = δM2 (10) shown in Fig. 1(a), is 15.5 ft (4.72 m) and its cross-section
dimensions were calculated using relations provided in this
where Mc is factored moment amplified for the effects paper equal to 26 x 26 in. (660 x 660 mm). The effective
of member curvature; and δ is magnification factor for length factor was calculated using alignment charts given
non-sway frames as given in Code Section 6.6.4.5.2 in Fig. R6.2.5.1 in ACI CODE-440.11-22 (that is, Fig. 2),
which depend on the relative stiffnesses of columns to the
_
C m floor members at column top and bottom joints. In this
δ = Pu ≥ 1.0 (11)
_ design example, the column frames into beams at the top
1 − 0.75P
c joint, whereas it frames at bottom in a two-way slab. It was
assumed that the columns in the stories above and below
where Cm (factor relating actual moment diagram to an equiv- had the same cross-section dimensions. The gross moment
alent uniform moment diagram) shall be in accordance with of inertia of the column was equal to 38,080 in.4 (15.8 × 109
6.6.4.5.3a and 6.6.4.5.3b, given herein as Eq. (12a) and (12b): mm4), and the effective moment of inertia calculated as per
(a) For columns without transverse loads applied between Table 3 was equal to 15,232 in.4 (6.3 × 109 mm4).
supports As stated in ACI CODE-440.11-22, Section R6.6.3.1.1, it
is sufficiently accurate to take the gross moment of inertia of
_M1 a T-beam equal to twice that of its web. Using this approach,
Cm = 0.6 − 0.4 M (12a)
2
the moment of inertia of T-beams framing into the column
at the top joint was calculated equal to 81,000 in.4 (33.7 ×
(b) For columns with transverse loads applied between
109 mm4), and the effective moment of inertia as per Table 3
supports
was equal to 12,150 in.4 (5 × 109 mm4). Similarly, the
moment of inertia at the lower joint was calculated for the
Cm = 1.0 (12b)
slab framing into the column. The width of the slab in the
transverse direction was considered equal to 14 ft (4.3 m)
EXAMPLES OF COLUMN DESIGN AND
and its thickness equal to 7 in. (178 mm), as given in the ACI
DISCUSSION
318-19 Design Handbook.13 Its gross moment of inertia was
The required strength for the two columns subjected to
calculated equal to 4802 in.4 (2 × 109 mm4), which reduced to
lateral and gravity loads was checked using the factored
720 in.4 (0.3 × 109 mm4) when calculating effective moment
load combinations in Chapter 5 and analysis procedures in
of inertia as per Table 3.
Chapter 6 of the ACI CODE-440.11. The calculated values
The relative stiffness at the top and bottom joints was calcu-
of axial load, moment, and shear demands used in this study
lated as per Eq. (5) of this paper (as given in ACI CODE-
are given in Table 4 as originally available from the sources of
440.11-22), which were found equal to 3 for top and 30 for
the steel-RC cases.12,15 It should be noted that for simplicity, a
bottom joints, respectively. Using relative stiffness factors in
single combination of ultimate axial load and moment (that is,
alignment charts given in Fig. 2, the effective length factor
Pu and Mu) for each of the two columns was adopted, whereas
was calculated equal to 2.8 and radius of gyration equal to
in practice, the demand of several combinations of loads and
7.5 from ACI CODE-440.11-22, Section 6.2.5.2. The values
moments must be satisfied. Also, the design was carried out
of effective length factor, unsupported length of column,
by keeping the reinforcement amount as close as possible to
and radius of gyration were used in Eq. (6.2.5.1a). It was
the minimum requirements of the Code (that is, 1%Ag).
observed that slenderness effects cannot be neglected; hence
the column should be designed by considering the second-
Sway column using low- and high-modulus GFRP
order effects.
bars
The second step after slenderness is to investigate if
The column was designed using the Code-referenced
the column should be analyzed as sway or non-sway. ACI
low-modulus (Ef = 6500 ksi [44,815 MPa]) GFRP bars, with
CODE-440.11-22, Section 6.6.4.3 states that a column can
a concrete strength of 5000 psi (35 MPa). The columns were
be analyzed as part of a non-sway frame if: (a) column end
designed as per procedure provided in ACI CODE-440.11-
moments due to second-order effects do not exceed 5%
22, Fig. R6.2.5.3. The column in this design example is
of the first-order end moments; and (b) the stability index
@seismicisolation
@seismicisolation
compressive force corresponding to zero eccentricity (Mn = the capacity (9%) for the column subjected to same demand
0 and x = +¥), which is the uppermost point in the interaction (that is, as a low-modulus GFRP-RC column given in
diagram. Point B′ on the design domain limit represents the Table 5) was noticed. The column cross section, reinforce-
case of maximum compressive force usable in design. The ment details, and interaction diagrams developed for both
two modes of failure (tension and compression-controlled low- and high-modulus GFRP-RC columns are shown in
modes) are separated by the “balanced failure” shown by Fig. 3.
point C′, representing FRP rupture (note: the guaranteed
strength of GFRP bars is replaced by 0.01Ef, as specified in Non-sway column using low- and high-modulus
the Code) and concrete crushing simultaneously (x = xb). If GFRP bars
the neutral axis shifts beyond xb, the failure mode shifts from The column example taken from the textbook by Wight
compression to tension. The lowermost point in the interac- and Macgregor15 was redesigned with GFRP reinforcement,
tion diagram (E′) corresponds to maximum tensile force (Mn considering it as part of a non-sway frame as the stability
= 0 and x = –¥), and maximum strain in the reinforcement. index (Q = 0.04) and magnified moment were within the
Any combination of ultimate axial load and moment (that limits stated in Section 6.6.4.3 for a non-sway frame. The
is, Pu-Mu, shown by a black dot) laying within the interac- magnification factor calculated by analysis was 0.7; hence,
tion curve represents safe (and outside, an unsafe) column a minimum magnification factor of 1.0 was used to calcu-
design. late the magnified moment. The P-M diagram developed for
The GFRP-RC column subjected to same ultimate loads non-sway columns showing the nominal and design capacity
required larger cross section compared to steel-RC as the curves is shown in Fig. 4. The (Pn, Mn) curve shows the
axial strength of GFRP reinforcement is not considered capacity before and (ΦPn, ΦMn) after the application of the
in resistance calculations and is replaced with equal area strength reduction factors. The column failed with dimen-
of concrete. Also, the higher magnification factor resulted sions of the steel-RC section (14 x 14 in. [356 x 356 mm])
in a very large, magnified moment, as given in Table 5, and reinforcement consisting of four No. 8 (M25) bars, as
together with the limits on the maximum strength of GFRP shown in Fig. 4(a), where the demand (Pu-Mu) shown by the
bars, the GFRP-RC column required bigger cross-sectional black dot lies outside the design capacity curve. Therefore,
dimensions than a steel-RC column. For example, a column the cross-section dimensions were increased to 18 x 18 in.
designed with GFRP reinforcement failed with dimensions (460 x 460 mm) and reinforcement consisting of four No. 9
similar to that of the steel-RC (24 x 24 in. [610 x 610 mm]); (M29) bars to satisfy the demand. As shown in Fig. 4(b),
therefore, the column size was increased to 26 x 26 in. (660 x the demand is within the design capacity curve, representing
660 mm) to augment its load-carrying capacity to exceed a safe column. It was observed that GFRP-RC columns
the demand. The minimum reinforcement depends on the require bigger cross sections and reinforcement areas than
gross area of the cross section (Afmin = 0.01Ag); therefore, steel-RC. The effective length factor, moment magnification
a GFRP-RC column required more reinforcement area than factor, cross-sectional dimensions, and required reinforce-
steel-RC. For example, a column designed with GFRP-RC ment are shown in Table 6.
required eight No. 9 (M29) bars, whereas that with steel-RC Similar to the case of sway frames, an effort was made
required eight No. 8 (M25) bars. The values of effective to investigate the impact of using high-modulus bars
length factor, moment magnification factor, cross-sectional on non-sway column design. The high-modulus bars
area, and longitudinal reinforcement for a column in sway showed a 12% increase in the column capacity compared to
frame are provided in Table 5, which includes values for low-modulus GFRP-RC column subjected to same demand
steel-RC taken from ACI 318-19 Design Handbook.13 (that is, the low-modulus GFRP-RC column given in Table 6).
This study also investigated the effect of high-modulus The column cross section, reinforcement details, and interac-
(Ef = 8700 ksi [60,000 MPa]) GFRP reinforcement on column tion diagrams developed for non-sway column are shown in
design. The compressive strength of GFRP reinforcement is Fig. 4.
not considered in resistance calculations when the GFRP
is in compression (that is, the area of GFRP replaced with PARAMETRIC STUDY
concrete); hence, the column dimensions remained same as A parametric study was carried out by varying concrete
with low-modulus GFRP bars (26 x 26 in. [660 x 660 mm]). compressive strength fc′, aspect ratio, and applied load eccen-
Further, due to a limit on maximum GFRP tensile strain (that tricity to evaluate implications on the design of GFRP-RC
is, 0.01) by Code Section 10.3.2.1, only a slight increase in column cross sections. For comparison, steel-RC sections
@seismicisolation
@seismicisolation
Note: 1 kip = 4.44 kN; 1 kip-ft = 1.35 kN-m; 1 in.2 = 645 mm2.
were also designed by changing the parameters stated 22 x 22 in. (560 x 560 mm), which decreased to 18 x 18 in.
previously. The yielding strength of steel used was 60 ksi (460 x 460 mm) and 16 x 16 in. (406 x 406 mm) as concrete
(414 MPa) and modulus of elasticity 29,000 ksi (200 GPa), compressive strength increased to 7500 and 10,000 psi (52
whereas the GFRP reinforcement used was compliant with and 70 MPa) respectively. Similarly, the required reinforce-
ASTM D7957 as referenced by ACI CODE-440.11-22. ment decreased from six No. 9 (M29) bars at fc′ = 5000 psi
To compare results with steel-RC, both sections (steel-RC (35 MPa) to four No. 9 (M29) and four No. 8 (M25) at
and GFRP-RC) were subjected to same demand (that is, fc′ = 7500 and 10,000 psi (52 and 70 MPa), respectively
no magnification factors were applied). Therefore, a cross (note: reinforcement used at all three concrete strengths
section of 20 x 20 in. (508 x 508 mm) was used and varied is 1.2%Ag). It was further noticed that, at higher concrete
as required. The reinforcement ratio was kept as close to strength, GFRP-RC sections performed similar to steel-RC.
minimum required 1% as possible throughout the parametric For example, as shown in Table 7, at fc′ = 5000 psi (35 MPa)
study. the required dimensions for GFRP-RC section are 22 x 22 in.
(560 x 560 mm), whereas those for steel-RC are 20 x 20 in.
Design with different fc′ values (508 x 508 mm). However, when concrete strength increased
Four different values of concrete strength (fc′ = 2500, to 7500 psi (52 MPa) and above, the required dimensions for
5000, 7500, and 10,000 psi [18, 35, 52, 70 MPa]) were used. both RC sections are the same.
The cross sections (steel-RC and GFRP-RC) were subjected In contrast, when concrete strength was decreased to
to ultimate axial compressive load of 789 kip (3510 kN) and 2500 psi (18 MPa), the required cross sections signifi-
ultimate moment of 2367 kip-in. (267 kN-m) (no magnifi- cantly increased to satisfy the demand. Similar effects were
cation factors applied). The reinforcement area was kept as observed in the case of steel-RC; however, unlike steel-RC,
close to 1% of gross concrete area as possible. As expected, GFRP-RC dimensions and reinforcement area increased
RC cross-section dimensions significantly decreased more rapidly. For example, the steel-RC section satisfied
with increasing concrete strength. For example, as shown the demand with cross-sectional dimensions equal to 26 x
in Table 7, at concrete strength of 5000 psi (3 MPa), the 26 in. (660 x 660 mm) with eight No. 9 (M29) bars (ρ =
required GFRP cross section to satisfy the demand (Pu = 0.011Ag), whereas GFRP-RC required 28 x 28 in. (710 x
789 kip [3510 kN] and Mu = 2367 kip-in. [267 kN-m]) is 710 mm) with eight No. 9 (M29) bars (ρ = 0.01Ag). It has
@seismicisolation
@seismicisolation
Note: 1 kip = 4.44 kN; 1 kip-in. = 113 kN-mm; 1 in. = 25.4 mm; 1 in. = 645 mm .
2 2
Note: 1 kip = 4.44 kN; 1 kip-in. = 113 kN-mm; 1 in. = 25.4 mm; 1 in.2 = 645 mm2.
Note: 1 in. = 25.4 mm; 1 kip = 4.44 kN; 1 in.2 = 645 mm2.
(√ )
_ ffr 25.4.9.1 states that the development length in compression
d b _ − 340 shall be conservatively taken the same as that for tension as
fc
_____________
ld = cb ω ( US units) in Code Section 25.4.2.1.
_
13 . 6 + The minimum overlap of tie bar ends shall be greater of
db
(23a) 20db or 6 in. (152 mm), as in 25.7.2.3.1.
( 0 . 083 √fc )
_ ffr _
db − 340 Code Section 10.7.6.2 states that the bottom tie shall be
__________________ located not more than one-half the tie spacing above the top
ld = _ cb ω ( SI units)
13 . 6 + of footing or slab; similarly, the top tie shall be located not
db
more than one-half the tie spacing below the lowest hori-
in which cb/db shall not be taken greater than 3.5, and where zontal reinforcement in the slab, drop panel, or shear cap.
cb is the lesser of: (a) the distance from center of a bar to If beams frame into all sides of column, the top tie shall be
nearest concrete surface, and (b) one-half the center-center located not more than 3 in. (76 mm) below the lowest hori-
spacing of bars being developed; db is the nominal bar zontal reinforcement in the shallowest beam.
@seismicisolation
@seismicisolation
DETAILING OF REINFORCEMENT FOR COLUMN The development length equation for GFRP bars in tension
EXAMPLES (Eq. (25.4.2a)) was conservatively adopted for this case as
The Code permits the use of the same equation for devel- well. Also, the bar location modification factor (ω = 1.5)
opment lengths of GFRP bars in compression and tension was used to calculate the development length. The term cb/
(Eq. (25.4.2)). Therefore, the development length equation db in this column resulted in 1.83, which is well below the
for GFRP bars was conservatively adopted as given in Code permissible limit of 3.5. The development length calculated
Section 25.4.2. Also, there is more than 12 in. (300 mm) for non-sway column as per Code Section 25.4.2 resulted
of fresh concrete to be placed below the longitudinal bar in 75 in. (1900 mm). Because the tensile bar stress reaches
being developed in a sway column; hence, the bar location its full capacity at a point (limited by maximum strain 0.01)
modification factor (ω = 1.5) was also used. The Code spec- which is greater than 0.5ffu, and the ratio of area of reinforce-
ifies a maximum limit for the term cb/db as 3.5, which in ment provided to area of reinforcement required is less than
this column resulted as 1.83, well below the permitted limit. 2.0, therefore, Class B lap splices were selected. The splice
The development length was calculated using three equa- length was calculated as given in Code Section 25.5.2.1,
tions mentioned in Code Section 25.4.2 and Eq. (25.4.2.1(a)) which resulted a value equal to 98 in. (2490 mm).
governed, which resulted in 64 in. (1625 mm). Code Section 10.7.6.2 explains the distribution of ties in
The tensile bar stress at a point reaches its maximum a beam-column joint. The first tie was conservatively placed
value (limited by maximum strain 0.01), which is greater at 3 in. (76 mm), as required by Code Section 10.7.6.2,
than 0.5ffu; also, the ratio of area of reinforcement provided followed by the required spacing. C-shaped ties were used
to the area of reinforcement required in this example is less in the non-sway column, with an overlap as stated in Code
than 2.0 (that is, 1.18); therefore, Class B lap splices were Section 25.7.2.3.1, which resulted in 10 in. (254 mm) when
used. The splice length was calculated as given in Code using No. 4 (M13) ties.
Section 25.5.2.1 (that is, the greater of 1.3ld, 20db, and 12 in.
[300 mm]), which resulted 84 in. (3134 mm). CONCLUSIONS AND RECOMMENDATIONS
The lap splice calculated for the steel-RC column was In this study, a sway column example originally
33 in. (840 mm), which shows that GFRP-RC columns designed with steel reinforcement was taken from the ACI
require very large splice lengths (2.5 times greater than Design Handbook,13 a companion to ACI 318-19,14 and a
steel-RC). Unlike steel, GFRP bars cannot be bent on site non-sway column example from the textbook by Wight and
and together with more GFRP reinforcement required, they Macgregor.15 These two columns were redesigned with glass
increase the complexity in cage preparation. It is the sole fiber-reinforced polymer (GFRP) reinforcement to show the
responsibility of the contractor to splice column reinforce- implications of ACI CODE-440.11-2211 with both low- (Ef =
ment cages. As shown in Fig. 7, the bottom cage bars are 6500 ksi [44,815 MPa]) and high-modulus (Ef = 8700 ksi
shown slightly tilted just after the start of the beam to differ- [60,000 MPa]) GFRP bars. A limited parametric study
entiate from top bars. was carried out to evaluate the effects of changing values
Code Section 10.7.6.2 explains the distribution of ties in of fc′, cross-section aspect ratio, and eccentricity. Based on
a beam-column joint. It states that the bottom tie shall be the outcomes of this study, the following conclusions were
located not more than one-half the tie spacing above the top drawn:
of footing or slab; similarly, the top tie shall be located not • The stiffness values for GFRP reinforcement result in
more than one-half the tie spacing below the lowest hori- higher moment magnification factors for GFRP-rein-
zontal reinforcement in the slab, drop panel, or shear cap. In forced concrete (RC) compared to steel-RC columns.
this example, the first tie was placed at 3 in. (76 mm) from • The advantage of high modulus/strength of new-gen-
the floor top, followed by others at the required spacing. eration GFRP bars can be beneficial to resist condi-
Two pieces of C-shaped stirrups with minimum overlap as tions of large eccentricities. However, due to limits on
the greater of 20db and 6 in. (152 mm) were used as per maximum tensile strain (0.01 in./in. [0.01 mm/mm]) to
Code Section 25.7.2.3.1. In current column design, No. 4 control column curvature, these benefits are not fully
(M13) bars are used; hence, an overlap of 10 in. (254 mm) used.
is provided. • The compressive strength and stiffness of GFRP rein-
The development length for non-sway GFRP-RC forcement is replaced with an equal area of concrete;
column (18 x 18 in. [460 x 460 mm]) was calculated using hence, bigger cross sections are typically required for
Code-specified Eq. 25.4.2.1(a), 25.4.2.1(b), and 25.4.2.1(c). GFRP-RC columns when compared to steel-RC.
@seismicisolation
@seismicisolation
@seismicisolation
@seismicisolation
https://www.concrete.org/educatorsandresearchers/aciinyourclassroom.aspx
ACI STRUCTURAL JOURNAL TECHNICAL PAPER
Title No. 120-S85
Cover concrete plays an important role in the torsional behavior of showed that cracked concrete exhibits a weaker compres-
reinforced concrete members because the resulting shear stresses sive response than what would be observed during a
are concentrated in these areas. Modeling its behavior is difficult cylinder test, a phenomenon known as compression soft-
due to: 1) compression softening; and 2) the possibility of spalling ening.10 Its discovery led to the issue of spalling in torsion
at high loads. Traditional approaches, which only consider one
being revisited by Hsu and Mo11 in 1985. They suggested
effect or the other, are limited in their ability to model the ultimate
that the tendency of the DCFT to overestimate the torsional
strength and torque-twist response of members over a wide range
of cover thicknesses. This paper presents a mechanics-based model strength of a member when the unspalled dimensions were
which can predict when torsional spalling occurs and quantify its used was not due to cover spalling. Instead, they argued that
effect on a member’s strength and stiffness. Its application within it was because the compression-softening effect—which
a nonlinear analysis framework and a design procedure based on would weaken the response of the concrete under torsional
ACI 318-19 is shown. Using the proposed model together with stresses, but had not yet been discovered in 1974—was not
existing compression-softening models results in improved strength considered in the model. They proposed a new model which
predictions of 187 pure torsion tests found in the literature. was similar to the DCFT but, like analytical tools for shear
modeling, neglected spalling and accounted for compression
Keywords: design; reinforced concrete; softening; space truss; spalling;
torsion.
softening. Their model showed good agreement with exper-
iments whose cover thickness were within a narrow range
INTRODUCTION of values.
Cover spalling has commonly been observed in exper- Despite further improvements in understanding how
iments of reinforced concrete beams subjected to pure reinforced concrete resists shear and torsion, there is
torsion, like those shown in Fig. 1, or torsion in combina- still a lack of consensus on the role of spalling in torsion
tion with other actions.1,2 Spalling, which usually affects the behavior and how to address it. Most analytical models
corners of the cross section but can propagate into the side for torsion developed since the 1980s have followed Hsu
cover as well, has important implications for design.3 For and Mo’s “softened approach” by neglecting spalling and
example, design codes such as ACI 318-19,4 CSA A23.3:19,5 considering compression softening,12-23 while most codi-
and AASHTO LRFD6 have special detailing requirements to fied design provisions4-6 adopt a “spalled approach” that
avoid premature failures caused by torsional spalling, and considers spalling but neglects compression softening like
neglect the cover concrete when determining a member’s the original DCFT. There are weaknesses to both methods.
torsional strength. The strength predictions made by “spalled approaches” can
Early analysis tools for torsion, such as the Diag- be very conservative for members that have thick covers or
onal Compression Field Theory7 (DCFT) developed by contain small quantities of reinforcement, and the stiffness is
Mitchell and Collins in 1974, assumed that the entirety of a usually underestimated. Spalled approaches are also incon-
member’s cover would spall to the depth of the hoop rein- sistent with modeling practices for flexural or shear behavior
forcement at failure. This is because the authors found that because: 1) the assumed section geometry neglects the
if the unspalled section geometry and concrete stress-strain cover; and 2) compression softening is neglected when the
response obtained from a cylinder test were used together, concrete is subjected to biaxial or triaxial stress states. On
the torsional strength would be systematically overesti- the other hand, softened approaches can seriously overesti-
mated. To address this issue, they recommended using the mate the torsional strength of members with thick covers.24
fully spalled dimensions because it improved the quality of A weakness shared by both methods is that spalling is treated
the strength predictions and appeared to be consistent with as a binary issue, with the cover being entirely considered
their experimental observations. The DCFT still forms the or entirely neglected, when experimental evidence—such as
basis for the torsion provisions in the previously mentioned the beams shown in Fig. 1—suggests that the reality is some-
design codes, which share its assumption that the cover fully where in between. The shortcomings of these approaches
spalls at failure.8,9
ACI Structural Journal, V. 120, No. 5, September 2023.
In the 1980s, research investigating the shear behavior of MS No. S-2022-311.R1, doi: 10.14359/51738839, received April 21, 2023, and
reviewed under Institute publication policies. Copyright © 2023, American Concrete
reinforced concrete led to further improvements in under- Institute. All rights reserved, including the making of copies unless permission is
standing how reinforced concrete members resist torsion. obtained from the copyright proprietors. Pertinent discussion including author’s
closure, if any, will be published ten months from this journal’s date if the discussion
Experiments performed on panels subjected to pure@seismicisolation
shear is received within four months of the paper’s print publication.
@seismicisolation
RESEARCH SIGNIFICANCE
Traditional methods of accounting for cover spalling in
reinforced concrete members can be inadequate for members
with very thick or very shallow covers. This paper presents
a simple model which can determine when torsional spalling
occurs and quantify its effect on a member’s strength and
stiffness. Implementing this model into existing analysis
tools and design procedures improves their ability to predict
the torsional strength and torque-twist response of members
over a wide range of cover thicknesses. The improved accu-
racy can be particularly beneficial for the evaluation of
torsion-critical members in existing structures.
PROPOSED MODEL
Consider the cracked concrete member subjected to pure
torsion shown in Fig. 2(a). If the torsion is assumed to be
primarily resisted by circulating shear stresses and not
warping torsion, then the member can be represented by a thin
tube which carries a uniform shear flow around its perimeter.
If tensile stresses in the cracked concrete are neglected, the
shear flow will be the result of diagonal compressive stresses
Fig. 1—Pure torsion experiments by Mitchell and Collins
in the concrete, which are equilibrated by tensile stresses in
showing effect of cover thickness on degree of spalling.1
the longitudinal and hoop reinforcement. These compressive
require engineers to use their judgement when addressing stresses, which vary through the thickness of the tube, can be
issues of torsion because a more general method is not yet represented using an equivalent rectangular stress block7,11
available. with an average stress of f2,av = α1βcsfc′ that acts over a width
Accounting for the effects of both compression softening ao, where fc′ is the concrete strength obtained from a cylinder
and spalling on torsion behavior can address the weaknesses test, βcs is a factor accounting for compression softening, and
of these traditional approaches. However, it is shown later in α1 is a stress block factor which is a function of the prin-
this paper that simply considering both at the same time is cipal compressive strain on the surface, ε2s. If the member
too conservative. The main challenge is that spalling, which is unspalled, the line of action of the shear flow will be at
can but does not always influence the torsional response, is a depth of 0.5ao beneath the outside surface of the cross
difficult to account for because the underlying mechanism section. This is shown in Fig. 2(b).
is not fully understood. One attempt at solving this problem Case 1 in Fig. 2(c) shows a corner detail of a member with
was made by Rahal and Collins, who formulated a model
25 a very thin cover as it resists an applied torque. A simple
to predict when spalling due to torsion would occur and representation of the tensile forces in the hoop reinforce-
implemented it into an analytical tool which also considered ment, T, and the compressive forces which make up the
compression softening. Although their model showed good shear flow, C, is also shown. To satisfy equilibrium at the
agreement with a selection of torsion experiments, their corner where these forces change direction, an additional
spalling check was empirically derived using a small set diagonal force must be present. This force will be compres-
of torsion tests and is relatively complex to use for design. sive because the tensile forces in the hoops will be closer to
Furthermore, it is unable to predict how much concrete is the surface than the compressive forces from the shear flow.
lost after spalling took place, which, as seen in Fig. 1, is not In this situation, the cover is not expected to spall because
always the full cover. there are no tensile stresses in the concrete.
This paper attempts to improve design and analysis prac- Case 2 in Fig. 2(c) again shows the corner detail of a
tices for torsion by presenting a general methodology which member subjected to torsion. Here, the cover is thicker than
considers both compression softening and spalling. A central in Case 1 and hence, the compressive forces which make up
part of the methodology is a new model which derives the the shear flow will be closer to the outside surface than the
mechanism of torsional spalling from first principles. In tensile forces in the hoops. The diagonal force in the corner
addition to having the ability to predict when spalling occurs, will now be tensile, and if the resulting tensile stress exceeds
the model also considers and quantifies the resulting loss of the cracking stress of the concrete, diagonal cracking causing
concrete, which is something that has not been investigated spalling will occur. For a 90-degree corner with a uniform
@seismicisolation
@seismicisolation
Ah pcp
_
ρh = (4)
Acp
s
= ft′(1 − 0.8
)
_ f2
fcr (5) Fig. 3—Simplified model for section loss following cover
f2,max
spalling.
Ah fh _____
In Eq. (5), fcr is the cracking strength; ft′ is the uniaxial 2 _ scrs ≥ 4 √ fc′ (1 − 0.8 α1); fc′ in psi (8b)
tensile strength; f2 is the accompanying principal compres-
sive stress; and f2,max is the peak compressive stress after
If Eq. (8) is satisfied and the centroid of the shear flow
considering compression softening, equal to βcsfc′. Note
falls outside of the hoop reinforcement, the corner of the
that unlike Kupfer and Gerstle’s original model, the soft-
cross section will spall off and no longer contribute to the
ened compressive strength is used instead of the uniaxial
member’s torsional resistance. An idealization of the spalled
compressive strength. This is because the concrete at the
member is shown in Fig. 3. The distance between the failure
corner is in a triaxial state of stress instead of being biaxially
planes that define the spalled concrete is the crack spacing,
loaded as in Kupfer et al.’s experiments.28
which matches the assumption used to derive Eq. (8), and
For a member subjected to torsion, f2—and hence, fcr—
within the cover, the failure plane is assumed to extend
will not be constant though the thickness of the cover.
along the centerline of the hoop reinforcement due to micro-
However, the average cracking stress, fcr,av, can be obtained
cracking. These L-shaped failure patterns at the corners of
by using the average compressive stress in the rectangular
the section lead to the following reductions to the overall
stress block, f2,av. Substituting the definitions of f2,av and f2,max
cross-section geometry
into Eq. (5) gives
1 _
psp = pcp− ∑scr,i(1 − 2 √ 2 ) ≥ ph
n _
The uniaxial tensile strength can be taken as (10)
i=1
_____
ft′ = 0.33 √fc′ ; fc′ in MPa (7a)
where Asp and psp are the reduced area and perimeter after
_____ spalling, respectively; and n is the number of corners that
ft′ = 4 √fc′ ; fc′ in psi (7b) have spalled off, which can be conservatively taken as four
for rectangular sections. Asp is generally larger than the area
Spalling will then occur when the applied stress exceeds enclosed by the hoop reinforcement, Aoh, but will approach
the cracking stress, which results in the following condition Aoh for members that have widely spaced cracks relative
when Eq. (2), (6), and (7) are combined to the size of the cross section. Note that Eq. (9) and (10)
assume that the failure plane forms at a 45-degree angle,
Ah fh _____ which is a reasonable assumption for members subjected to
2 _ scrs ≥ 0.33√fc′ (1 − 0.8 α1) ; fc′ in MPa (8a) pure torsion. Instances where the failure plane is unlikely
@seismicisolation
@seismicisolation
flow are distributed over the depth of compression instead depth of compression, td. This can be done by revisiting
of only being concentrated at the centerline, so some tensile some of the assumptions employed by the DCFT in its orig-
stresses in the concrete will always be needed to turn the inal formulation. First, the longitudinal reinforcement is
corner, even if the resultant shear flow does not appear to assumed to be uniformly distributed around the perimeter of
require any. Therefore, all external or salient corners are the member, so it will remain straight as it twists and hence,
@seismicisolation
@seismicisolation
cover thicknesses. Despite the thicker cover, PT6 failed especially well. It should be noted the initial stiffness of both
at the same torque as PT5 and exhibited severe spalling at beams was underestimated because tensile stresses in the
failure. The proposed model correctly predicted that PT5 cracked concrete (that is, tension stiffening) were neglected
did not experience spalling. This is because the shear flow in the analysis.
was contained within the hoop reinforcement for the entire The results of a more detailed study comparing the four
analysis, so diagonal tension did not develop at the corners. modeling approaches against a broader variety of variables
The ultimate strength and ductility were somewhat under- (such as concrete strength, reinforcement strength, aspect
estimated, though this may be because the confinement ratio, and so on) can be found in the first author’s doctorate
provided by the hoop reinforcement was not considered thesis.54
when modeling the concrete’s compressive response. The
ultimate strength and torque-twist response of PT6 were SIMPLIFICATION FOR DESIGN
predicted with excellent accuracy. Like the experiment, the Although the proposed spalling model improves the
cover was predicted to spall, and the calculated post-spalling predictive capabilities of the DCFT, its full formula-
response followed the experimentally observed behavior tion is too complex to use in design. In this section, a
@seismicisolation
@seismicisolation
( √ 2 )
_Asp
_ Tu psp
ao = psp 1 − 1 − 4 (24)
α1fcʹ Asp
_1
Ao = Asp
− 2 pspao (25)
@seismicisolation
@seismicisolation
Recently, hybrid reinforcement by combining steel with fiber- govern the design of FRP-reinforced concrete (RC) members
reinforced polymer (FRP) bars has emerged as a new system in in most instances. ACI 440.1R-15 (ACI Committee 440
reinforced concrete (RC) constructions. This reinforcement system 2015) and CSA S806-12 (2017) permit using glass FRP
can effectively overcome the ductility and serviceability challenges (GFRP), carbon FRP (CFRP), and aramid FRP (AFRP) bars
of FRP-RC structures. A total of 11 full-scale bridge-deck slabs
in concrete constructions.
were constructed and tested. The test parameters were reinforce-
Basalt fibers have been introduced as a promising addition
ment type, ratio, arrangement, and slab thickness. Moreover, a
comparison between the experimental and predicted deflections to the current types of FRPs. Basalt FRP (BFRP) bars have
from design provisions was carried out to verify the efficiency of relatively greater strength and modulus, comparable costs,
the models for hybrid RC sections. Based on test results, hybrid and higher chemical resistance than GFRP bars (Wu et al.
RC slabs exhibited ductility leading to an ample warning before 2015). Thus, using BFRP bars with a relatively high elastic
failure rather than brittle shear failure observed for FRP-RC slabs. modulus, compared to GFRP bars, would significantly
In addition, hybrid RC slabs displayed good stiffness, service- decrease the amount of reinforcement required and reduce
ability, and load-carrying capacity. Furthermore, test results give the crack width (Elgabbas et al. 2016). Moreover, BFRP-RC
an average bond-dependent coefficient, kb, of 1.27, close to the beams exhibited acceptable deformability when investigated
1.2 recommended by ACI CODE-440.11-22. In addition, some in flexure and shear (Duic et al. 2018).
modifications were proposed to shear equations available in
different design codes to be valid for hybrid RC members without
LITERATURE REVIEW
shear reinforcement.
In FRP-RC members, deeper cracks reduce the contri-
Keywords: basalt fiber-reinforced polymer (BFRP) bar; concrete bridge; bution of uncracked concrete to the shear stress due to the
hybrid reinforcement; shear behavior. lower concrete depth in compression. Moreover, in the trans-
verse direction, FRP bars have lower strength and stiffness,
INTRODUCTION which led to wider cracks and lower aggregate interlock
Bridge-deck slabs are the most critical infrastructure and dowel action supplement to the tensile reinforcement
exposed to harsh environments (deicing salts, humidity, compared to that of an equivalent steel area (El-Sayed et al.
freezing-and-thawing cycles, and chlorides) that make these 2006a). Finally, the total shear strength of FRP-RC members
structures very susceptible to corrosion of steel reinforce- is lower than that of steel-RC members. However, traditional
ment. The associated deterioration can accelerate such failure stirrups are not feasible for constructing slab bridges; conse-
or reduce the expected life span of the structure. Within the quently, the mode of failure may be dominated by shear
past two decades, the most effective way to diminish mainte- (Abdul-Salam et al. 2016). Slabs reinforced with GFRP
nance costs and extend the life span of structures has been to or CFRP bars failed in diagonal tension failure, while the
use fiber-reinforced polymer (FRP) composites as an alter- steel-reinforced slabs failed in ductile flexure mode by steel
native to traditional steel reinforcement in structural compo- yielding followed by concrete crushing (El-Salakawy and
nents, especially where steel corrosion is a major concern. Benmokrane 2004). In addition, increasing the reinforce-
In addition to corrosion resistance, FRP composites have ment ratio significantly improved the shear strength and the
many characteristics over steel reinforcement, such as a post-cracking stiffness (El-Sayed et al. 2005; Matta et al.
high strength-to-weight ratio, excellent fatigue resistance, 2013). Thus, the reinforcement type and its axial stiffness
and nonmagnetic and nonconductive nature, which can be can be confirmed to have a pronounced effect on the shear
used in harsh environments for civil structures. Unfortu- strength of the RC sections.
nately, FRP composites have some drawbacks: a low elastic Gradual failure can be attained by using both FRP and steel
modulus compared to steel (Ef/Es = approximately 0.25) and reinforcements. Therefore, steel reinforcement improves the
linear-elastic behavior up to failure without presenting any ductility by the yielding of steel reinforcement and enhances
yielding plateau, which resulted in the brittle collapse of the serviceability by decreasing the deflection and crack width,
member (Goldston et al. 2016). Many investigations have while FRP reinforcement maintains the load-carrying
been performed to study the overall performance of concrete
ACI Structural Journal, V. 120, No. 5, September 2023.
members reinforced with FRP bars. FRP-reinforced slabs MS No. S-2022-357.R1, doi: 10.14359/51738840, received April 27, 2023, and
reviewed under Institute publication policies. Copyright © 2023, American Concrete
had larger deflections and wider crack widths and depths Institute. All rights reserved, including the making of copies unless permission is
compared to steel-reinforced slabs (Michaluk et al. 1998; obtained from the copyright proprietors. Pertinent discussion including author’s
closure, if any, will be published ten months from this journal’s date if the discussion
Ferrier et al. 2015). Therefore, serviceability criteria often is received within four months of the paper’s print publication.
@seismicisolation
@seismicisolation
@seismicisolation
@seismicisolation
Note: 1 mm = 0.0394 in.; 1 MPa = 145 psi; 1 GPa = 0.145 ksi. @seismicisolation
@seismicisolation
Fig. 2—Stress-strain relationships of reinforcements. D 9SCFB-250 1.40 0.46 0.64 0.15 1.60
9SCFB-120 1.14 1.59 0.38 3.98
of the hybrid deck slabs. The last group comprised three
slabs (9SFCB-200, 9SFCB-250, and 9SFCB-120) with three cracking load, as listed in Table 3. As the load increased,
different thicknesses (200, 250, and 120 mm) to study the more flexural cracks began to develop below or between the
effect of different slab thicknesses. Table 2 summarizes the point loads. With further loading, flexural-shear and shear
details of the slabs. Figure 3 shows the cross sections of all cracks appeared and spread in the shear spans, demonstrating
deck slabs. the shear stresses in the shear span. Once the steel yielded,
the number of cracks increased and propagated, accompa-
Test program and instrumentation nied by the widening of the existing cracks.
The slabs were tested under four-point bending loading Slab 8B showed a larger number of shear cracks, and
up to failure. Figure 4 provides the test setup and schematic the crack length propagated quickly due to the low stiff-
diagram of the slabs. A steel spreader beam was used to ness modulus of the BFRP bar. In addition, Slab 8B yielded
transform the two concentrated loads 900 mm apart, yielding higher crack spacing, more severe cracking, and a lower
a shear span (a) of 900 mm on both sides with a clear span number of cracks than Slab 5B3S, thus tending to suggest
length (L) of 2700 mm. Two half-cylinders with 100 mm substituting BFRP bars with steel bars to improve the axial
diameters were used for loading; however, two full cylinders stiffness of the slab. Increasing the reinforcement ratio
were used to support the slab specimens, as shown in Fig. 4. increases the number of major and minor cracks, improves
The load was applied with displacement control at a constant the crack distribution length, and reduces the average minor
rate of 0.6 mm/min by the hydraulic jack of 500 kN capacity and major crack spacing. This result is clarified by better
to measure the applied loads, with a displacement sensor to bond strength as the number of longitudinal bars increases
measure the corresponding deflection. The load gradient was (Nguyen et al. 2020). Moreover, Slab 5B7S displayed low
established at 5 kN up to an applied load of 100 kN to detect average crack spacing at the same load level compared to
crack widths. After the applied load reached 100 kN, the load Slab 10B6S. The same phenomenon was reported by Ge
gradient was increased to 10 kN until the failure of the slab. et al. (2015); the average crack spacing diminishes with the
Three linear variable differential transformers (LVDTs) were reduction of Af/As.
positioned at the bottom midspan and the two loading points
of each slab, and two LVDTs were installed at the supports Modes of failure
to offset their settlements. Five electrical strain gauges were The combined effect of high shear force and bending
attached to the slab surface to measure the concrete strain moment leads to a spatially high-stress area; hence, failure
along the depth of the slab, and four strain gauges were also occurs in the shear span. Four different failure modes were
attached to the surface of the tensile reinforcement. The observed in the experimental tests and are listed in Table 3.
applied loads, deflections, and strain readings at a frequency Mode I: Diagonal tension failure (DTF)—This mode was
of 10 Hz were automatically recorded using a data acquisi- observed only for Slab 8B reinforced with BFRP bars owing
tion system. The crack width was measured using a hand- to the low stiffness of BFRP bars, as shown in Fig. 5. At
held readout microscope with a magnification factor of 40× a high load level, the shear cracks continued to widen due
with an accuracy of 0.01 mm. to the absence of shear reinforcement until failure occurred.
The diagonal shear crack occurred at 100 mm far from the
EXPERIMENTAL RESULTS AND DISCUSSION support, making it approximately 42 degrees with the hori-
Crack patterns and propagation zontal, then extended toward the loading point and widened
Figure 5 shows the crack distribution of the tested deck and propagated, leading to slab collapse. Moreover, the
slabs upon failure. The first crack was initiated in the pure failure was accompanied by local bending of steel bars in
bending moment zone. Generally, the cracking load was the compression zone keeping the slab intact as one part,
recorded at a similar load level for slabs of the same thick- as shown in Fig. 6. This phenomenon led to improved
ness, while increasing the slab thickness increased the ductility and integrity of the slab. In contrast, a previous
@seismicisolation
@seismicisolation
Fig. 6—Bending of steel bars in compression zone: (a) 8B; and (b) 10B6S.
due to flexural cracks, which reduced the moment of inertia. yielding, a secondary stiffness was detected only for hybrid
Slab 8B exhibited a bilinear load-midspan deflection curve deck slabs. This observation proves the concept of the
and degraded faster than hybrid RC slabs due to the lower hybrid section that the significant role of FRP bars brightens
elastic modulus of the BFRP bars, as shown in Fig. 8(a). after steel yielding. For Group B, increasing the reinforce-
Conversely, hybrid RC slabs exhibited trilinear load-deflec- ment ratio increased the post-cracking stiffness, and hence
tion curves owing to the presence of steel bars. As shown decreased the deflection at similar load levels (El-Sayed
in Fig. 8, after steel yielded, a pronounced reduction in the et al. 2006b).
slope of the load-deflection curves as the load increased to Regarding the Af/As, there was a slight influence of the
high levels means that steel reinforcement cannot resist any ratio of Af/As on the stiffness after cracking. However, a
additional load, and only the BFRP reinforcement carried the significant enhancement in stiffness was noticed after steel
load upon failure. In this stage, the deflection of hybrid deck yielding for Slab 5B7S in comparison to Slab 10B6S; thus,
slabs was lower than the deflection of Slab 8B, attributed as Af/As increases, the deflection decreases (Safan 2013). For
to the efficiency of the BFRP bars restricting the excessive example, at the load of 294 kN, the deflection of Slab 10B6S
deflection even after steel yielding. Hence, using hybrid bars was 30% lower than the deflection of Slab 5B7S because the
to reinforce the concrete bridge-deck slabs keeps the load former had a high Af/As, which enhanced the slab rigidity
growing with a reasonable deflection value. by restricting the excessive deflection after the steel yielded.
In Group A, Slab 8B suffered a larger deflection than The foregoing results proved the significant influence of the
Slab 5B3S under the same load level, as shown in Fig. 8(a). Af/As on the post-elastic strength of bridge-deck slabs with
This is attributed to the high elastic stiffness of steel bars sufficient deformability and stiffness.
compared to BFRP bars, which increased the rigidity of the Considering the reinforcement arrangement in Group C,
hybrid slabs. At the yielding load of Slab 5B3S, the deflec- both slabs reinforced in a double layer developed larger
tion decreased by 43% compared to Slab 8B. After @seismicisolation
steel deflection than the slabs arranged in a single layer, as shown
@seismicisolation
Fig. 8—Load-deflection relationships of tested slabs: (a) Groups A and B; (b) Group C; and (c) Group D. (Note: 1 kN =
0.225 kip; 1 mm = 0.0394 in.)
@seismicisolation
@seismicisolation
Vc = 2.5βϕcfcrbwdv
βmod =[ 1 + 1500εx,mod][1000 + sze]
____________
0.4 _ 1300
CSA S6:19
Mu
_
+ Vu
dv
_____________ _ 35sz
εx,mod = ≤ 0.003; s ze = 15 + a ≤ 0.85sz
2(EfAf+ EsAs) g
@seismicisolation
@seismicisolation
This paper presents the characteristics of concrete girders in aggressive service environments. Furthermore, owing to
prestressed with carbon fiber-reinforced polymer (CFRP) tendons restricted budgets and resources, the demand for sustainable
when subjected to compression- and tension-controlled failure materials is commonplace in federal, state, and municipal
modes. To understand the full-range behavior of these girders agencies.
beyond the boundary of design specifications, progressive failure
Carbon fiber-reinforced polymer (CFRP) tendons may
is simulated using an advanced computational approach, agent-
be employed for prestressed concrete application due to
based modeling. Five bulb-tee girder sections are adopted with a
variable amount of CFRP tendons and their flexural responses are several benefits: corrosion resistance, high strength and light
examined until the intended failure modes are accomplished. The weight, nonconducting and nonmagnetic characteristics, low
rate of capacity reductions in the compression-controlled sections is relaxation, and tailorable formation.6 Findings from labora-
governed by the degree of concrete crushing in the upper flange and tory research were integrated into design specifications6,7
the depth of the girder, whereas the rate in the tension-controlled and numerous highway bridges have been erected with
girders is dominated by the sequential rupture of CFRP without CFRP-prestressed concrete girders around the world.8-10 The
demonstrating size dependency. When the girder concrete cracks, failure mode of CFRP-prestressed concrete girders plays a
locally unstable responses are observed in the compression-con- crucial role in the context of safety. Because both tension-
trolled sections, which are not noticed in the other sections. As far and compression-controlled sections are allowed,6,7 the
as deformability is concerned, both girder configurations are satis-
selection of a certain failure mode is largely dependent on
factory. Upon initiation of the progressive failure processes, the
practitioners’ discretion. Technically speaking, if concrete-
level of safety varies differently depending upon the girder type and
the arrangement of the tendons. The tension-controlled sections crushing governs, the girder suffers progressive failure in
require more activation energy, representing a transition rate from conjunction with sequential crushing across its section until
the initial to damaged states, than their compression counterparts. force equilibrium is not achieved; on the other hand, if CFRP
rupture dominates, the girder collapses as soon as tensile
Keywords: carbon fiber-reinforced polymer (CFRP); failure mode; full- strains exceed the ultimate strain of the tendons. The pros
range behavior; prestressed concrete; safety. and cons of those sections are currently inconclusive, and
the research community has been debating the effectiveness
INTRODUCTION of the design approaches over decades.11 From a cursory
Civil infrastructure constitutes the backbone of a nation’s point of view, compression-controlled sections appear to be
economy and is instrumental in the operation of modern safer than tension-controlled sections; contrarily, contem-
society. Aligning with the practical needs for controlled plating the high tensile strength of CFRP, expected safety in
quality, low maintenance, and efficient geometry, prestressed the tension-controlled sections may be compatible with that
concrete bridges are prevalent and account for a significant of the compression-controlled sections. These unapparent
portion of highway systems. For example, the National aspects are the sources of the foregoing arguments on the
Bridge Inventory of the Federal Highway Administration implications of the failure modes.
reports that 66.5% and 57.9% of constructed bridges in Assessments on the flexural response of CFRP-prestressed
Florida and New York were built with prestressed concrete, concrete members, particularly for compression-controlled
respectively.1 Despite such favorable advantages and popu- sections, are circumscribed by articles that stipulate the
larity, durability remains one of the most critical problems maximum usable strains of concrete (εcu = 0.003 and 0.0035
when managing prestressed concrete members, epitomized in ACI 440.4R-046 and SIMTReC,7 respectively), which
by recent statistics stating that an annual budget of $84,500 were empirically assumed for design purposes12 and do not
to $111,800 was spent per bridge for the last 10 years in the mean the actual failure of the members. To expand the scope
United States.2 Because hot-rolled steel outperforms cold- of investigations over the physical collapse of CFRP-pre-
formed steel from a material standpoint,3 prestressing strands stressed concrete members, an in-depth understanding of
comprising multiple wires are more susceptible to failure in full-range behavior is indispensable beyond the prescribed
comparison with reinforcing bars.4 The collapse of a precast
bridge in Lakeview Drive, Washington, PA, was attributed ACI Structural Journal, V. 120, No. 5, September 2023.
MS No. S-2022-361.R1, doi: 10.14359/51738841, received April 25, 2023, and
to a combination of excessive spalling and corrosion damage reviewed under Institute publication policies. Copyright © 2023, American Concrete
in the bottom flange of a box girder. Attention should thus
5 Institute. All rights reserved, including the making of copies unless permission is
obtained from the copyright proprietors. Pertinent discussion including author’s
be paid to the deterioration of prestressed concrete girders closure, if any, will be published ten months from this journal’s date if the discussion
is received within four months of the paper’s print publication.
@seismicisolation
@seismicisolation
Note: H is girder depth; df is effective depth of tendons; Mn is nominal flexural capacity; compression-controlled section failure mode is crushing of concrete; tension-controlled
section failure mode is rupture of CFRP; 1 mm = 0.0394 in.; 1 kN·m = 0.738 kip∙ft.
Fig. 4—Schematic of sequential failure modes: (a) crushing of concrete; and (b) rupture of CFRP.
in a layer of the agents was greater than the threshold limit, on average (Fig. 5(b)). For the replication of possible live
the layer was deemed to fail and the next layer was reset to load intensities during the course of girder failure, four live
be the ultimate compression fiber with a reduced height of load factors (αL) were adopted from published literature,26,27
the girder. This procedure was repeated until the compres- which were then coupled with the performance levels elabo-
sion zone reached the haunch depth of the girder. For the rated in NCHRP 44028: Fully Operational (FO) at αL = 1.75,
tension-controlled sections (Fig. 4(b)), the same procedure Operational (OP) at αL = 1.00, Life Safety (LS) at αL = 0.50,
was applied, except when a tensile strain of CFRP went over and Near Collapse (NC) and αL = 0.25. The physical inter-
the ultimate strain (εf ≥ εfu): the CFRP layer was regarded pretation of these performance indicators is as follows: FO
to fail without concrete crushing (εc < εcu, where εc is the is a full design load, OP is the typical service state, LS is
compressive strain of concrete), which proceeded until all an extreme event, and NC is a critical situation. Figure 5(c)
CFRP tendons ruptured. provides the distributions of the factored load effects (ME)
versus the factored resistance (MR) of the girders, where
Characterization of load effects the dead load factors of αD = 1.25 (structural components)
To evaluate the performance safety of the benchmark and αW = 1.50 (wearing surface) and the strength resistance
girders under the aforementioned failure scenarios, a rela- factor27 of ϕ = 1.0 were associated with the predefined αL
tionship was characterized between variable load effects and factors. The average ratios of the applied loads and resis-
girder resistance. Taken from actual bridge design projects in tance are charted in Fig. 5(d): Case I with the factored dead
Colorado, Illinois, Michigan, Pennsylvania, and Texas,25 the loads is intended to focus on repair design, whereas Case II
unfactored dead and live load moments (Fig. 5(a)) and the with the unfactored dead loads can be used to infer the load
nominal resistance (Fig. 5(b)) of prestressed concrete girders effects (fractions of girder resistance) representing specific
were plotted. The span length of the sampled girders ranged performance levels for safety assessment.
from L = 12 to 42 m (39 to 138 ft). The dead load of the
girders (all structural components plus wearing surface: MDL RESULTS
and MDW, respectively) generated a 327% higher moment The ramifications of compression- and tension-con-
than the live load (MLL) involving impact factors (IM) and trolled sections are delineated for the full-range behavior of
dynamic load allowance (DLA), on average (Fig. 5(a)). CFRP-prestressed concrete girders. Emphasis is placed on
The nominal resistance of the girders (Mn) was greater than failure particulars, flexural responses, safety appraisals, and
the unfactored load effects (MDL + MDW + MLL) by 53.6%, activation energy with regard to the transition of a girder state.
@seismicisolation
@seismicisolation
Fig. 7—Failure of tension-controlled sections: (a) flexural capacity with four CFRP layers; (b) normalized capacity for BT84;
(c) neutral axis depth with four CFRP layers; and (d) CFRP-area reduction.
@seismicisolation
@seismicisolation
Fig. 11—Normalized sectional response of compression- and tension-controlled girders with moment development: (a) neutral
axis; (b) curvature; (c) concrete strain; and (d) CFRP strain.
of the normalized concrete strains was virtually unrelated to Safety assessment
the failure modes (Fig. 11(c)). Likewise, the strains of the Mensuration—The safety of the girders was estimated
most-tensioned CFRP near the bottom of the girders were using Eq. (3), indirectly expressing the probability of
alike from Mnor = 0 to 0.6 (Fig. 11(d)), implying that there failure,31 together with the extent of damage (Eq. (4))
should be no concerns about the premature failure of these
girder configurations under service loading. In excess of Mnor LN(MN /ME )
__________________
_________________
= 0.6, the CFRP strain of the compression-controlled girder
β =
(3)
√ (
COVR )2+ (COVE )2
forged ahead (εf < εfu) relative to that of the tension-controlled
girder suffering the rupture of the tendons (εf = εfu).
Mn0
_ − MnD
Ω = M (4)
n0
@seismicisolation
@seismicisolation
@seismicisolation
@seismicisolation
QUESTIONS?
E-mail any questions to Journals.Manuscripts@concrete.org.
@seismicisolation
@seismicisolation
ACI STRUCTURAL JOURNAL TECHNICAL PAPER
Title No. 120-S88
This paper reports the results of a comprehensive analytical Therefore, employing effective methodologies to predict the
study implemented to develop deflection prediction methodol- deflection of GFRP-RC members with high accuracy is of
ogies for curvilinear reinforced concrete (RC) members with great importance.
glass fiber-reinforced polymer (GFRP) reinforcement, focusing Two common approaches can be employed to calculate the
on precast concrete tunnel lining (PCTL) segments. The first step
immediate deflection of flexural RC elements: 1) using the
involved modifying the procedures for estimating elastic deflection,
general assumptions of elastic deflection calculation along
cracking moment, and cracked moment of inertia, which were then
introduced for use with curvilinear members. In the next step, three with the effective moment of inertia (Ie); and 2) integration
methodologies of effective moment of inertia, integration of curva- of curvature along the length of the member. Ie considers the
ture, and integration of curvature considering tension stiffening were effective transition between the gross moment of inertia (Ig)
developed for curvilinear members. Then, the analytical results in uncracked regions of a member to the cracked moment
were compared to the experimental database, and a novel method of inertia (Icr) in the cracked part considering the effect of
was developed for predicting deflection in curvilinear GFRP-RC tension stiffening. Branson8 originally assumed the rigidi-
members. In the third and final step, a procedure was developed to ties of the cracked and uncracked parts of a RC element as
adapt the presented methodologies for use with a tunnel segment springs in parallel. His assumption can be written in the form
under real load and boundary conditions. The results indicate that of Eq. (1) as a general model to predict Ie in RC members
the proposed method could predict the deflection of curvilinear
GFRP-RC members with high accuracy.
Ie = k1Ig + k2Icr ≤ Ig (1)
Keywords: curvilinear reinforced concrete (RC) members; deflection;
effective moment of inertia; glass fiber-reinforced polymer (GFRP); precast where k1 and k2 are functions of the ratio of cracking to
concrete tunnel lining (PCTL) segments; reinforced concrete (RC). applied moment (Mcr/Ma), which has been empirically
proposed. Equation (1), based on Branson’s recommenda-
INTRODUCTION tions, can predict the deflection of simply supported straight
Curvilinear reinforced concrete (RC) elements are widely rectangular concrete beams reinforced with typical amounts
used in many types of structures, such as tunnels, bridges, of steel reinforcement with reasonable accuracy. Such
water tanks, and culverts.1 Precast concrete tunnel lining models, however, underestimate the deflection of FRP-RC
(PCTL) segments number among the most frequently used elements, as it was correlated for beams with Ig/Icr smaller
curvilinear RC elements. When a tunnel is bored with a than approximately 4.0, while Ig/Icr in FRP-RC members
tunnel boring machine (TBM), such segments are placed generally ranges between 5 and 25.9 Besides, supposing
sequentially as the boring advances.2 Corrosion is one of the parallel springs for the rigidities of uncracked and cracked
major problems associated with RC structures reinforced sections in Eq. (1) is an incorrect assumption because they
with conventional steel reinforcement. Such issues are exac- are series springs.10,11 By neglecting such wrong assump-
erbated in the corrosive environment of tunnels.3,4 Replacing tions, various researchers tried to modify the values of k1 and
steel reinforcement with glass fiber-reinforced polymer k2 based on the experimental results of FRP-RC beams.12-21
(GFRP) reinforcement is recognized as a viable solution for Bischoff9 developed a new form of equation (Eq. (2)) for
dealing with corrosion issues.5 Serviceability often governs Ie based on the true assumption of series springs for the
the design of flexural members reinforced with GFRP bars cracked and uncracked rigidities in a flexural member
either through cracking, deflection, or stress verification.6
In general, curvature, loading distribution, span length,
and boundary conditions do not play a considerable role in
cracking control and stress verification procedures as they
are mainly related to sectional properties. These issues must
be considered in predicting deflection. Due to commercial ACI Structural Journal, V. 120, No. 5, September 2023.
MS No. S-2022-397.R1, doi: 10.14359/51738842, received May 1, 2023, and
GFRP bars having lower moduli of elasticity than steel rein- reviewed under Institute publication policies. Copyright © 2023, American Concrete
forcing bars, deflection in GFRP-RC flexural members at Institute. All rights reserved, including the making of copies unless permission is
obtained from the copyright proprietors. Pertinent discussion including author’s
service load is generally greater than in steel-RC members.7 closure, if any, will be published ten months from this journal’s date if the discussion
is received within four months of the paper’s print publication.
@seismicisolation
@seismicisolation
Curvilinear longitudinal No. 5 15.0 199 55.1 ± 1.25 1115 ± 60 2.0 ± 0.1
GFRP bars No. 6 20.0 284 52.9 ± 0.6 1068 ± 49 2.0 ± 0.1
No. 5 15.0 199 53.5 ± 1.1 1283 ± 42 2.4 ± 0.1
U-shaped GFRP bars*
No. 6 20.0 284 53.2 ± 2.9 1131 ± 35 2.1 ± 0.0
U-shaped and closed
No. 4 3.0 129 55.6 ± 1.6 1248 ± 74 2.2 ± 0.1
GFRP ties*
*
Reported values are based on applying tension to straight bars manufactured with same process as bent bars.
Note: 1 mm = 0.0394 in.; 1 mm2 = 0.00155 in.2; 1 GPa = 145 ksi; 1 MPa = 145 psi.
@seismicisolation
@seismicisolation
To obtain the vertical deflection, dr, each point should be Cracked moment of inertia
multiplied by sin(θ) in the integral. Therefore, the following Calculating the cracked moment of inertia is essential
equation is obtained to calculate the vertical deflection in a in all the deflection-calculation methods presented. In the
polar coordinate system in a curvilinear member following, the procedures to calculate the cracked moment
of inertia for GFRP-reinforced NSC, HSC, and FRC tunnel
Δθ = ∫(∫ − ϕθRdθ + C1)Rsin(θ)dθ + C2 (7) segments were presented. The contribution of reinforcing
top bars is neglected in the presented procedures as a simpli-
where C2 can be obtained using boundary conditions. When fying and conservative assumption.
it is aimed to calculate the deflection at a certain location, GFRP-reinforced NSC and HSC curvilinear members—
deflection can be obtained using a virtual work method with When there is no axial load, the cracked moment of inertia
this equation does not depend on the applied bending moment on a section
in which a linear stress-strain relationship is considered for
mθMθ concrete in compression. When axial load is present, both
Δ = ∑∫_
E I Rdθ (8)
θ θ the axial load and bending moment in the section affect the
cracked moment of inertia. The cracked moment of inertia
where mθ is the moment induced because of a unit dummy in such conditions can be calculated with the following
load applied at the point where deflection is being calcu- equation when the contribution of the top reinforcement is
lated. It should be noted that the deflections resulting from neglected.
shear and axial forces are neglected as they are generally
insignificant when the span length-to-depth ratio is large.30 Icr = (bd3/3)ka3 + nfAf d2(1 – ka)2 (9)
According to the presented methods, an equation for calcu-
lating the elastic deflection of the test specimens was where Eq. (10) can be used to calculate ka
developed, as reported in Appendix A.*
__________________________
√
N2− 2Nωnρf+ ω2nρf(nρf+ 2) − nρfω + N
_____________________________________
Cracking moment ka =
ω (10)
Cracking moment (Mcr) is one of the most influential
parameters in estimating the deflection in an RC member. where ω = Ecεcbd. As can be inferred from Eq. (10), the
Therefore, predicting the exact cracking moment is of great uncracked depth depends on the level of axial load and the
importance in accurately estimating deflection. Equalizing maximum concrete compressive strain. Equation (11) pres-
the maximum tensile stress in the uncracked section______ to the ents the relationship between the bending moment, axial
maximum tensile capacity of concrete (0.62√fc′ according load, ka, and εc in a section.
to ACI 318-1931) yields the theoretical cracking moment
(Mcr,theo). Shrink restraint in an RC member might lead to
pre-existing tensile stresses in the member, which reduce the
_
(
3 − ka _h _
) (
ka d
Ma = nρfωd(1 − ka ) 9 + N 2 − 3 ) (11)
cracking moment.32 Bischoff and Gross22 reported a range of
0.48 to 1.44 with a median of 0.85 for the ratio of theoret- Inserting ka from Eq. (10) into (11) yields an equa-
ical to experimental cracking moment based on the exten- tion with εc as its unknown variable. Due to complexity,
sive data from the literature for FRP-RC flexural members. however, there is no closed-form solution for that equation.
ACI 318-1931 recommends multiplying the theoretical In such situations, the value of εc can be found by trial and
cracking moment by 0.67 in the deflection-calculation proce- adjustment. Subsequently, the values of ka and Icr can be
dure. The ratio of theoretical to experimental cracking moment calculated. The presented equations are based on the linear
in the tested specimens was 0.7 ± 0.02, 0.58 ± 0.04, and 0.61 concrete stress-strain assumption in compression, which is
± 0.02 for NSC, HSC, and FRC specimens, respectively. valid until approximately 0.7fc′.33 According to a prelim-
There is a need for an extensive study to determine the value inary comparison conducted by the authors for the tested
of cracking moment in different concrete types for FRP-RC specimens, neglecting the effect of axial load led to an error
members. Given the lack of such study, the author’s study of approximately 10% in the value of ka for a given value of
recommends taking Mcr equal to 0.7Mcr,theo for NSC GFRP- Ma when the axial load was below ±0.0045fc′Ag. When the
reinforced tunnel segments and 0.6Mcr,theo for HSC and axial load was increased, neglecting such contributions led
FRC GFRP-reinforced tunnel segments according to the to considerable errors.
Curvilinear GFRP-reinforced FRC members—Finding the
cracked moment of inertia in an FRC section considering
*
The Appendix is available at www.concrete.org/publications in PDF format,
appended to the online version of the published paper. It is also available in hard copy the contribution of fibers requires assuming a stress-strain
from ACI headquarters for a fee equal to the cost of reproduction plus handling at the model for FRC in compression and tension. The authentic
time of the request. @seismicisolation
@seismicisolation
(
Ma = Ffd + Fcr
_
3 )
ka d + 2 dcr
+ Fp ( 2 )+ N _
dcr
_ + h
()
h
( )ka d
− Fc _
2 3
Δ = Δg + δΔcr (16)
(13)
where Δg and δΔcr can be calculated with these equations
(refer to Appendix B)
After obtaining εc and ka, the cracked moment of inertia
considering the contribution of fibers can be calculated with _ 1
the following equation
Δg = EI ∫mθMθ Rdθ (17)
c g
Ma ka d
_
Icr
= Eε (14) n η m M
δ Δcr = ∑ _ i ∫_
c c
θ θ Rdθ (18)
i =1 Ec Icr,i
Deflection prediction using effective moment of
inertia These equations are valid for all types of loading and
The effective moment of inertia proposed by Bischoff9 boundary conditions, provided that no settlement or move-
(Eq. (2)) can be replaced with the value of Iθ in the methods ment has occurred in the supports and that the relationship
presented to calculate deflection when the relationship between the applied load and internal forces can be deter-
between the applied loads and internal forces can be spec- mined. Because variations in the axial-load level changes
ified. This method is referred to herein as Ie. This procedure the value of Icr, using a constant value for Icr might be a
does not consider the effect of load and boundary conditions source of errors in the calculation procedure. To account for
or variations in the axial load in the member. It does consider the effect of axial load on deflection calculation, the cracked
the effect of tension stiffening by using the tension-stiff- sections along the member can be divided into a reasonable
ening factor β. This factor theoretically varies between 0 number of parts. The values of ηi and Icr,i for each part can
and 1 depending on the level of bending moment. Bischoff9 be calculated by obtaining δΔcr and summing δΔcr,i along
recommended using Mcr/Ma as the tension-stiffening factor. the cracked section. When the level of axial load and its
As this method supposes a constant effective moment of variation along the member are not significant, however, the
inertia along the member, it is simple to use, but its accu- minimum value of Icr along the member expected in a section
racy depends on the types of loading, boundary conditions, with the greatest bending moment and the lowest axial load
reinforcement ratio, and level of bending moment. The@seismicisolation
best
@seismicisolation
[0 ]
⎢ ⎥
_π _
π
⎡ ⎤
2 −α 2 −α
η
λ = 2(1 − η) ∫ Cθ Rdθ+ (1 − η )∫ Cθ Rdθ (20)
_ _π
2 −α
θg (1 − η) ∫ mθMθ Rdθ+η
_ 2 0
Δ = EI
π
_
−α
(24)
c cr 2
where Cθ = m(π/2)–αMθ. Equation (20) is valid for those types ∫ mθMθ (1 − β(Mcr
/Mθ ))Rdθ
of load and boundary conditions where the distribution of ⎣θg ⎦
bending moment is symmetrical, and the uncracked section
starts at the supports of an angle θg, followed by a cracked Equation (19) can be used to calculate the midspan deflec-
section from θg to midspan. In addition, Icr of the critical tion once the value of λ has been determined with this equa-
section (often located at midspan) was used as the cracked tion (refer to Appendix B)
moment of inertia for all the cracked sections in Eq. (15),
[0 ]
because the level of axial load was not significant. The angle _π
2 −α
_π
2 −α
θg corresponds to the angle from the support to the point at
λ = 2(1 − η) ∫ Cθ Rdθ+ η ∫ ( Cθ − Cθ g) Rdθ (25)
which Mθ = Mcr. For the load and boundary conditions of the θg
θ____________________
g = Evaluation of presented methods with
( )
experimental data
cot (α ) (2 Cb − Ca ) − √(
2
2PR) + 4 Ca Cb − ( Ca )
______________________________________
2
sin−1 − α Table 3 presents the ratio of the theoretical to the exper-
2PR( cot2(α ) + 1) imental midspan deflection of the tested specimens. In
(21)
addition, Fig. 4 and 5 compare the experimental and
where Ca = 2Pfcot(α) + Pl – 4Mcr; and Cb = PRcot(α). It analytical moment-deflection curves of the specimens. The
should be noted that Mcr is dependent on the level of axial moment-deflection curves were drawn up to 50% of the
load. When the variation in axial load is not significant, experimental bending-moment capacity of the specimens.
however, Mcr obtained from the minimum axial load along Table 3 provides a comparison of the midspan deflection
the member can be used as a conservative assumption for at the loads corresponding to the three reference points of
simplicity. For the tested specimens, the integrals in Eq. (20) 2000 µε, 1.1Mcr, and 0.3Mn, when applicable. The deflection
were obtained and are reported in Appendix B. in FRC specimens was obtained according to two scenarios
of considering or neglecting the contribution of the fibers.
Deflection prediction using integration of The average and standard deviation are presented __
sepa-
curvature considering tension stiffening rately for each concrete type. In this study, 3 320√ fc′ + 6900
Due to the effect of tension stiffening, the stiffness in was used as the concrete modulus of elasticity according to
the cracked parts of an element is greater than the cracked ACI 363R-10.37
moment of inertia.22 By supposing Iθ = Icr,θ/[1 – ηθβθ According to Table 3, applying the Ie and Integ. methods
(Mcr,θ/Mθ)] in the calculations related to the cracked parts of overestimated the deflection of NSC tunnel segment speci-
the section, the effect of tension stiffening can be considered. mens by 20% and 50%, respectively, on average, for different
In such situations, deflection can be obtained with reference points. In contrast, considering tension stiffening
based on the Integ. TS method underestimated deflection by
mθMθ 33%, on average. In the HSC specimens, the overestima-
Δ = ∑ ∫_ EI Rdθ
Uncrackd c g tion yielded by the Ie and Integ. methods was 3% and 33%,
mθMθ (1 − ηθβθ( Mcr,θ
/ Mθ ) ) respectively. In contrast, considering the tension stiffening
+ ∑ ∫_______________________
Rdθ (22) with the Integ. TS method underestimated the deflection by
Cracked EcIcr,θ
54%. Neglecting the contribution of the fibers in FRC tunnel
This method is referred as Integ. TS herein when β is segment specimens resulted in a significant overestimation
considered as Mcr/Mθ according to the recommendation of of the deflection (42% for Ie and 81% for Integ. methods).
Bischoff and Gross. To simplify Eq. (22), Eq. (16) can be
22 Considering the contribution of the fibers by considering
their effect on the cracked moment of inertia reduced the
@seismicisolation
@seismicisolation
overestimation of the Ie and Integ. methods by 12% and to Table 3. Therefore, although using effective moment
43%, respectively. The Integ. TS method underestimated the of inertia is simpler than the other methods, it yielded the
deflection by 30% and 41%, respectively, on average, when most accurate results with acceptable conservativeness.
neglecting or considering the contribution of the fibers. As mentioned previously, however, the method’s accuracy
Table 3 and Fig. 4 and 5 reveal that the accuracy of these depends on the load and boundary conditions.22
deflection prediction methods depends on the reinforcement
ratio and concrete type. For instance, Ie underestimated the Proposed model
deflection of 7G No. 5 by 5%, while it overestimated the The method based on integration of curvature considering
deflection of 7G No. 6 and 13G No. 5 by 2% and 49%, the tension-stiffening factor equal to Mcr/Mθ seems to be
respectively. Generally, Integ. TS significantly underesti- theoretically correct. As reported in the preceding section, it
mated the midspan deflection, which was more pronounced considerably underestimated the midspan deflection. Under-
at the lower reinforcement ratio. The approach to calculating estimation with such methods is consistent with some studies
deflection yielded relatively more reasonable results for FRC in the literature.22,38 The underestimation was greater at the
specimens (Fig. 5). In general, according to the average for reference points 2000 µε and 1.1Mcr, especially in the NSC
all the specimens and reference points, the ratio of theoret- specimens with lower reinforcement ratios. To help demon-
ical to experimental deflection was 1.11, 1.41, and 0.55 for strate the reason for such underestimation, Fig. 6 presents
the Ie, Integ., and Integ. TS methods, respectively, according the theoretical midspan moment-curvature diagrams of 7G
@seismicisolation
@seismicisolation
of curvature upon cracking was not considered and the βm = m (
fc′ ρf _
______
fc,n′ρfb ( )
) M where m
Mcr
θ
(
fc,n′ρfb)
fc′ ρf
______
≤ 1.0 (26)
@seismicisolation
@seismicisolation
ηn = 1 − n(
_____
()
fc′ρf ) Ig
where n (
fc,n′ ρfb _ Icr
fc′ρf )
fc,n′ ρfb
_____
≥ 1.0 (27) Table 4—Proposed coefficients to be used in
proposed model
where fc,n′ is 80 MPa (11.6 ksi) for HSC and 40 MPa (5.8 ksi) Concrete type m n p q
for NSC and FRC. The modification constants m, n, p, and
NSC 0.14 4.0 1.12 1.0
q were obtained using regression analysis according to
the results for different types of concrete and are reported HSC 0.65 5.0 0.26 0.8
in Table 4. The effect of the reinforcement ratio on the FRC 0.30 7.7 0.70 2.3
tension-stiffening characteristic was considered using the
ratio of ρf to ρfb. Yost et al.20 and Mousavi and Esfahani17 with closed ties and U-shaped ties, identical constants were
used a similar approach to consider the effect of the rein- proposed for both tie configurations. The modified deflection
forcement ratio on the deflection of GFRP-RC beams. In values can be obtained by replacing the values of η and β in
addition, the ratio of fc,n′/fc′ was added to the proposed model the equations presented to calculate the deflection consid-
to consider the effect of variations in concrete strength in ering the tension-stiffening factor. Equation (28) presents the
the specimens. Furthermore, as the difference was minimal modified value of curvature in each section in the proposed
between the stiffness behavior of the specimens reinforced modified model.
@seismicisolation
@seismicisolation
Fig. 7—Comparison of curvature response of different deflection-calculation methods: (a) moment-curvature response; and
(b) curvature along specimen. (Note: 1 mm = 0.0394 in.; 1 kN·m = 0.7376 kip·ft.)
_
(
Mθ
) Mcr
_
ϕθ = EI 1 − ηnβmM
c cr θ
(28) which reduces the value of Mcr/Mθ according to the rein-
forcement ratio, concrete strength, and the modification
factors. Figure 7(b) presents the curvature value along the
Appendix C describes the procedure used to obtain the
tunnel segment specimen obtained using the Ie, Integ., Integ.
modified equations. Figure 7(a) shows the main concepts
TS, and Mod. methods (note that the specimens’ centerline
of the assumed moment-curvature response in the modi-
was selected as the x-axis for a better view). As can be seen,
fied (referred to as Mod. herein) and Integ. TS methods (the
the methods based on the integration of curvature yielded
graphs are for specimen 7G No. 5). Tension stiffening at a
minimal curvature in the uncracked sections. In the cracked
certain bending moment refers to a change in curvature (ΔΦ)
sections, using Icr in the Integ. method eventually increased
relative to the curvature of the cracked member obtained
the curvature right after passing θg; the curvature increased
using Icr. ΔΦmax is the maximum possible tension stiffening
linearly up to the midspan. In contrast, using the effective
at cracking. This is the tension-stiffening value considered
moment of inertia according to Eq. (2) led to a gradual
in Integ. TS right after the formation of the first crack. This
increase in curvature. In the Mod. method, an increase in
can be the main reason for the significant underestimation
curvature after passing from the uncracked region is consid-
of the integration-of-curvature method considering tension
ered, and the curvature is modeled to gradually increase up
stiffening in calculating deflection for the bending moments
to midspan.
near the cracking load. An indirect method was used to
Figure 6 presents a comparison of the moment-
consider the curvature increase when the first crack appeared
curvature obtained using the proposed equations and the
by modifying the value of η by ηn (refer to Appendix C). By
experimental results for 7G No. 5 and 13G No. 6. As can
using ηn, the maximum change in curvature will be limited to
be seen, the modified model fitted well with the exper-
ΔΦm. Increasing the bending moment decreases the tension-
imental moment-curvature of the sections. The moment-
stiffening effect. The ratio of change in curvature at a certain
deflection relationships obtained with the modified model
bending moment (ΔΦ) to the maximum change in curva-
for different specimens appear in Fig. 4 and 5. As shown, the
ture at cracking is known as the tension-stiffening factor.
modified model was quite consistent with the experimental
Because using Mcr/Mθ overestimated the tension-stiffening
results for all the specimens. In addition, Table 3 gives the
effect, the modified tension-stiffening factor βm is proposed,
@seismicisolation
@seismicisolation
Fig. 9—Schematic view of proposed deflection-calculation procedure for GFRP-reinforced PCTL segments under real load and
boundary conditions.
A schematic arbitrary external load is shown, which varies 4. Calculate the cracking moment (Mcr,i) for each point
according to the loading conditions. Note that the displayed according to the axial-load level.
springs and external load do not play a role in the calcula- 5. For the points where M ≥ Mcr,i, calculate the cracked
tion procedure, which will be described in the following and moment of inertia (Icr,i).
presented solely to show a semi-real condition of a segment. 6. For deflection calculation based on the effective moment
The analysis in Step 1 is supposed to yield the forces and of inertia, consider the moment of inertia for each point (Ii)
moments as well as the deformations and rotation of the according to Eq. (32). For the other methods, where M ≥
joints. The following procedure is proposed to obtain the Mcr,i, calculate the moment of inertia based on the deflection-
rotation and deflection diagrams of GFRP-reinforced PCTL calculation method to be used. Equations (33) and (34)
segments. In addition, Fig. 9 shows a hypothetical schematic represent the moment of inertia based on the Integ. and Mod.
view of the proposed procedure. methods, respectively. When M ≤ Mcr,i, consider Ii equal to
1. Divide the tunnel segment into certain number of Ig.
elements (even number) referred to as nseg. In such situa-
tions, Δθ = 2θmax/nseg. Icr,i
______________
Ii = (32)
2. Designate each element as ni, i = 1, 2, 3,…, nseg. The
angle between starting point of the element ni and the starting
point of the segment is iΔθ.
Mcr,i
_
1 − ηiβi M ( )
θ,i
3. Specify the axial load (Ni) and bending moment (Mi) at Ii = Icr (33)
the points corresponding to α + iΔθ, where i = 1, 2, 3,…, nseg
based on the results in Step I (refer to Fig. 5 for the defini- Icr,i
________________
Ii = (34)
( )
tions of α and θmax). Mcr,i
_
1 − ηn,iβm,i M
@seismicisolation θ,i
@seismicisolation
@seismicisolation
@seismicisolation
1 REGISTER
This paper presents analytical investigations into the behavior configurations are frequently adopted for experimental and
of a reinforced concrete column with and without carbon fiber- theoretical investigations.4,5
reinforced polymer (CFRP) confinement when subjected to Numerous buildings built prior to the enforcement of
earthquake and fire loadings. A data set of 100 ground motions seismic codes, particularly before the 1970s, are consid-
covering short and long durations is collected and integrated
ered deficient and may not safely accommodate lateral sway
with 0 to 3 hours of fire exposure. Two strengthening categories
precipitated by earthquakes. Non-seismically designed
are implemented: 1) one to six CFRP layers; and 2) six layers of
CFRP with a 40 mm (1.6 in.) thick insulation. A computational members thus encounter increased vulnerability to failure
platform incorporating autonomous discrete entities is used for and, in lieu of costly reconstruction, rehabilitation is
the simulation of heat transfer, while static pushover and nonlinear preferred to address assorted issues induced by the inade-
dynamic analyses predict the seismic response of the unconfined quate capacity of those nonconforming ones.3 Because the
and confined columns. Thermal gradients are generated across the primary objective of seismic retrofitting is to ameliorate the
column section to identify the physical and mechanical properties strength and ductility of structural elements against exces-
of constituents at elevated temperatures, which are linked with sive drift ratios,6 relevant rehabilitation strategies would
the static and dynamic models. The CFRP-confined column with mitigate the risk of physical impairment. Among prevalent
insulation outperforms its unconfined counterpart from a behav- strengthening techniques that enhance the capacity of rein-
ioral standpoint, specifically for axial capacities, flexural failure,
forced concrete columns, such as enlarging cross sections
energy dissipation, and deformability. The implications of the
and adding steel plates,7 confinement with carbon fiber-
seismic-fire-combined loadings are remarkable in terms of
degrading the load-resisting ability of the columns compared with reinforced polymer (CFRP) sheets is regarded as a convenient
those of the uncoupled actions. The duration of the ground motions and efficient method8 and is widely implemented around the
dominates the development of a relationship between the spectral world.9 Specifically speaking, the efficacy of CFRP strength-
acceleration and drift ratio of the columns. Design recommenda- ening is remarkable for seismic upgrading associated with
tions are rendered to address the limitations of current practice. the ductile failure, inelastic rotations, and energy dissipa-
tion of substandard columns.10 Another notable benefit of
Keywords: carbon fiber-reinforced polymer (CFRP); column; earthquake the CFRP application is that it raises the strength of existing
duration; fire endurance; modeling; seismic performance.
columns without changing stiffness, which is desirable for
preserving the magnitude of seismic forces.11
INTRODUCTION
Seismic events may ravage energy lines and electricity
The seismic failure of structural elements is contingent
networks in built environments and can prompt fires. A clas-
upon the degree of resistance to earthquake intensities and
sical instance is found in the 1906 San Francisco earthquake,
wave directions. Sequential earthquakes (also known as
which was accompanied by multiple fires: the disaster devas-
aftershocks) aggravate the deterioration of load-bearing
tated the city, and 498 people died and more than 225,000
components, thereby reducing the overall performance of
people became homeless.12 Accordingly, extensive research
building structures. Previous studies enunciated that the
has been conducted for the last several decades to figure
amplitude, frequency, and duration of oscillating ground
out the implications of earthquakes, fires, and a combina-
motions are crucial factors controlling the detrimental impact
tion thereof.13,14 Most cases, however, focused on their own
of seismic waves.1,2 Technical regulations allow a certain
individual consequences,8,15 and insufficient efforts were
extent of seismic damage within the boundary of preventing
expended to understand interactions between these hazards,
the collapse of buildings by developing plastic hinges that
which can bring about a significantly high level of damage
alleviate externally generated excitations.3 As far as the
in constructed facilities.16 It is worth noting that the seismic
stability of a structural system is concerned, columns play
resistance of fire-damaged members is not comparable to that
an important role in maintaining force equilibrium. When
of intact members, and the likelihood of collapse increases
seismic loading is applied to a column, the adequacy of axial
capacity needs to be combined with appropriate ductility ACI Structural Journal, V. 120, No. 5, September 2023.
so that the physical failure of the member is retarded until MS No. S-2023-014.R1, doi: 10.14359/51738843, received April 26, 2023, and
reviewed under Institute publication policies. Copyright © 2023, American Concrete
substantial energy is dissipated. The behavior of columns Institute. All rights reserved, including the making of copies unless permission is
near the base where connections are made with floors obtained from the copyright proprietors. Pertinent discussion including author’s
closure, if any, will be published ten months from this journal’s date if the discussion
or footings is of interest and, that being so, cantilevered is received within four months of the paper’s print publication.
@seismicisolation
@seismicisolation
@seismicisolation
@seismicisolation
Title PGA, DS5-75, Title PGA DS5-75 Scale Title PGA, DS5-75, Title PGA, DS5-75, Scale
No. (station) g seconds (station) (g) (seconds) factor No. (station) g seconds (station) g seconds factor
Chalfant
Valparaiso Petrolia Tohoku
1 Valley 0.444 2.20 0.437 27.6 1.31 26 0.498 2.30 1.044 67.5 1.50
(Llol.) (Petrolia) (Koho.)
(CZR)
Coyote
Hokkaido La Habra Tohoku
8 Lake 0.249 1.40 0.132 25.3 1.05 33 0.361 0.60 0.299 85.4 1.01
(Haya.) (Walnut) (Kawa.)
(CLD)
Loma
Maule Whittier Tohoku
10 Prieta 0.542 1.70 0.684 30.2 1.29 35 0.392 1.50 0.900 71.3 1.68
(Angol) (LA116) (Miyak.)
(GSF)
Loma
Maule Landers Tohoku
15 Prieta 0.356 1.60 0.462 51.4 1.34 40 0.171 21.4 0.254 70.4 1.25
(Talca) (DHS) (Iwanu.)
(GGCS)
San
Kocaeli El Mayor Tohoku
17 0.253 9.70 0.244 23.6 0.82 42 Fernando 1.171 5.40 0.199 69.1 0.20
(Ambarli) (Chih.) (Higash.)
(SP)
El
El Mayor Coalinga Tohoku
18 Mayor-Cu- 0.183 15.4 0.244 26.8 1.29 43 0.109 6.10 0.120 81.0 1.42
(Chih.) (2W) (Kamin.)
capah (3)
earthquakes excerpted from Chandramohan et al.27; after- complete set of the manipulated ground motions is visible
ward, their spectrally equivalent motions were identified in Fig. 3.
in tandem with 5%-damped pseudo-acceleration spectra.3
Each of the 50 pairs between the 100 short- and long-dura- HEAT TRANSFER
tion earthquakes revealed minimum squared errors and opti- Conduction and thermal properties
mized scale factors for the individual sets listed in Table 1. The governing equation of heat transfer in the benchmark
By using these spectrally equivalent ground motions, unnec- column may be expressed by Fourier’s law
essary distractors arising from the morphological dissimi-
larity of the paired durations are eliminated (Fig. 2(c)). A ∂ T(x, t)
_ ∂ T2(x, t)
_
@seismicisolation
@seismicisolation = α (2)
∂t ∂ x2
172 ACI Structural Journal/September 2023
Fig. 3—Fifty spectrally equivalent earthquake motion sets.
where T(x,t) is the temperature at location x and time t; others was invariant (Fig. 4(h) and (i)), except for the thermal
and α is the thermal diffusivity (α = λ/(Cpρ), in which λ is decomposition of CFRP at 400°C (752°C). Figures 4(j) to
the thermal conductivity; Cp is the specific heat; and ρ is (l) summarize the thermal diffusivity of the concrete, CFRP,
the density). In compliance with the recommendations of and insulation. Equation (3) reproduces the ASTM E119
preceding research,28,29 the thermal contribution of rein- standard fire,36 as graphed in Fig. 5(a)
forcement in the concrete column was ignored. Figure 4 plots _ _
the thermal properties of the constituent materials predicted T = 750(1 − exp(− 3.79553√th))+ 170.41√th + T0 (3)
by previously reported equations.30-32 While the conduc-
tivity of the concrete and CFRP descended with tempera- where T is the applied temperature in Celsius; th is the
ture (Fig. 4(a) and (b)), that of the insulation was constant heating time in hours; and T0 is the reference temperature
(Fig. 4(c)). The abrupt drop in the CFRP’s conductivity at (T0 = 20°C [68°F]).
400°C (752°F) was attributed to the breaking of chemical
bonds in the resin, accompanied by a chain scission process Formulation and validation
that lowered the residual mass of effective molecules in the A heat transfer model was developed using a computa-
polymeric composition and precipitated interfacial failure tional platform built with discrete entities, which is called
between the fibers and the matrix.33,34 For modeling conve- agent-based modeling. This nontraditional approach is
nience,35 the mechanical resistance of CFRP was disregarded often employed in social science to study a reciprocal
beyond the thermal decomposition temperature of 400°C relationship between autonomous individuals,40 and the
(752°F) that caused the malfunctioning of the composite. concept was useful to predict the interactive behavior of
Unlike the stable case of the concrete (Fig. 4(d)), the specific the column components during a fire. An open-source code,
heat of the CFRP and insulation was temperature-dependent NetLogo, formed the basis of heat transfer in conjunction
(Fig. 4(e) and (f), respectively). The endothermic reactions with the aforementioned material properties. Further expla-
of these materials involving the conversion of molecular nations on the background and implementation of the code
kinetic energy to chemical energy were responsible for such are available elsewhere.41 Figures 5(b) and (c) compare
erratic behavior. 36-38
The variation in the concrete’s density the theoretical temperature of concrete members with and
was marginal up to 1063°C (1945°F) on account of changes without CFRP strengthening against experimental data
in its mineralogical composition and free water contents at obtained from the literature.22,41 The compressive strength
elevated temperatures (Fig. 4(g) ); however, the density of
39
of the concrete spanned between 28 and 39 MPa (4061
@seismicisolation
@seismicisolation
Fig. 5—Validation of heat transfer model: (a) ASTM E119-2019; (b) plain concrete (test: Weerasinghe et al.42); and (c) confined
concrete with insulation (test: Williams et al.22).
and 5656 psi). The thickness of the CFRP and insulation Appendix A* contains detailed information. As shown in
layers was 1 mm (0.0394 in.) and 38 mm (1.5 in.), respec- Fig. 6(a), both the strength and stiffness of the concrete
tively, and their thermal properties were the same as those declined with the increased thermal load. The yield strength
in the present study (Fig. 4). Temperatures were recorded of the reinforcing bars was sustained up to 400°C (752°F),
at 45 mm (1.8 in.) from the concrete surface42 and in the after which noticeable degradation was rendered (Fig. 6(b)).
vicinity of concrete and insulation levels.22 The computed CFRP sheet—For the temperature-dependent mechanical
and measured values agreed across the board. properties of CFRP, Eq. (4) and (5) may be used43
where Ef(T) and ffu(T) are the elastic modulus and tensile
strength of CFRP at temperature T in Celsius. The hyper- 2fc′
______
bolic responses of CFRP dwindled until thermal decomposi- εt′ = E− E
(10)
c 2
tion occurred (Fig. 6(c) and (d)). According to experimental
observations,44 the linearity of stress-strain in CFRP can be where fc and εc are the stress and strain of the confined
preserved in fire
concrete, respectively; and Ec is the elastic modulus of the
core concrete. The compressive strength of the confined
ff(T) = Ef(T)εf (6)
concrete (fcc′) and its maximum strain (εc,max) are shown in
Eq. (11) and (12), respectively10
where ff(T) is the stress of CFRP at temperature T; and εf is
the CFRP strain.
fcc′ = fc′ + ψf3.3κa fl (11)
Confined concrete—Pursuant to ACI 440.2R-17,10 the
confining pressure of core concrete (fl(T)) in the column is
calculated by
εc,max ≤ εccu ≤ 0.01 (12)
fl(T) = 2Ef(T)ntfεfe/D (7)
εccu = εc′(1.50 + 12κb(fl/fc′)(εfe/εc′)0.45) (13)
where n is the number of the confining layers; and εfe is the
effective strain of CFRP (εfe = 0.55εfu). When more layers
where ψf is a reduction factor (ψf = 0.95); κa and κb are the
were applied, the degree of confinement was raised and its
geometry and efficiency factors (κa = κb = 1.0 for circular
sensitivity to elevated temperatures increased (Fig. 6(e)).
columns); and εc′ is the compressive strain of the core
This fact points out that the ramifications of fire are signifi-
concrete (εc′ = 0.002). Substituting Eq. (7) and (11) into
cant, justifying the demand for insulation to retain the effec-
Eq. (8) with the Eurocode model for the unconfined concrete
tiveness of CFRP strengthening until an intended fire rating
yields the temperature-dependent constitutive relationship
is achieved. The stress-strain curve of the confined concrete
of the confined concrete (Fig. 6(f)).
is written as10
Static pushover analysis
{ for εt′ ≤ εc ≤ εc,max}
Ecεc− ( Ec− E2 )2/ (4fc′) for0 ≤ εc ≤ εt′
fc = A plastic hinge model was formulated with the following
fc′+ E2 εc
assumptions3,21,45 to establish a moment-curvature relation-
(8)ship (Fig. 7(a)):
@seismicisolation
@seismicisolation
Hysteretic simulation
Conceptual development—The benchmark column was
modeled as a single-degree-of-freedom (SDOF) system. This
archetypal approximation for dynamic analysis is justified
by the fact that the plastic hinge region (Fig. 7(a)) dominates
the lateral displacement of the column mass.48 A simpli-
fied trilinear load-displacement curve was then constructed
(Fig. 7(b)). In line with ACI 374.2R-13,49 the initial stiffness
of the curve (Ke) was characterized using the point at which
a lateral load equaled 75% of the maximum load (0.75F0),
and the fundamental period of the column was expressed
to be T1 = 2π(P/(gKe))0.5, in which P is the axial load. The
load and displacement of the column at yielding (Fy and uy,
respectively, and Fy = Keuy) were calculated as instructed in
ASCE/SEI 41-17.3 When the maximum displacement (um)
coincided with the value at a post-peak load of 0.75F0, tanta-
mount to a capacity loss of 25%, the column failed.49
Framework—Because the hysteretic approach proposed
by Ibarra et al.50 adequately models in-cycle deterioration
along with seismic energy dissipation that is essential for
elucidating structural damage,7 it was adopted to simu-
Fig. 7—Hysteretic model: (a) idealized behavior of canti- late the behavior of the benchmark column. The extent of
levered column; (b) trilinear envelope curve; and (c) cyclic cyclic degradation can be quantified through a change in the
deterioration. amount of energy dissipation50
• The shear deformation of the benchmark column is
βi = (Ei/(Et− ∑ ij=1Ej))
c
negligible. (18)
• Plane sections remain plane before and after bending.
• Strain compatibility is valid between the constituents of where βi is the in-cycle degradation parameter at the i-th
the column. alternation (one full load reversal comprises the i-th and
• The interfacial slip of CFRP is negligible at the i+1-th alternations [i ≥ 1], as retraced in Fig. 7(c)); Ei and
concrete-surface level. Ej are the dissipated hysteretic energy values (∫F(u)du, 1 ≤
At fire-exposure time th, the lateral capacity and displace- j ≤ i) at the current and previous alternations (that is, ΣEj
ment of the column (F(th) and u(th), respectively) are deter- is the cumulative hysteretic energy from the first to the i-th
mined by alternations); Et is the reference energy (Et = γFyuy, in which
γ is an empirical constant: γ = 120 for UC columns51 and a
F(th) = Mn(th)/h (14) value should be found for confined ones); and c is the rate
of deterioration (c = 1.0 is used for most reinforced concrete
columns51). With the progression of cyclic loadings, the
u(th) = ue(th) + up(th) = (15) degradation parameter (βi) is updated every load reversal
1/3ϕy(th)h2 + (ϕu(th) – ϕy(th))lp(h – 0.5lp) so that the in-cycle deterioration of the column is computed
from an envelope generated from the pushover analysis
where Mn(th) is the ultimate moment; ue(th) and up(th) are the (Fig. 7(c))
elastic and plastic deformations, respectively; ϕy(th) and ϕu(th)
are the yield and maximum curvatures, respectively; and lp is Fi = (1 – βi)Fi–1 (19)
the plastic hinge length, which can be estimated by46,47
where Fi is the load of the post-peak envelope at the i-th
lp,unconf(T) = 0.08h + 0.022db fy(T) (16) alternation commencing from the initial load of the envelope
without cyclic degradation (F0).
Fig. 9—Thermal conduction across benchmark column: (a) temperature field of unconfined column exposed to fire for 3 hours;
and (b) development of temperature at level of concrete surface.
305 mm (12 in.) and h = 2000 mm (6.5 ft), which were Incremental dynamic analysis
confined with CFRP sheets (tf = 1 mm [0.0394 in.], Ef = For the evaluation of seismic performance through the
70.6 GPa [10,240 ksi], and ffu = 849 MPa [123 ksi]). Irre- foregoing time-history model, incremental dynamic analysis
spective of axial load, the static moment envelopes created (IDA) was carried out.55 The intensity of ground motions was
by the pushover model were in agreement with experi- increased in a piecewise manner until the columns collapsed.
mental responses (Fig. 8(a)). For the hysteretic simulation A spectral acceleration at the fundamental frequency of
of the confined column, the energy constant γ was calibrated the columns was labeled as Sa(T1), which was equivalent
against the measured data (Fig. 8(b)). The converged value to the magnitude of an input motion. This numerical tech-
of γ = 240 was then employed to generate full cyclic curves nique repeatedly updated input motions using the recorded
(Fig. 8(c)). Because others reported that the γ constant was data (all short- and long-duration motions in Table 1), and
an invariable property with respect to elevated tempera- ensuing drift ratios were figured out (the execution algo-
tures,17 γ = 240 was taken for all thermally loaded columns. rithm is delineated in Vamvatsikos and Cornell55). When
the maximum drift ratios of the columns reached the preset
NONLINEAR DYNAMIC ANALYSIS limits from the pushover analysis (Fig. 7(b)), incremental
Time-history analysis iterations were terminated and IDA curves (Sa(T1) versus
The SDOF hysteretic model was expanded to conduct a drift ratio) were drawn.
time-history analysis with the elastic stiffness, fundamental
period, damping ratio, and cyclic degradation detailed previ- RESULTS
ously. Nonlinear solutions were sought using the constant Thermal gradient
average acceleration method, also known as the Newmark Figure 9(a) exhibits a temperature field in the UC column
method,53 and the modified Newton-Raphson iteration.54 exposed to the ASTM E119 standard fire for 3 hours. The
A convergence criterion was chosen at a tolerance limit of distribution of temperature was uneven across the section:
10–5 mm (3.9 × 10–7 in.). Although complete methodological the cover concrete outside the closed tie was thermally loaded
procedures are omitted due to the page limit, one can refer to over 1000°C (1832°F), the longitudinal reinforcing bars were
dynamics texts.48 Predictions included time histories for the subjected to a range of 400 to 550°C (752 to 1022°F), and the
displacement, velocity, and acceleration of the unconfined inner side of the column underwent below 400°C (752°F).
and confined columns. The rise in temperature at the level of the concrete surface
is demonstrated in Fig. 9(b). The incipient stage of the UC
@seismicisolation
@seismicisolation
Fig. 13—Median IDA curves as to variable earthquake durations combined with thermal loading (EQ. is earthquake): (a) UC
column; and (b) C6 column.
of disparity was reduced to 11% at 3 hours. This fact corrob- when subjected to the 50 pairs of the spectrally equivalent
orates that the elevated temperatures mitigated the discrep- short and long ground motions under the exposure periods
ancy between the short- and long-duration earthquakes. of 0 and 3 hours. The strength of the columns waned as
Reliance on characteristic duration—Figures 14(a) to (d) the characteristic durations lengthened, which implies that
demonstrate the distribution of cyclically degraded strengths persistent earthquake loadings need to be a consideration for
with DS5-75 for the unconfined and confined columns at failure seismic design. The capacities of the UC and C6 columns,
@seismicisolation
@seismicisolation
Fig. 15—Seismic capacity degradation with earthquake duration in fire (log-scale plots): (a) UC column; and (b) C6 column.
represented by the peak spectral accelerations, are appraised the earthquake duration (Fig. 15(a)). An opposite trend
in Fig. 14(e) to (h). At the threshold duration of 25 seconds, was, however, noticed for the confined column by virtue of
there was a stepwise decrement in the median capacities of the insulation layer (Fig. 15(b)): the less susceptible vari-
the columns (Fig. 14(e) to (h)) and the importance of the ations of the 1- to 3-hour lines in comparison with that of
durations was once again emphasized. the 0-hour line are ascribed to the hysteretic energy dissipa-
tion explained earlier. For a practical appraisal under seis-
DESIGN RECOMMENDATIONS mic-fire-combined loadings at the preliminary design phase,
Shown in Fig. 15 are the sets of strength degradation the quantity of the strength degradation was equated with a
data belonging to the short- and long-duration earthquakes capacity reduction factor (namely, the nominal capacity of a
consolidated with the fire-exposure time (100 ground column is multiplied by the factor to estimate its deteriorated
motions collated in Table 1 were used for each temporal capacity without going through rigorous dynamic inves-
category). Regression lines were added to definitize a rela- tigations). Table 3 arranges those factors in the context of
tionship between the degraded strength and external attri- column types, earthquake durations, and fire ratings, which
butes (dependent and independent variables, respectively). were rounded for the simplification of the exact degrada-
Aligning with Fig. 13(a), the exposure time was more influ- tion amounts. The applicable boundary of the proposal is an
ential in degenerating the strength of the UC column than
@seismicisolation
@seismicisolation
Unconfined without insulation DS5-75 < 25 seconds 0.95 0.85 0.80 0.75
(before strengthening) DS5-75 ≥ 25 seconds 0.90 0.80 0.75 0.70
CFRP-confined with insulation DS5-75 < 25 seconds 0.95 0.95 0.90 0.90
(after strengthening) DS5-75 ≥ 25 seconds 0.90 0.90 0.85 0.85
*
Applicable range: 0 seconds ≤ duration ≤ 85 seconds.
earthquake duration of 0 to 85 seconds in agreement with the established on cyclically accumulated damage plus the
range of the sampled seismic events (Table 1). widened usable drifts that enabled high-level deform-
ability, were a metric to gauge the functionality of the
SUMMARY AND CONCLUSIONS CFRP system.
This paper has discussed analytical procedures to compre- • The intensity of ground motions and earthquake dura-
hend the behavior of reinforced concrete columns with and tions dominated the pattern of conformational loops in
without carbon fiber-reinforced polymer (CFRP) confine- the hysteretic curves. Besides, the seismic-fire-com-
ment in a multi-hazard scenario, comprising earthquake bined loading was a critical factor for dissipating the
and fire loadings. A data set of 100 short- and long-dura- energy of the UC and C6 columns. Upon examining
tion earthquakes was gleaned and partitioned at 25 seconds, the features of the characteristic duration (DS5-75), the
which was amalgamated with a fire-exposure period of 0 to prominence of the persistent earthquake loadings
3 hours. The strengthening scheme involved two categories: was recognized.
1) one to six layers of CFRP sheets; and 2) six CFRP layers • The proposed reduction factors addressed the limita-
with a 40 mm (1.6 in.) thick insulation. Heat transfer was tions of current practice,56 which is reliant on risk-
modeled on a discrete entity platform, and the mechanical targeted maximum considered earthquake responses
response of the columns was predicted by static pushover and (MCER) without allowing for the duration of ground
nonlinear dynamic methods (time-history analysis and incre- motions. The factors can facilitate the design of UC and
mental dynamic analysis [IDA]). After linking the seismic CFRP-confined columns under the combined loadings
approaches with the distribution of temperature from the with an earthquake duration of 0 to 85 seconds and a fire
conduction model, the performance of the unconfined (UC) rating of up to 3 hours.
and confined columns was studied. Technical interests lay
in capacity degradation, load-moment interactions, lateral AUTHOR BIOS
resistance, a relationship between spectral acceleration and Ju-Hyung Kim is a Postdoctoral Fellow in the Department of Civil Engi-
neering at the University of Colorado Denver, Denver, CO. He received his
drift ratio, and hysteretic energy dissipation. Findings were BE, MS, and PhD from the Department of Architecture and Architectural
integrated to suggest performance-based design guidelines. Engineering at Seoul National University, Seoul, South Korea. His research
The following conclusions are drawn: interests include seismic design, performance evaluation, and statistical
analysis of reinforced concrete structures.
• The initial development of temperature was rapid at
the surface of the column concrete (233°C [451°F] Yail J. Kim, FACI, is President of the Bridge Engineering Institute, An
at 1 minute) and, as uneven thermal distributions International Technical Society, and a Professor in the Department of
Civil Engineering at the University of Colorado Denver. He is Chair of
progressed across the section, the growth rate became ACI Subcommittee 440-I, FRP-Prestressed Concrete, past Chair of ACI
stable up to 1063°C (1945°F) at 180 minutes. Whereas Committee 345, Bridge Construction and Preservation, and a member of
the degree of heat transfer was delayed in part by ACI Committees 342, Evaluation of Concrete Bridges and Bridge Elements;
377, Performance-Based Structural Integrity & Resilience of Concrete
wrapping the column with CFRP, insulation was Structures; and 440, Fiber-Reinforced Polymer Reinforcement; and Joint
required to reach a 3-hour fire rating (351°C [664°F] at ACI-ASCE Committee 343, Concrete Bridge Design. He received the
180 minutes). Chester Paul Siess Award for Excellence in Structural Research in 2019.
His research interests include advanced composite materials for rehabilita-
• The thermally degraded properties of the constituent tion; structural informatics; complex systems; and science-based structural
materials were responsible for lowering the pure axial engineering, including statistical, interfacial, and quantum physics.
capacity of the UC column and accelerating its flexural
ACI member Jun Wang is a Postdoctoral Fellow in the Department of
failure. The influence of the fire exposure was marginal Civil Engineering at the University of Colorado Denver. She received her
on the behavior of the confined column with the insu- BS and MS from Northeast Forestry University, Harbin, Heilongjiang,
lation layer (C6), leading to reasonable retention of the China, and her PhD from the University of Colorado Denver. She is a
member of ACI Committee 345, Bridge Construction and Preservation. Her
capacity with a 6.0% loss at the design service load of research interests include multi-object interaction, concrete structures, and
0.2fc′Ag. advanced modeling.
• Through the pushover model, it was construed that the
seismic-fire-combined distress was more detrimental ACKNOWLEDGMENTS
This research is in part supported by the U.S. Department of Transpor-
than the uncoupled actions. The hysteretic analysis tation through the Mountain-Plains Consortium, for which the authors are
confirmed the increased vulnerability of the UC column grateful. Proprietary information is not included to avoid commercialism.
to stability failure when the exposure period was Technical contents presented herein are based on the opinion of the writers
and do not necessarily represent that of others.
enlarged. The bilinear IDA curves of the C6 column,@seismicisolation
@seismicisolation
@seismicisolation
@seismicisolation
@seismicisolation
@seismicisolation
@seismicisolation
@seismicisolation