Research Article  Open Access
Ductility and Ultimate Capacity of Prestressed Steel Reinforced Concrete Beams
Abstract
Nonlinear numerical analysis of the structural behaviour of prestressed steel reinforced concrete (PSRC) beams was carried out by using finite element analysis software ABAQUS. By comparing the loaddeformation curves, the rationality and reliability of the finite element model have been confirmed; moreover, the changes of the beam stiffness and stress in the forcing process and the ultimate bearing capacity of the beam were analyzed. Based on the model, the effect of prestressed force, and Hsteel to the stiffness, the ultimate bearing capacity and ductility of beam were also analyzed.
1. Introduction
Nowadays, the prestressed concrete structural components with their excellent performance have a wide range of applications in housing construction, bridges, water, and other civil engineering projects [1–5]. Prestressed steel reinforced concrete (PSRC) structure is a new type of structure which improves the mechanical properties with the application of the prestressed technology. PSRC structure has double advantages of steel reinforced concrete (SRC) structure and prestressed concrete structure; the builtin steel can improve the structural load carrying capacity and seismic performance, while the application of prestressed technology can ameliorate the performance of the structure under service load conditions (such as crack control and deformation control), so PSRC structure is easy to achieve large span, heavy load structure [2, 6, 7].
As a result, it was quite necessary to make indepth analysis on the static and dynamic characteristics of PSRC beam. A large number of engineering practices show that the finite element method was a quite effective method in numerical analysis methods [8]. In this paper, analysis of the structural behaviour of PSRC beams was carried out using the ABAQUS.
2. Finite Element Simulation
2.1. Introduction of the Finite Element Model
In order to verify the effectiveness and rationality of finite element analysis, first of all, test components in the papers of Wenzhong et al. [9] were used as research objects, and a twospan PSRC beam finite element model was established. The beam crosssection is shown in Figure 1, and the designed strength grade of concrete was C40; the prestressing steel was 15 steel strand, the tensile strength of which was N/mm^{2}, the inner diameter of corrugated duct was 50 mm, and initial tensile stress ; bonded tendons were used; the material grade of Hsteel ( 100 mm × 100 mm × 6 mm × 8 mm) was Q235, and Hsteel was symmetrically arranged in beam crosssection; stirrups are 12@100 and material was HPB235; the tendon profile is shown in Figure 2.
2.2. Prestressing Effect Simulation
The equivalent load method and physical reinforcement method are two kinds of prestressing simulation methods; the equivalent load method is that the prestressed is converted to an equivalent load applied to the structure. But this method uses simplified assumptions, which cannot simulate the prestressed spatial effect on the structure and cannot simulate the complex force structure; physical reinforcement method contains entity segmentation method, the node coupling method, and constraint equation method, in which prestressing is simulated by cooling and initial strain method. Cooling method is relatively simple and can simulate stress losses; initial strain method usually ignores the prestress loss; otherwise, each of the real constants is not equal, which is a larger workload.
In this paper, cooling method was used to simulate the prestressing effect [9]. At first, we set an initial temperature value and a temperature drop value for reinforced unit to make reinforced unit generate contraction deformation; this initial strain will let the reinforced steel pretensioned as prestress in the finite element model. Steel temperature drop value was calculated as follows:where is the value of temperature drop; is the elastic modulus of prestressing steel; is the linear expansion coefficient; is the effective prestressing values.
2.3. The Connection between Steel and Concrete
In this paper, the bond relationship between steel and concrete is simulated through embedded region method in ABAQUS. Embedded region method is usually designated to deal with one or a group of cells located in other units arrested problems; the method can handle tendons and steel mesh [10]. The interaction between cracked concrete and reinforcement that is known as “bondslip” is indirectly included by modifying the postpeak behaviour of concrete which is modeled by “tension stiffening” effect for concrete in tension [11].
2.4. Material Properties
In this article, the concrete constitutive relationship is referenced from Chinese Code for Design of Concrete Structure (GB500102010). The stressstrain curve is shown in Figure 3. The uniaxial tensile stressstrain curve equation of concrete iswhere is the parameter value of curve falling section, obtained according to the specification of different concrete strength; is the uniaxial tensile strength of concrete obtained by material test; is the peak tensile strain corresponding to , the uniaxial tensile strength of concrete; is the uniaxial tensile damage evolution coefficient of concrete.
Uniaxial compressive stressstrain curve equation of concrete iswhere is the parameter value of curve falling section, obtained according to the specification of different concrete strength; is the uniaxial compressive strength of concrete obtained by material test; is the peak tensile strain corresponding to , the uniaxial compressive strength of concrete; is the uniaxial compression damage evolution coefficient of concrete; is the elastic modulus of the concrete.
In the finite element model, the constitutive relationship of Hsteel, ordinary steel, and prestressing tendon is perfectly plastic model, which is shown in Figure 4.
2.5. Element Types and Boundary Conditions
ABAQUS offers a variety of elements; the literature [12], using the “EulerBernoulli” beam elements and “Timoshenko” beam elements to simulate concrete beam and using keyword Rebar to simulate the nonprestressed steel, but using beam element to simulate the prestressed concrete beam, is an approximate simulation method; there are still shortcomings. In this paper, solid element C3D8R was used to stimulate the concrete; Truss element T3D2 is the linear component which can only withstand tensile and compressive loads in space, which is used to simulate nonprestressed steel and prestressed steel; Hsteel girder was modeled by using fournode reducedintegration shell element S4R, which allows shear deformation change with the thickness variating in the thickness direction, its solving method will automatically obey the thick shell theory or shell theory. S4R is finite strain shell element considering the finite membrane strains and rotation is suitable for large strain analysis [13]. Simple support conditions simulating hinge and roller supports were applied to the ends of the beam. The finite element model is shown in Figures 5 and 6.
There are two load steps in the model; the first step only counts the response of beam under the effective prestressed force and its own weight; the second step is gradually applying displacement on the mid of the beam span. In order to prevent the stress concentration from appearing at the bottom of the beam in the loading process, steel pads were applied on the beam. Before applying the vertical displacement, the translational degrees of freedom on the vertical load direction of all nodes on the surface of loading surface were coupled. Considering geometric nonlinear effects the analysis step used “Static General” algorithm.
3. Comparative Analysis of Finite Element Results
After calculation, loaddeflection curves of the midspan of PSRC beam obtained from FE analysis and experimental results are compared in Figure 7.
Figure 7 shows that the maximum capacity of beam () in finite element model is slightly larger than the experimental values. The more results can be calculated by using the energy equivalent method [13]; the yield load () of experimental beam and finite element model is 234.10 kN and 245.52 kN, respectively, the difference value of which is 4.87%; the yield displacement experimental beam is 6.58 mm, the yield displacement of finite element model is 6.86 mm, the difference value is 4.26%, and both error values are controlled within 5%. The stiffness and bearing capacity differences between the finite element simulations with the experimental results are mainly caused by the model not considering the bond slipping between steel and concrete, general program limitations, and other factors. Furthermore, the loaddeflection curve in Figure 7 can be divided into four sections, AB is before concrete cracked, BC is from the concrete cracked to the tensile nonprestressed steel, Hsteel yield, CD is from the steel yielded to reaching the ultimate capacity, and DE represents the plastic deformation stage.
4. Section Stiffness and Stress Analysis of PSRC Beam
Beam deflection consists of three parts: first, the deflection of the beam generated by its own weight; second, invertedarch deflection due to prestressed force; third, the deflection caused by the load. In order to investigate the change of crosssection neutral axis and the section stiffness components in the loading process, the horizontal stress results of a series of points across the beam crosssection under different loading displacement were extracted. Figure 8 shows the horizontal stress of a series of nodes on the midspan section under different displacements (0.5 mm, 2.0 mm, 4.0 mm, 8.0 mm, 13 mm, 17 mm, 25 mm, and 40 mm), in which represents the distance from the nodes to the beam bottom.
As Figure 8 shows, PSRC beam deformation can be also broadly divided into four stages.
First Stage. At this point load displacement is small, the compressive stress generated by prestressed force is not offset. As a result, the stress of entire crosssection is compressive stress, and affected by the prestressed force the neutral axis position moved down under the previous stage loading.
Second Stage. With the load displacement increasing, the concrete at the bottom of the midspan section begins to bear tensile stress; however, because the load displacement is smaller, the concrete at the section bottom does not reach the ultimate tensile strength, so there is no crack in the entire section, and the section stiffness can be regarded as the stiffness consisting of prestressed concrete and Hsteel.
Third Stage. With the load displacement increasing, the concrete at the section bottom cracked and gradually stopped bearing tension force; the figure shows that when the deflection increases from 4.0 mm to 13 mm, the neutral axis constantly moves; at this stage the beam section stiffness consists of confined concrete core, perimeter concrete (including tendons conversion section), and Hsteel stiffness.
Fourth Stage. The bottom flange of Hsteel will be yielded with the increase of the loading deflection, which results in the decrease of the section stiffness and the gradual rise of neutral axis. As can be seen from the figure, when the deflection increases from 13 mm to 25 mm, the distance between the neutral axis and beam bottom increases from 0.13 m to 0.18 m. At this point, the beam section stiffness consists of confined concrete core, a part of surrounding concrete, and Hsteel stiffness. In addition, with the loading displacement increasing before the concrete on the top of the beam is crushed, the neutral axis position no longer moves up.
5. Comparative Analysis of PSRC Beam Capacity
In this section, in order to investigate the effect of prestress on the yield strength, ultimate strength and ductility of RC (reinforced concrete) beam, SRC (steel reinforced concrete) beam, and the changes of mechanical properties of beam after builtin Hsteel, four same size numerical models were made, which are prestressed steel reinforced concrete (PSRC) beams, steel reinforced concrete (SRC) beams, prestressed reinforced concrete (PRC) beams, and reinforced concrete (RC) beams, and their mechanical properties were analyzed.
It can be seen from Figure 9 that the loadmidspan deflection curves of SRC beam generally have a twofold feature, which is different from the threeline feature of RC beam; it shows that the SRC beam has relatively high capacity. Compared with the PRC beam, the ultimate bearing capacities of the SRC beam and PSRC beam increase with the addition of internal Hsteel, and the SRC beam and PSRC beam still have high bearing capacity in the plastic stage; however, the bearing capacity of PRC beam after the ultimate load decreases sharply. Compared with SRC beam, there is invertedarch deflection before loading in PSRC beam, and the cracking load is improved, so the cracking load of PSRC beam is 68 kN; the cracking load of SRC is only 26 kN; it can be confirmed that the prestressed force can improve load carrying capacity of the beam under service load conditions.
Ductility is an important index to reflect the deformation and energy dissipation of composite beams; in this paper, the displacement ductility coefficient is used to evaluate the ductility of the beam, the calculation formula of which is as follows:
It can been seen from Table 1 that, compared with the PRC beam, the ductility and capacity of the PSRC beam have been significantly increased, the ductility increased by 56.49%, and the ultimate capacity increased by 35.89%; being affected by the builtin Hsteel, the improvement of the ductility was more obvious; moreover, the comparison of PRC beam and RC beam shows that adding prestressed steel increases the cracking load and section stiffness of beam, but ductility deteriorated. It is clear that prestressed steel can increase the bearing capacity of beam under normal use state, but the ductility of the beam may be reduced.

6. Conclusion
(1)Using ABAQUS finite element model can accurately simulate the stress and deformation of experiment components; under the premise of rational boundary conditions and accurate material constitutive relations, the results have higher credibility.(2)Loadmidspan deflection curves of SRC beam generally have a twofold feature, which is different from the threeline feature of RC beam. The bearing capacities of PSRC beam and SRC beam keep at a high level; even at the plastic stage, it can be seen that PSRC beam and SRC beam have higher capacity and ductility.(3)PSRC beam has invertedarch deflection before loading, which was advantageous to the normal use of beam. Compared with SRC beam, PSRC beam cracks latter and has larger section stiffness and good deformation recovery performance.(4)The prestress can improve the crack resistance of the beam but reduce the ductility. Adding Hsteel can improve the ductility and ultimate strength of beam; furthermore, the increase of ductility is more significant.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors gratefully acknowledge the supports from the Program for Natural Science Foundation of Zhejiang Province (LY16E080006) and the Program for Science and Technology Foundation of Guangzhou City (201704020148).
References
 G. Manfredi, G. Fabbrocino, and E. Cosenza, “Modeling of steelconcrete composite beams under negative bending,” Journal of Engineering Mechanics, vol. 125, no. 6, pp. 654–662, 1999. View at: Publisher Site  Google Scholar
 H. Xiao, A. Li, and D. Du, “Experimental study on ultimate flexural capacity of steel encased concrete composite beams,” Journal of Southeast University (English Edition), vol. 21, no. 2, pp. 191–196, 2005. View at: Google Scholar
 C. Fu, Y. Li, X. Sun, and J. Xu, “Experimental study on seismic performance of prestressed and nonprestressed steel reinforced concrete frames,” Journal of Building Structures, vol. 31, no. 8, pp. 15–21, 2010. View at: Google Scholar
 A. Zona and G. Ranzi, “Finite element models for nonlinear analysis of steelconcrete composite beams with partial interaction in combined bending and shear,” Finite Elements in Analysis and Design, vol. 47, no. 2, pp. 98–118, 2011. View at: Publisher Site  Google Scholar
 N. Zhang and C. C. Fu, “Experimental and theoretical studies on composite steelconcrete box beams with external tendons,” Engineering Structures, vol. 31, no. 2, pp. 275–283, 2009. View at: Publisher Site  Google Scholar
 Z. Chen, X. Fu, X. Yang, S. Chen, and L. Gu, “Study on steel reinforced concrete Vierendeel portal structure—portal structure design of Shenzhen University Science and Technology Building,” Journal of Building Structures, vol. 25, no. 2, pp. 64–71, 2004. View at: Google Scholar
 X. Hu and W. Xue, “Progress in research and application of PSRC beams,” Architecture Technology, vol. 37, no. 5, pp. 336–339, 2006. View at: Google Scholar
 P. Fanning, “Nonlinear models of reinforced and posttensioned concrete beams,” Electronic Journal of Structural Engineering, vol. 1, no. 2, pp. 111–119, 2001. View at: Google Scholar
 Z. Wenzhong, J. Wang, B. Han, H. Ye, and Y. Zhang, “Experimental research on mechanical behavior of continuous prestressed composite concrete beams with encased Hsteel,” Journal of Building Structure, vol. 31, no. 7, pp. 23–31, 2010. View at: Google Scholar
 K. Chen, J. Song, and S. Zhang, “Simulation and analysis for externally prestressed concrete bridge based on ANSYS,” Advanced Materials Research, vol. 243, pp. 1737–1742, 2011. View at: Publisher Site  Google Scholar
 K. Baskar, “Finiteelement analysis of steelconcrete composite plate girder,” Journal of Structural Engineering, vol. 128, no. 9, pp. 1158–1168, 2002. View at: Google Scholar
 ABAQUS, Abaqus Analysis User's ManualVersion 6.8, ABAQUS, 2008.
 M. Murthy, D. R. Mahapatra, K. Badarinarayana et al., “A refined higher order finite element for asymmetric composite beams,” Composite Structures, vol. 67, no. 1, pp. 27–35, 2005. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Chengquan Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.