Design Models For Shear Strengthening of Reinforced Concrete Beams With Externally Bonded FRP Composites: A Statistical Vs Reliability Approach
Design Models For Shear Strengthening of Reinforced Concrete Beams With Externally Bonded FRP Composites: A Statistical Vs Reliability Approach
1 INTRODUCTION
Vr = Vc + Vs + Vf (1)
1
FRPRCS-8 University of Patras, Patras, Greece, July 16-18, 2007
where Vc and Vs, may be calculated according to provisions existing on current design codes,
independently of the FRP strengthening system adopted. The methodology to estimate the design
value of the FRP contribution in shear, Vfd, according to each of the aforementioned design proposals
is briefly described in Table 1. Figure 1 summarizes the adopted notation defining the geometric
properties of a generic beam reinforced in shear with externally bonded FRP.
Fig. 3: Adopted notation to define the main geometric properties of an FRP shear reinforcement
⎛ w ⎞
Vfd = 0.9 ⋅ ε fed ⋅ Ef ⋅ ρf ⋅ bw ⋅ d ⋅ ( cot θ + cot β ) ⋅ sin β Vfd = φ ⋅ψ f ⋅ ⎜⎜ 2 ⋅ tf ⋅ f ⋅ ffe ⋅ (sin β + cos β ) ⋅ d f ⎟⎟
⎝ sf ⎠
2 ⋅ tf ⋅ w f 2 ⋅ tf ⋅ sin β
ρf = (strips) ; ρf = (cont.) φ =0.85 ; ψ f =0.95 (O) ; ψ f = 0.85 (U,S)
bw ⋅ sf bw
ffe = E f ⋅ε fe
0.8 ⋅ ε fe
ε fed = ; γ f =1.2 / 1.3 / 1.35
γf i) Full wrapping configuration (O):
⎪ ⎛ ⎞
2
fcm 3
⎪ 0.17 ⋅ ⎜⎜ ⎟⎟ ⋅ ε fu df − Le d − 2 ⋅ Le
⎩⎪ ⎝ Ef ⋅ ρf 1000 ⎠ k2 = (U) ; k 2 = f (S)
df df
Notation:
ε fed - design value of effective FRP strain; φ - shear strength reduction factor [10];
ε fe - mean value of effective FRP strain; ψ f - additional reduction factor for FRP;
ε fu - FRP ultimate tensile strain; kv - bond reduction coefficient;
γ f - partial factor for FRP reinforcement; k1 - modif. factor regarding the concrete strength;
ρf - FRP reinforcement ratio; k 2 - modif. factor regarding the FRP configuration;
Ef - elasticity modulus of FRP reinforcement; Le - effective bond length of FRP reinforcement;
fcm - concrete average compressive strength; fck - concrete characteristic compressive strength;
2
FRPRCS-8 University of Patras, Patras, Greece, July 16-18, 2007
γ Rd - partial factor for the resistance model (1.2); φR - reduction factor due to local stress in corners;
f - average concrete tensile strength; λ - normalized maximum bond length;
ctm
ffed - design value for the FRP effective stress; Df - stress distribution factor;
ffd - design value for the ultimate FRP stress; ffd ,max - maximum design stress in FRP;
ffdd - design value for the FRP debonding stress; ffu - ultimate FRP tensile stress;
Gfk - bonded joint specific fracture energy; hfe - effective height of the bonded reinforcement;
3
FRPRCS-8 University of Patras, Patras, Greece, July 16-18, 2007
χ=1 χ=1
250 250
ACI design model − Vfd (kN)
fib design model − Vfd (kN)
200 200
χ = 1.5 χ = 1.5
150 150
100 100
50 50
0 0
0 50 100 150 200 250 300 0 50 100 150 200 250 300
Vf,exp (kN) Vf,exp (kN)
300 300
χ=1 χ=1
250 250
CNR design model − Vfd (kN)
200 200
χ = 1.5 χ = 1.5
150 150
100 100
50 50
0 0
0 50 100 150 200 250 300 0 50 100 150 200 250 300
Vf,exp (kN) Vf,exp (kN)
Fig. 4 - Vf,exp vs Vfd scatterplots regarding fib, ACI, CNR and CT design models
A large scatter is observed in the experimental vs predicted design values for all of the considered
analytical formulations. Table 2 summarizes the main descriptive statistical measures regarding the χ
factor, namely minimum (MIN) and maximum (MAX) values, the average (AVG) that represents a
global safety factor associated with the design procedure, the standard deviation (STD) and the
coefficient of variation (COV) that are indicators of accuracy. The first quartile (Q1) that cuts off the
lowest 25% of data, the median (MED) corresponding to the 50th percentile and the third quartile (Q3)
that cuts off the highest 25% of data are also included.
4
FRPRCS-8 University of Patras, Patras, Greece, July 16-18, 2007
The obtained results show that the fib design model presents, in average, the lowest safety factor
while the safest predictions are attained with CNR. The largest scattering is attained by the CNR
model (COV=0.73) while the least scattered model is fib (COV=0.55). The CT model globally presents
a good performance with an average value of χ =1.43 and COV = 0.58.
However, from the structural safety point of view, a classification system based only on the main
descriptive statistics measures regarding the behavior of the χ factor may not provide enough
information to assess the reliability of a design proposal, considering that for structural purposes
having χ=0.5 is worst than χ=2.0, which is not taken into account on the statistical analysis.
To overcome this limitation a weighed penalty classification system was applied, based on the
“Demerit Points Classification” (DPC) model proposed by [14], where a penalty (PEN) is assigned to
each range of χ ratios according to Table 3, and the total of penalties determines the performance of
each design model.
From Table 3 it can be noticed that the fib design model presents the weakest performance, with
the highest number of penalty points corresponding to 40% of Predictions Against Safety (PAS, χ<1),
while the best results are attained by the CNR design proposal with the lowest of number of PAS
(20%). The CNR model also provides the highest number of extremely conservative values (32%).
5
FRPRCS-8 University of Patras, Patras, Greece, July 16-18, 2007
300
0.8 6
χ=1
LRM - Linear Regression Model Normal Density
NLR - Non Linear Regression Model Est. Desnsity
250
5
0.6
fib design model − Vfd (kN)
200 4
Relative Frequency
χ = 1.5
χ fib
150 0.4 3
100 2
0.2
1
50
0.0 0
0
0 50 100 150 200 250 300 0 1 2 3 4 5 6 O S U
Vf,exp (kN) χ fib Reinf. Configuration
300
0.8 6
χ=1
LRM Normal Density
NLR Est. Desnsity
250
5
0.6
ACI design model − Vfd (kN)
200 4
Relative Frequency
χ = 1.5
χ ACI
150 0.4 3
100 2
0.2
1
50
0.0 0
0
0 50 100 150 200 250 300 0 1 2 3 4 5 6 O S U
Vf,exp (kN) χ ACI Reinf. Configuration
300
0.8 6
χ=1
LRM Normal Density
NLR Est. Desnsity
250
5
0.6
CNR design model − Vfd (kN)
200 4
Relative Frequency
χ = 1.5
χ CNR
150 0.4 3
100 2
0.2
1
50
0.0 0
0
0 50 100 150 200 250 300 0 1 2 3 4 5 6 O S U
Vf,exp (kN) χ CNR Reinf. Configuration
300
0.8 6
χ=1
LRM Normal Density
NLR Est. Desnsity
250
5
0.6
CT design model − Vfd (kN)
200 4
Relative Frequency
χ = 1.5
χ CT
150 0.4 3
100 2
0.2
1
50
0.0 0
0
0 50 100 150 200 250 300 0 1 2 3 4 5 6 O S U
Vf,exp (kN) χ CT Reinf. Configuration
Fig. 5 – Analysis results with the RDB (from top to bottom: fib, ACI, CNR and CT design models)
6
FRPRCS-8 University of Patras, Patras, Greece, July 16-18, 2007
The values in Table 4 show that, despite the global improvement in the design models
performance with the RDB, the results follow the same trend as for the IDB analysis, thus ratifying the
consistency of the collected data.
Figure 6 plots the safe (PSS) vs unsafe (PAS) predictions diagrams for both the IDB (a) and the
RDB (b). From their analysis, it is mandatory to emphasize that all the studied design models show a
poor performance taking into account the large amount of unsafe predictions for the design value of
the FRP contribution in shear.
0.75 0.75
0.60
0.74 0.68 0.70
0.80 0.78
0.87 0.94
0.5 0.5
0.25 0.25
0.40
0.26 0.32 0.30
0.20 0.22
0.13 0.06
0 0
FIB ACI CNR CT FIB ACI CNR CT
a) b)
3.4 Influence of other parameters not explicitly considered in the analytical formulations
Such poor performance showed by the aforementioned analytical formulations indicates that the
relative influence of the considered parameters is simulated deficiently and the effect of others, not
explicitly taken into account, should not be neglected.
Figure 7 presents some relationships that are supposed to affect the performance of the analytical
models namely, the shear force gain ratio, Vf,exp ⁄ (Vr,exp−Vf,exp), the influence of the longitudinal
reinforcement percentage, ρsl, and the influence of shear steel reinforcement presence are
investigated.
All the studied analytical formulations seem to show an increase of the χ factor with beams where
the global shear force gain is higher. Such trend is observed for both discrete (DISC) and continuous
(CONT) reinforcement arrangements, being more obvious for the fib and CT models while a more
diffuse pattern is observed within the ACI and CNR models.
From the interaction between the χ factor and the longitudinal steel reinforcement it is found that χ
tends increase with the increment of the Esρs ⁄ Efρf ratio suggesting a major interaction between the
FRP and longitudinal reinforcement. On the plotted diagrams, beams with conventional shear
reinforcement are set aside from those without stirrups decoupling the interactions between these two
phenomena. It proved that both kinds of beams (with and without stirrups) follow the same trend
regarding the χ vs Esρs ⁄ Efρf relationship.
The plots of the χ vs Eswρsw ⁄ Efρf relation suggest that beams without stirrups may have higher χ
factors with the studied analytical models than those with higher shear reinforcement ratio.
Nevertheless it is the authors’ belief that the pattern found in the χ vs Eswρsw ⁄ Efρf scatter when Asw>0
may also be influenced by the longitudinal reinforcement ratio. In a future approach these phenomena
should be investigated, making several analysis of χ vs Eswρsw ⁄ Efρf for different clusters of ρsl.
7
FRPRCS-8 University of Patras, Patras, Greece, July 16-18, 2007
6 6 6
CONT Asw=0
DISC Asw>0 LRM
LRM LRM
5 5 5
4 4 4
χ fib
χ fib
χ fib
3 3 3
2 2 2
1 1 1
0 0 0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 0 20 40 60 80 100 120 140 0 2 4 6 8 10 12
Shear Force Gain Esl ρsl / Ef ρf Esw ρsw / Ef ρf
6 6 6
CONT Asw=0
DISC Asw>0 LRM
LRM LRM
5 5 5
4 4 4
χ ACI
χ ACI
χ ACI
3 3 3
2 2 2
1 1 1
0 0 0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 0 20 40 60 80 100 120 140 0 2 4 6 8 10 12
Shear Force Gain Esl ρsl / Ef ρf Esw ρsw / Ef ρf
6 6 6
CONT Asw=0
DISC Asw>0 LRM
LRM LRM
5 5 5
4 4 4
χ CNR
χ CNR
χ CNR
3 3 3
2 2 2
1 1 1
0 0 0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 0 20 40 60 80 100 120 140 0 2 4 6 8 10 12
Shear Force Gain Esl ρsl / Ef ρf Esw ρsw / Ef ρf
6 6 6
CONT Asw=0
DISC Asw>0 LRM
LRM LRM
5 5 5
4 4 4
χ CT
χ CT
χ CT
3 3 3
2 2 2
1 1 1
0 0 0
0.0 0.5 1.0 1.5 2.0 2.5 3.0 0 20 40 60 80 100 120 140 0 2 4 6 8 10 12
Shear Force Gain Esl ρsl / Ef ρf Esw ρsw / Ef ρf
Fig. 7 – Influence of shear gain, longitudinal and shear reinforcement on the models performance
8
FRPRCS-8 University of Patras, Patras, Greece, July 16-18, 2007
4 CONCLUSIONS
Based on the information available in the literature regarding the shear strengthening of RC beams
with externally bonded CFRP, a comprehensive database was assembled containing experimental
results of 212 beams. The results obtained from a statistical analysis carried out on such database
demonstrate that none of the analytical formulations predicts with enough accuracy the contribution of
the EBR CFRP systems for the shear strengthening of RC beams. A large scatter of the χ=Vf,exp/Vfd
was found within all the design models, even when a reduced database (RDB) was used in the
analysis.
Using the RDB the average of the χ factor varies between 1.4 (fib) and 2.9 (CNR) and the
coefficient of variation is comprehended between 43% (fib) and 57% (CNR). From a statistical point of
view the CT model can be pointed as the one with the best performance since it always combines an
appropriate global safety factor (AVG χ = 1.67) with the one of most least scattered behaviors (COV χ
= 47%). Although the large amount of calculated Vfd values that are against safety suggest that all of
the aforementioned models are still not robust enough for generalized practical design purposes.
A reliability analysis and classification based on structural safety was also implemented. Among
the studied formulations, the fib design model presented the most unsafe results from all, while the
safest results were attained with the CNR design code provisions. CNR also provided the largest
amount of extremely conservative predictions especially for the S type strengthening configuration.
The influence of some parameters not explicitly considered on the analytical models was
assessed, proving that the performance of the aforementioned design models is subordinated to the
global attained shear force gain. Furthermore, the influence of conventional steel reinforcement
(longitudinal and transversal) proved to be significant, and none of the studied analytical models
explicitly considers these parameters to determine the FRP contribution to shear. Other parameters
not taken into account in the analytical formulations may also influence the behavior of strengthened
beam. The collected database provides a significant source for data mining techniques in order to
decouple the interactions between all the phenomena involved. Thus, more investigation in this field is
needed in order to improve the existing design models.
ACKNOWLEDGEMENTS
The study reported in this paper forms a part of the research program “CUTINSHEAR -
Performance assessment of an innovative structural FRP strengthening technique using an integrated
system based on optical fiber sensors” supported by FCT, POCTI/ECM/59033/2004. The first author
acknowledges the support provided by the grant in the ambit of this research project.
REFERENCES
[1] Triantafillou, T. “Shear strengthening of reinforced concrete beams using epoxy bonded FRP
composites” ACI Structural Journal 11, 9 (March-April 1998), 107-115.
[2] Khalifa et al “Contribution of externally bonded FRP to shear capacity of RC flexural member”
Journal of Composites for Construction ASCE 2, 4 (1998), 195-202.
[3] Monti, G., and Liotta, M. “FRP-strengthening in shear: tests and design equations”. em 7th
International Symposium on Fiber Reinforced Polymer (FRP) Reinforcement for Concrete
Structures (FRP7RCS), 2005, ACI Symposium Publication 230.
[4] fib “Bulletin 14 - Externally Bonded FRP reinforcement for RC structures” Technical report,
2001, Task Group 9.3 FRP (fibre reinforced polymer) reinforcement for concrete structures.
[5] ACI “440.2R-02: Guide for the design and construction of externally bonded FRP systems for
strengthening of concrete structures”, 2002, Reported by ACI Committee 440.
[6] CNR-DT200 “Guidelines for design, execution and control of strengthening interventions by
means of fibre reinforced composites”, 2004, National Research Council.
[7] CIDAR - “Design guideline for RC structures retrofitted with FRP and metal plates: beams and
slabs” Draft 3 - submitted to Standards Australia, 2006, The University of Adelaide.
[8] Chen, J.F. e Teng, J.G. - “Shear Capacity of FRP Strengthened RC Beams: FRP Rupture”
Journal of Structural Engineering, ASCE 2003: 129(5): 615–625.
9
FRPRCS-8 University of Patras, Patras, Greece, July 16-18, 2007
[9] Chen, J.F. e Teng, J.G. - “Shear Capacity of FRP Strengthened RC Beams: FRP Debonding”
Construction and Building Materials 2003: 17(1): 27– 41.
[10] ACI 318 “318-02/318R-02: Building code requirements for structural concrete and
commentary”, 2002, Reported by ACI Committee 318.
[11] Bousselham, A., and Chaallal, O. “Shear strengthening reinforced concrete beams with fiber
reinforced polymer: assessment of influencing parameters and required research''. ACI
Structural Journal 110, 2 (March-April 2004), 219-227.
[12] Aprile, A., and Benedetti, A. “Coupled flexural-shear design of R/C beams strengthened with
FRP” Composites Part B: Engineering 35, 1 (January 2004), 1-25.
[13] Lima, J.L.T. “Assessment of the effective strain in FRP laminates used in shear strengthening
of reinforced concrete beams”, MSc thesis (in preparation) 2006, Univ. Minho, Portugal.
[14] Collins, M.P. - “Evaluation of shear design procedures for concrete structures”, 2001, A Report
prepared for the CSA technical committee on reinforced concrete design.
10