Abstract
This paper shows a theoretical model for reinforced concrete beams subjected to torsion and strengthened with composite of carbon fibers (CFRP). It is based on researches of Michael Collins and collaborators and approaches the theory of diagonal compressive field. This model was inserted in a mathematic software, due to the need of interactions to obtain the torsion and torsion angles. The validation of theoretical values was realized through the results of tests obtained from the literature and it was understood that the model is valid and conservative for the ultimate torsion in most of the studies carried out. However, the beams that had the geometric ratios width and height b/h < 0.5 and the ratio between the cover and the smaller dimension of the d/cob element ≤6 had not so accurate results.
Introduction
Today the building industry, universities and researchers are being challenged to make new developments by means of improved technologies and methodologies that make maximum use of the available resources and minimize the cost of projects. Carbon fiber reinforced strengthening of concrete has been widely used during the last years in several countries and can be considered a well-established construction procedure. During the last few years some important contributions by researchers have introduced new design models into this reinforcing system.
This is an extremely flexible technique of strengthening structures, enabling adaptation of the material to different needs in accordance with the given architectural, structural and economic requirements. The technique of strengthening reinforced concrete elements with external bonded composite material is both economic and effective. The composite materials technology offers a wide range of different products, which vary in type, diameter and fiber orientation, and in their epoxy properties. The arrangement of the strengthening is another fundamental parameter that influences beam resistance.
The reinforcement concept to increase the performance of reinforced concrete beams submitted to torsion is not very widespread. Some studies on the torsional behavior of concrete beams reinforced with carbon fiber composites (CFRP) have been recently developed, however they remain inconclusive. The main objective of this study is to examine the parameters that influence the ultimate torsion of reinforced concrete beams strengthened with externally bonded CFRP. Drawing from existing experimental results a new approach for the ultimate torsion is investigated herein.
From the 1990s the study of reinforcement with externally bonded fiber composites on reinforced concrete elements has increased. The investigations, in general, were aimed at analyzing the strength increase of the beams bending or shear (Panchacharam and Belarbi, 2002).
In the early 2000s, were started some tests with beams submitted to torsion with Ghobarah et al. (2002), Panchacharam and Belarbi (2002), and Zhang et al. (2001). Throughout the decades of 2000 and 2010 some studies have revealed experimental results and analytical models, in which the behavior of this type of reinforcement was analyzed, allowing a better understanding of the interaction between the three materials: concrete, steel and carbon fiber composites. These studies differ in several assumptions adopted, such as: theoretical models for the torsion of the reinforced concrete beam, different orientations of the fibers (90°, 45°, or 0°), stirrups spacing, CFRP confinement effect and number of layers involving the section.
This paper proposes a theoretical model for the analysis of strengthened reinforced concrete beams with CFRP subjected to torsion, which is compared with experimental results obtained in the literature (Silva, 2017).
Theoretical models
The Federation Internacionale du Beton furnish empiric expression for strain that solve the ultimate torsion of reinforced beams strengthened with CFRP when completely involved or in U format (lateral and bottom faces), and the CFRP maximum strain for only 90° orientated fibers with the longitudinal beam direction (Fib-Bulletin 14, 2001). Stands out the publications (Ameli et al., 2004; Ghobarah et al., 2012; Mohammadizadeh et al., 2009; Mohamed and Benmokrane, 2016; Silva Filho, 2007; Tudu, 2012), for beams tests of quadrangular and rectangular sections, which are the most used in reinforced concrete structures. However, in the sizing can be considered the contribution of the slabs, thus it has beams with T sections. Chalioris (2006b), Deifalla and Ghobarah (2010), Abduljalil and Sarsam (2012) and Sure (2013) presented their studies for this case. Also tests of beams of box section were realized (He et al., 2014; Hii and Al-Mahaidi, 2006; and format L Deifalla et al., 2013).
With the aforementioned publications it is concluded that the 90° directed reinforcement, which involves the entire section along the element, is more effective than in the case of 0° and for reinforcement in U-shape. The fibers at 45° are more efficient in increasing the strength of the beam when compared to those oriented at 90°. The researchers have been increasingly validating theoretical models through computational modeling (Finite Element Method − MEF), as examples are the researches of Ameli et al. (2004), Ganganagoudar et al. (2016), Neale et al. (2005) and Tudu (2012).
Mitchell and Collins (1974) presented the Generalized Diagonal Compression Field Theory (GDCFT) model for torsion, which was developed by the equilibrium conditions of the generalized spatial truss model, by an analogy with Wagner’s model for steel beams of slender web (Onsongo, 1978). This model is again presented in Collins and Mitchell (1980) in order to emphasize the design for reinforced a prestressed concrete beam submitted to shear and the torsion. Vecchio and Collins (1986) added to the GDCFT the concrete tensile strength by studying the local stresses in the crack openings, and the softening of the concrete compressive strength associated with the concrete tensile strain. The theoretical model proposed in this paper does not consider the concrete tensile strength.
The author’s theoretical model includes the tensile strength of the CFRP in the cited models in order to increase the ultimate torsion. The CFRP maximum stress considered in the model is reduced, because the experimental results found in the literature show that the debonding of CFRP occurs before this material reaches the ultimate strength. The formulation of Chen and Teng (2003) is used to obtain the maximum value of the CFRP at the debonding.
Generalized diagonal compression field theory
This theory assumes that after the cracking, the concrete no longer has the capacity to resist to tension. The tensile forces are only resisted by the longitudinal and transverse reinforcements, and the compressive forces are resisted by the diagonal compression field with the strut angle θ, which is assumed equal to the crack angle.
The strength contribution of the CFRP reinforcements is added in the diagonal compression field model. These reinforcements only resist longitudinal and transverse tensions along with the steel reinforcements. Figure 1 shows all parameters of generalized spatial truss model for a reinforced concrete beam.

Generalized Spatial Truss Model: (a) spatial truss model for reinforced concrete, (b) longitudinal and strut forces, (c) corner forces, (d) strut force, and (e) reinforced concrete beam submitted to torsion and strengthened with CFRP.
The equilibrium of the CFRP forces is performed in an analogous way to steel reinforcements (Figure 1(e)), but differing transversely by the spacing between the CFRP stirrups,
where
Considering a unit concrete cracked panel with thickness

Cracked panel with the longitudinal and transversal strains, principal strains and strain Morh’s circle.
The strain compatibility considers an effective bond between steel-concrete and CFRP- concrete, then
Where
The compression strain is considered positive.
Vecchio and Collins (1986) proposed a new constitutive relationship, that is, considered the reduction of the compressive strength by the tensile strain increase that occurs perpendicularly to the compression strain (Figure 3(a)). This model considers the tension softening versus strain. In the study developed by Rahal and Collins (1996) the softening of the compression strain is also considered (Figure 3(b)).

Concrete compression strength versus ratio between tensile and compression strains (Vecchio and Collins, 1986) and concrete constitutive relationship (Rahal and Collins, 1996).
Figure 3 shows the material softening that is not only due to the terms associated to the concrete compression strength. The parameter of tensile strain was included in order to reduce the limits proposed.
The equations for the concrete constitutive relationship are:
where
The elastic-plastic steel constitutive relationship is gives by (Figure 4(a)),
where

Steel (a) and CFRP (b) constitutive relationship.
The CFRP constitutive relationship is linear, and the ultimate strain is linked with ultimate stress
where
The equation (9) is used to both directions of CFRP reinforcements, where the index
An epoxy resin, to guarantee the efficiency of the link, carries out the bond between the CFRP reinforcements and the concrete. Therefore, according to the tests results obtained in the literature, CFRP usually does not reach its maximum tensile strength, since the concrete substrate material detachment occurs. In Chen and Teng (2001) are described some fracture mechanics models complemented in Chen and Teng (2003) in order to determine the axial CFRP effective stress, which considers the debonding, then:
This model is based on the effective bond length, ultimate bond strength, and the stress distribution factor, respectively, given by:
where
Equation (13) is adopted for when the normalized bond length factor is
where
For λ < 1, see Chen and Teng (2003). When
where
This model analyzes the stirrups CFRP reinforcement or lateral bonded strips, so the reinforcements with orientation of 0° are not included. Therefore, in this paper the reduction of the maximum stress by the bond for reinforcement with orientation 0° is not considered in the theoretical model.
The wall thickness of the full section

Stress rectangular block of the concrete strut: (a) strain, (b) stress, and (c) retangular stress block.
Collins and Mitchell (1980) give expressions the following expression adopting the stress rectangular block of the concrete strut, so
where
The wall thickness is a function of the equivalent thickness
The equations related to the perimeter and the area bounded by the center line of the shear flow are, respectively:
and the area of concrete and perimeter are given, respectively, by
where
With the expression for torsion angle gives by Onsongo (1978), and the expression of the equivalent thickness as a function of the torsion angle, follow
where
Collins and Mitchell (1980) report that the wall of the diagonal compression field should be calculated from center of the stirrups, because the covering resists until cracking occurs, then is neglected. However, this paper will use the entire concrete area, because other authors suggest this consideration, beside that the beam is reinforced externally with CFRP, which guarantees this integrity, then
The orthogonal strain to the direction of the strut diagonal compression field s given by Mohr’s circle, then:
The interactions of the model are calculated by flowchart, presented at Figure 6.

Flowchart of interactions.
The elastic limit of the beam is defined by the torsion and torsion angle at cracking. These parameters allow to consider the complete beam behavior on the theoretical model. The torsional cracking equation adapted for the International System of Units is taken from
Rahal and Collins (1996) furnish for the torsional cracking equation the following equation (SI units):
For gives the torsion angle of cracking, it is necessary to estimate the tube wall thickness
where
The latter was estimated according to the first strain of the DCFSS-CFRP model. Before using the flowchart, it is necessary to calculate the parameters
The model verifies the CFRP debonding of concrete, and if this happens the process is completed.
Analysis of theoretical results
The theoretical was implemented in the software Mathcad v.15 and has as answer the graph of torsion versus torsional angle. The beams used in the implementation resemble section shape, shape and type of reinforcement. The experimental data and results analyzed were obtained from the beams tested by Ameli et al. (2004), Chalioris (2006a), Ghobarah et al. (2002), Hii and Al-Mahaidi (2006), Mohammadizadeh et al. (2009) and Silva Filho (2007).
Experimental results from the literature refer to externally CFRP reinforced beams with stirrups, which vary in width, spacing and number of layers. The sections of the beams were approximately twice as high as the width, noting that only the beam tested by Hii and Al-Mahaidi (2006) was width bigger than the height. The beams had stirrups, upper and lower longitudinal steel bars, and in some samples the beams had lateral longitudinal reinforcement.
Silva Filho (2007) shows two sets of three beams tested, in each one the beams contained the same concrete, reinforcement and configuration of reinforcement. The choice of beams to be analyzed was according to the proximity of the average results obtained from the set of the beams.
The all beams differ in dimensions, reinforcements, concrete strength and steel mechanical and geometrical parameters (Tables 1 and 2).
Geometric and mechanics parameters of analyzed beams.
CFRP parameters of analyzed beams.
Reinforced longitudinally with four 15 cm wide strips.
The results presented in Table 3 are experimental and analytical for the studies of ultimate torsion and the torsion angle, related to that torque. The values vary due to the different beams characteristics (Tables 1 and 2). The objective of choosing these beams is to show how the theoretical model (ultimate value) applies and behaves in different situations.
The ratio between the experimental results and theoretical model results.
The proposed model, developed in Mathcad v.15, shows the behavior of the beams along the load history. It is identified the threshold at the time of the cracking that makes evident the limit of the section to consider the rectangular section as homogeneous, and, later, it is generated the stretch for which the theory of the field of diagonal compression for torsion is valid. The shear flow is resisted by the “wall” of the section, by the steel reinforcements and by the CFRP reinforcements.
Table 3 also shows the ratio between the experimental results and theoretical model results for torsional angle, torsion, and the mean, standard deviation and coefficient of variation values of these parameters.
The beams behavior is presented by the graphs

Beam VT1 of Silva Filho (2007).

Beam VTL3 of Silva Filho (2007).

Beam BCS1 of Mohammadizadeh et al. (2009).

Beam ACS1 of Mohammadizadeh et al. (2009).

Beam FS050D2 of Hii and Al-Mahaidi (2006).

Beam CFS of Ameli et al. (2004).

Beam C2 of Ghobarah et al. (2002).

Beam RaS-FS150 of Chalioris (2006a).

Beam RbS-FS200 of Chalioris (2006a).
It is understood that in the comparison between the experimental results and those of the model, for the behavior there was inconsistency for the C2 beams (Ghobarah et al., 2002), ACS1 and BCS1 (Mohammadizadeh et al., 2009) in relation to the other beams under study. So, it was identified that these beams have the geometric ratios width and height b/h < 0.5 and the ratio between the cover and the smaller dimension of the d/cob element ≤6.
The study carried out shows the relevance of the cover of the reinforced concrete beam in the solution of the proposed model. The model developed in this work considers for the calculation of geometric parameters of the section all the covering of the beam.
Therefore, a new analysis was made consisting of two arguments: do not consider the half of the cover or the total cover. They are considered in the calculation for the perimeters and areas. Then:
Substituting these parameters adapted by the corresponding parameters in the equations presented in item 3, we have the adjusted model for the consideration expressed in this new approach.
In order to complement the new model, we used the suggestion of Ameli and Ronagh (2007) for the CFRP strain, considering that the maximum strain is not on the outer face of the section, thus:
Following are the graphs of beam behavior from the experimental results and the theoretical model are show in Figures 16 to 21.

Beam BCS1 of Mohammadizadeh et al. (2009) not considering half of covering.

Beam BCS1 of Mohammadizadeh et al. (2009) not considering the total cover.

Beam ACS1 of Mohammadizadeh et al. (2009) not considering half of covering.

Beam ACS1 of Mohammadizadeh et al. (2009) not considering the total cover.

Beam C2 of Ghobarah et al. (2002) not considering half of covering.

Beam C2 of Ghobarah et al. (2002) not considering the total cover.
The results and the ratio between the experimental result and the theoretical model are shown in Tables 4 and 5:
Experimental result and analytical model not considering half of the beams coverings of BCS1 and ACS1 by Mohammadizadeh et al. (2009) and C2 by Ghobarah et al. (2002).
Experimental result and analytical model not considering total of the beams coverings of BCS1 and ACS1 by Mohammadizadeh et al. (2009) and C2 by Ghobarah et al. (2002).
Based on the beams with b/h ≥ 0.5 and d/cover > 6, it was found that the experimental and theoretical results are consistent, in addition to generating conservative values for the ultimate torsion (Table 6).
Experimental and theoretical model results for beams with b/h ≥ 0.5 e d/cover > 6.
Therefore, about the analysis that modifies the model in consideration of the thickness of the shear flow and disregards the half or fully the cover, it is concluded that the model initially should not consider, in some cases, all the contribution of the cover.
Table 7 shows the ratios of Table 6 and the ratios of Table 4 determining the statistical parameters for this set.
Experimental and model results.
According to Chen and Teng (2003) method, for the bond between CFRP-concrete substrate, is conclusive for beams with high torsion strength, since the CFRP-resisted portion must be of great proportion, so that can reach the tension provided in the model and, therefore, the material debonding.
Concluding remarks
Most tests results behavior is similar and above the curves obtained with the theoretical model.
The beam C2 experimental curve of Ghobarah et al. (2002) is below the curve of the theoretical model. This beam has asymmetry in longitudinal reinforcements and coverings.
The torsional cracking, in general, is lower for the experimental results.
Beams with two layers had similar behavior in both curves.
Beams in which the width and height ratio is less than 0.5 have the experimental results for torsion below the result of the theoretical model.
It is understood that the thickness of the cover influences the efficiency of the theoretical model. A study was carried out in order to adapt this model, where it is considered the contribution of half or total of the cover. It was observed that the curves for both hypotheses lie below the experimental curves. A better accuracy was also found for not accounting for half the cover.
The parameters ratios in general are superior to unity, which is desirable because it is in favor of safety.
For ratios less than unity, the results are generally close to that value.
The results provided by the theoretical model have better accuracy for the torsion and is not as effective for the torsional angle.
The mean of all ratios is greater than unity.
The ultimate torsion has a coefficient of variation below 20%, which is considered one acceptable limit, however the coefficient of variation of the torsional angle has exceeded this limit.
The standard deviation of the ratios varies between 0.20 and 0.30, but the beams analyzed in this paper have different characteristics.
Figures represent behavior during loading, but the key purpose of the model is found the ultimate values, which is explicit in the all table data. The coefficients of variation obtained for the ultimate values are all less than 15% for torsion and 25% for torsion angle, which validates the consistency of the results obtained with the theoretical model, despite the appearance of the figures.
The results of the theoretical model for CFRP-reinforced beams applied to torsion when compared to the experimental results searched in literature, it is concluded that: the efficiency of the theoretical model depends on the geometric ratios width and height of the beam b/h, and of the ratio between the smallest dimension and the covering of the element d/cover.
Footnotes
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
