Abstract

SAP2000 software was used to build the finite element model of a six-storey-three-span reinforced concrete (RC) frame structure. The numerical simulation of the seismic performance of the RC frame structure incorporating different levels of rebar corrosion was conducted using pushover analysis method. The degradation characteristics of the seismic performance of the corroded structure under severe earthquake were also analyzed. The results show that the seismic performance of the RC frame decreased significantly due to corrosion of the longitudinal rebars. And the interstory drift ratios increase dramatically with the increasing of the corrosion rate. At the same time, the formation and development of plastic hinges (beam hinges or column hinges) will accelerate, which leads to a more aggravated deformation of the structure under rare earthquake action, resulting in a negative effect to the seismic bearing capacity of the structure.

1. Introduction

The building seismic fortification requirements are gradually improving around the world, and the buildings in strict accordance with design and construction code also showed good seismic performance during earthquakes. But in recent years, the houses collapses, highway cracks [1], and casualties caused by earthquakes are still high, such as the Wenchuan earthquake, the India-Pakistan earthquake [2], and the Chilean earthquake [3]. The postdisaster survey data show that collapsed or severely damaged houses are mainly multistorey masonry buildings or concrete frame structures with relatively poor seismic performance [4].

Currently, multistory masonry houses have been gradually restricted to be built and used in large and medium-sized cities, while the reinforced concrete frame buildings are still widely used in schools, office buildings, residential buildings, street shops, and other buildings due to their flexible plan layout and strong adaptability. With the increase of service years, under the induction of external corrosion factors, the material of reinforced concrete frame structure will deteriorate [58], resulting in the durability damage (such as surface cracks [9], carbonization, spalling, and steel corrosion [10]). Among which, the corrosion of steel bar is regarded as the prime factor that affecting the durability of concrete structures [1113]. Corrosion leads to the degradation of geometric parameters and mechanical properties of the rebar and, to some extent, weaken the static bearing capacity of a concrete structure and increase its brittleness. Meanwhile, the seismic performance will be inevitably impaired [1417]. Therefore, it is of important theoretical significance and engineering guiding value to study the degradation law of seismic performance of corroded reinforced concrete frame structures. At the same time, it can also provide reference about the seismic performance evaluation and maintenance reinforcement for the old reinforced concrete frame structures in the seismic area.

In this paper, the finite element model of a six-storey-three-span corroded reinforced concrete (RC) frame structure was built using SAP2000 software [18, 19]. The numerical simulation of the seismic performance of the corroded RC frame structure was conducted using pushover analysis method [20] because pushover analysis method is the most widely used and convenient method for the seismic performance analysis of the middle-low level RC frame structure. Pushover method is based on the structural static elastic-plastic analysis theory, by increasing the lateral load on the inertial force center of each floor of the structure to obtain the relationship between internal forces and deformation response in this process. Finally, using the seismic demand spectrum and capacity spectrum (ATC-40 response spectrum or Chinese code response spectrum) to estimate the performance point index of the structure, one can easily evaluate the seismic performance of the structure [21]. The main advantage of the pushover analysis method is that the elastic-plastic response of the structure can be considered and the calculation results are stable. This paper also analyzed the degradation law of the seismic performance indicators of the corroded RC frame structure under severe earthquake. The results were expected to provide a reference for seismic safety analysis, reliability assessment, maintenance reinforcement, and reconstruction of corroded RC frame structures.

2. Mechanical Properties of Corroded Reinforced Concrete Materials

2.1. Degradation of the Mechanical Properties of Concrete Cover

A large number of experiments and theoretical analyses have shown that the concrete cover of a RC frame structure will be cracked and peeled off due to the role of rust expansion force after the steel corrosion. In addition, it weakens the bond strength between the rebar and concrete, hence reducing the service life of the structure [22, 23]. Considering that the core concrete binding force of an ordinary reinforced cement concrete member is mainly from the protective layer, the cracking and peeling of the concrete cover will reduce its restraint on the core concrete and further weaken the load and deformation capacity of the member. Therefore, the degradation of mechanical properties of the protective concrete cover should not be ignored during structural analysis. Due to the strong random characteristics of the deterioration of the cover concrete, it is difficult to use analytic methods to determine the degree and location of deterioration. In order to simplify the calculation, (1) is used here to calculate the strength of the cover concrete [24, 25].where is the peak value of the compressive strength of the concrete cover after the steel bar corrosion; is the correlation coefficient between the rebar surface shape and its diameter, and it is usually suggested to be 0.1; is the generalized cracking strain of concrete; is the original width of the member; is the cross-sectional width of the corrosion member; is the number of longitudinal rebars; is the total width of the corrosion crack; is rust oxidation products and the volume ratio coefficient before rust, and it can be taken value 2.0; is the width of the corrosion crack; is the depth of steel bar corrosion; is the corrosion loss rate of the rebar section; and is the radius of the steel rebar before corrosion, and it can be taken as the weighted average value when the reinforcement diameters are different.

2.2. Degradation of the Corrosion Steel Bar Mechanical Properties

Usually, because the distribution of rust factors (such as cracks, chloride ions, and sulfate ions [26]) has significant random characteristics, the corrosion status of reinforcement in RC members often appears as pitting corrosion. Numerical analysis usually adopts the method of equivalent uniform corrosion to deal with the degradation of mechanical properties of pit corroded steel bars, so as to improve the accuracy of numerical simulation or theoretical analysis [27]. It is assumed that the cross section of the rebar is in a uniform corrosion state, but its mechanical properties are degradated according to the pit corrosion equation. In order to consider the effect of pitting corrosion on the mechanical properties of corroded rebar, the nominal yield stress and elastic modulus of the corroded rebar are calculated by (2) [28], and the ultimate stress and strain [29] can be calculated by (3).where () is the nominal yield strength of the steel bar before (after) corrosion; () is the nominal modulus of elasticity of the steel bar before (after) corrosion; () is the nominal ultimate strength of the steel bar before (after) corrosion; () is the nominal ultimate strain of the steel bar before (after) corrosion; and stands for the mass loss rate of the corroded steel. The relationship between and can be seen in (4) [28].

3. Pushover Analysis Results of the Noncorroded RC Frame Structure

3.1. Engineering Background and Calculation Parameters

For the purposes of comparison, a model of noncorroded RC frame building was presented here. The building model has six stories with typical story height 3.6 m, and the overall plan area is 24 × 12 m. All beams are 300 × 550 mm. The columns in 1–3 storeys are 500 × 500 mm, while the columns in 4–6 storeys are 400 × 400 mm, rectangular. Figure 1 shows the typical structural layout. It can be seen that the building model is a typical biaxially symmetric structure. Therefore, to simplify the following computation, only the plane frame in axis ③-③ (Figure 1(b)) is selected to perform pushover analysis.

According to Chinese Standard GB50011-2010, Code for Seismic Design of Buildings [30], the building is located in class II site, with the design earthquake group of 2, and the design basic acceleration of ground motion 0.20 g. The loads acting on the building are in strict accordance with Chinese Standard GB50009-2012, Load Code for the Design of Building Structures [31], and the representative values of the gravity load are shown in Table 1.

The design strength of concrete was C30 grade, and the longitudinal reinforcement was Grade II rebar (HRB335), while Grade I rebar (HPB300) was used for the stirrups. Figure 2 shows the geometric parameters and the reinforcement layout of the beams and columns. The detailed design material parameters and reinforcement parameters of the beams and columns are shown in Tables 2 and 3, respectively.

3.2. Pushover Analysis Method

Pushover analysis is a static-nonlinear analysis method [20]. During pushover analysis, a structure is subjected to lateral load which continuously increases through elastic and inelastic behaviours until the ultimate condition is reached. The key thing here is that the plastic hinge theory plays an important role in pushover analysis. Plastic hinges are often used to simulate the nonlinear behaviour of the structure; the plastic development degree of the plastic hinges is used to reflect the plastic development degree of the structure. In the SAP2000 program, yielding and postyielding behaviour can be modelled using either default hinges or user-defined hinges. Each hinge represents concentrated postyield behaviour in one or more degrees of freedom. Default hinges include uncoupled moment (M), torsion (T), axial force (P), shear (V) hinges, and coupled P-M2-M3 (or P-M-M) hinge which yields based on the interaction of axial force and bending moments at the hinge location. For the P-M-M hinge, an interaction (yield) surface should be specified in the three-dimensional P-M2-M3 space that represents the first yielding location for different combinations of axial force P, minor moment M2, and major moment M3. In the SAP2000 program, the built-in default hinge properties for steel members are provided based on FEMA-356 criteria [32], and the built-in default hinge properties for concrete members are generally based on ATC-40 criteria [33]. For the user-defined plastic hinges, XTRACT [34] software is generally utilized to determine the moment-rotation curves for beam and P-M-M interaction curves for columns. When analyzing an RC frame structure, plastic hinge regions are usually assigned to the ends of the members. For the beam element, the plastic hinge is usually yielded only by the bending moment (M), and the plastic hinge is generally considered by the axial force and the two-way bending moment correlation (P-M-M) for the column element. When the default hinge is adopted to simulate the nonlinear behavior of the structure, only the influence of rebar corrosion on degradation law of the plastic hinge parameters can be considered. However, from the above analysis, it can be seen that the degradation of mechanical properties of the protective concrete cover should not be ignored during the structural analysis. The user-defined hinge can fully consider the effects of rebar corrosion and the degradation of mechanical properties of the protective concrete cover as well as the plastic hinge parameters degradation law, so we adopt a user-defined plastic hinges in the nonlinear analysis of the structure.

To determine the development of the plastic hinges on the frame structure, a force-displacement (moment-rotation) curve can be defined according to FEMA-356 criteria, which gives the yield value and the plastic deformation following yield. This is done in terms of a curve with values at five points, A-B-C-D-E, as shown in Figure 3.

As is suggested by Figure 3, there are four line segments (i.e., AB, BC, CD, and DE) in the skeleton curve, representing the elastic stage, the strengthening stage, the unloading stage, and the failure stage, respectively. Moreover, FEMA-356 criteria also set three performance points (i.e., IO, LS, and CP) between the characteristic points B and C. IO is the abbreviation for “Immediate Occupancy,” which means that the structure is in the serviceability limit state. LS is the abbreviation for “Life Safety,” which means that the structure is close to the safety limit state. CP is the abbreviation for “Collapse Prevention,” which means that the structure is close to the collapse limit state.

3.3. Pushover Analysis Results
3.3.1. Comparison of the P-M-M Curves between Default Hinges and User-Defined Plastic Hinges

In view of the biaxial symmetry of the cross section of the noncorroded frame column, Figure 4 presents only the correlation curve of the axial force P and the major moment M3 in the x direction of the column section. It can be seen that the P-M-M curves obtained by user-defined plastic hinge are in good agreement with those by the built-in default hinge in SAP2000 software.

Table 4 provides the parameters of the P-M-M curves of the noncorroded column section. It can be seen from Figure 4 and Table 4 that the computed key-point parameters of the P-M-M curves using user-defined hinge are generally consistent with those acquired by the built-in default hinge in SAP2000 software. The maximum discrepancy is within 5%. This indicated that the user-defined hinge has relatively higher precision and can be used in the subsequent analysis.

3.3.2. Comparison of the Pushover Analysis Results between Default Hinges and User-Defined Hinges

Considering that the load pattern may affect the pushover analysis results dramatically, it is important to consider at least two different lateral pushover cases to represent different sequences of response that could occur during dynamic loading [29]. Therefore, the pushover analysis in this paper should be performed using both of the following lateral load patterns: the uniform pattern, corresponding to uniform unidirectional lateral acceleration; and the 1st modal pattern. The 1st modal pattern is a pattern of forces on the joints that is proportional to the product of the fundamental mode shape times its circular frequency squared times the mass tributary to the joint. In view of the remarkable material nonlinear of concrete, the load pattern based on displacement control is chosen in this paper to avoid trouble converging. The pushover curves for the plane frame in axis ③-③ under different load patterns are shown in Figure 5.

From the simulation results in Figure 5, we can see that although the pushover curves of the frame are quite different under different lateral load patterns, the pushover curves obtained by using user-defined hinges almost exactly coincide with those by using the default plastic hinges in SAP2000 under the same lateral load pattern.

Table 5 gives the calculation results of the elastoplastic interstory drift ratios of the noncorroded plane frame under rare earthquake action. From Table 5, it could also be seen that the story drift ratios obtained by using user-defined hinges are generally consistent with those obtained by using the default hinges in SAP2000 under the same lateral load pattern. The maximum discrepancy is less than 5%, and the maximum story drift ratio is very close to 1/120, which can meet the code requirement [32].

4. Pushover Analysis Results of the Corroded Frame Structure

4.1. Plastic Hinge Properties of Corroded Beams and Columns

When considering rebar corrosion, the geometric parameters and mechanical properties of the rebar will degenerate, which will inevitably affect the plastic hinge properties of corroded beam and columns. To simplify the following computation, in this paper, the P-M-M interaction curves for different corroded RC columns are calculated using user-defined hinge with the help of XTRACT software, while the moment curvature (M3) for different corroded RC beams can be obtained using user-defined hinge by literature [35].

To understand the influence of different amounts of rebar corrosion on the seismic performance of the RC frame structure, four sets of operating conditions, including noncorroded, slight corrosion, moderate corrosion, and severe corrosion, were selected. According to the literature [35], the corresponding uniform corrosion rates of the rebar cross section under the four conditions above were 0%, 5%, 10%, and 18%, respectively. The corresponding concrete cover and reinforcement parameters under the four cases, including the strength of concrete cover, the yield strength of the rebar, the elastic modulus of the rebar, etc., are listed in Table 6.

4.1.1. Degradation of the P-M-M Curves of Corroded Column Using User-Defined Hinge

According to the parameters in Table 6, the P-M-M curves of the RC column under different degrees of corrosion can be obtained, and the results are shown in Figure 6. It can be seen that with an increased corrosion rate, the ultimate bearing capacity of compression, tension, and bending of the frame column section decreased after the rebar corrosion. However, comparatively speaking, the ultimate bending moment of the section deteriorates observably. Table 7 shows the detailed results of the degradation of maximum bending moments under cases with different corrosion rates.

4.1.2. Degradation of the Bending Strength of Corroded Frame Beam

In this paper, the bending strength for noncorroded beam hinge (M3) can be computed using the built-in default hinge properties in SAP2000. For the plastic hinge properties of corroded beams, the bending strength of the beam hinge (M3) can be obtained with the help of the calculation results of literature [35]. The ultimate bending moment of corroded beams with four different amounts of corrosion are shown in Table 8.

4.2. Pushover Analysis Results for Corroded RC Frame

Generally, the corrosion of rebar in RC frame structures is mainly caused by carbonisation, cracking, and spalling of the concrete cover. Therefore, it can be reasonably assumed that the corrosion rates of columns and beams are the same [12]. On the base of the parameters for plastic hinge properties of corroded beams and columns in Figure 6 and Tables 68, the static pushover analysis under rare earthquake action was carried out incorporating inelastic material behaviour for concrete and steel with different corrosion rates.

4.2.1. Plastic Hinges Formation and Development

Figures 7 and 8 show the distribution of plastic hinges for the corroded plane frame in axis ③-③ when the performance point (shown in Figure 3) is reached. As can be seen from the figures, with an increasing rebar corrosion rate, the number of plastic hinges in beams and columns increased gradually. Meanwhile, the hinges developed rapidly from the bottom to the higher storey with an increasing development level. Under slight corrosion condition (corrosion rate is less than 5%), the number and distribution of the plastic hinges for the corroded frame are approximately the same as that for the noncorroded frame (corrosion rate is equal to 0%). For moderately or seriously corroded frame, the number of plastic hinges increased significantly relative to the noncorroded structure, and the beam hinges have a large development rate than column hinges. By comparison, under the uniform lateral load pattern, the development degree of the plastic hinges for the same corroded frame is higher than that under the 1st modal load pattern.

As shown in Figures 7 and 8, there are 6 segments (i.e., B-IO, IO-LS, LS-CP, CP-C, C-D, and D-E) in the legend band. Each legend band is represented by a different colour, which can help the reader to judge the damage state of the structure or member. For example, it can be seen from the frame figure that there are pink points, blue points, and green points in the frame beams and columns, which are distributed in the B-IO segment, IO-LS segment, and LS-CP segment of the legend band, respectively. Therefore, according to the explanation shown in Figure 3, we can determine the stress states of these points as the serviceability stage, the safety stage, and the prevention of the collapse stage, respectively.

4.2.2. Degradation of the Seismic Performance of the Corroded Structure under Rare Earthquake Action

Figure 9 provides the pushover curves for the corroded frame under different degrees of corrosion. It is easy to see that the trend of the pushover curves and the degradation laws are basically the same despite the slight difference of the curves under different lateral load patterns. Table 9 lists the values of seismic bearing capacity of the frame under different corrosion rates. From the results in Figure 9 and Table 9, it can be concluded that the seismic bearing capacity of the corroded frame decreased significantly due to corrosion of its rebars. In addition, with increasing rebar corrosion rate, the degradation of the bearing capacity gradually increased.

4.2.3. Deformation Performance of the Corroded Structure under Rare Earthquake Action

In order to better understand the lateral deformation capacity of the corroded structure, this paper also provides the story drift ratios, as shown in Table 10. Similarly, we can see that the interstory drift ratios increase dramatically with the increasing of the corrosion rate of the structure. However, the increasing extent varies by floor position; for instance, there is a phenomenon of negative increase in the fifth and sixth floor, which is mainly caused by the redistribution of internal forces under different lateral load patterns.

Adding the interstory drift of each storey, one can get the storey displacements of the structure, which also represent the global behaviour of the corroded frame with stiffness and ductility. Figure 10 reveals the storey displacements of the corroded structure under different lateral load patterns. The results indicate that the storey displacement increases with the increasing of corrosion rate. In addition, the deformation diagram of the frame is consistent with its 1st modal shape.

5. Conclusions

In this paper, the finite element model of a six-storey-three-span corroded reinforced concrete (RC) frame structure was established using SAP2000 software. The numerical simulation of the seismic performance of the corroded RC frame structure with four different amounts of corrosion was conducted using pushover analysis method. The main conclusions are as follows:(1)For the noncorroded RC frame, the P-M-M curves of column section obtained by using user-defined hinges with the help of XTRACT software are generally consistent with those acquired by the built-in default hinges in SAP2000. For the corroded RC frame column section, the calculated P-M-M curves based on user-defined hinges present out the same change trends as the noncorroded one. Therefore, XTRACT software can be utilized to determine the properties of user-defined P-M-M plastic hinges.(2)With the increasing corrosion rate of rebars, the plastic hinges in beams and columns of the RC frame developed rapidly from the bottom to the higher storey with an increasing development level, which leads to a more aggravated deformation of the structure under severe earthquake.(3)The seismic bearing capacity of the corroded frame decreased significantly due to corrosion of its rebars. In addition, with increasing rebar corrosion rate, the degradation of the bearing capacity gradually increased.(4)The interstory drift ratios increase dramatically with the increase of the corrosion rate of the structure. However, the increase varies with by the location of storeys, and it has a strong correlation with the lateral load pattern.

Conflicts of Interest

The authors declare that there are no conflicts of interests regarding the publication of this paper.

Acknowledgments

This research has been funded by the Outstanding Young Talent Research Fund of Zhengzhou University (Grant no. 1521322004) and Guangdong Provincial Key Laboratory of Durability for Civil Engineering, Shenzhen University (Grant no. GDDCE 12-06) and Foundation for University Young Key Teacher by Henan Province (Grant no. 2015GGJS-151).