Numerical Analysis of the Degradation Characteristics of Bearing Capacity of a Corroded Reinforced Concrete Beam

,e static load carrying capacity of a noncorroded reinforced concrete (RC) simply supported beam is numerically simulated by ABAQUS software, and the reliability of the finite element model is verified by comparing with the test results. Based on the above model, the macroscopic mechanical properties of the beam under different degrees of corrosion are calculated. In the calculation, the degradation of the bond-slip performance andmechanical properties of corroded rebars and the coupling effect on the bearing capacity and ductility degradation of the beams are considered. ,e results show that, under conditions of slight corrosion, the degradation of bond-slip performance between the rebar and concrete has no significant influence on the bearing capacity of the beam, while the degradation of the corroded rebar had a significant effect. Under moderate and severe corrosion conditions, the bearing capacity and ductility degradation caused by bond-slip are dominant in the mechanical property degradation of the beam. Overall, the macroscopic mechanical properties of the corroded beam are influenced by the coupling effect of bond-slip degradation and the mechanical property degradation of the rebar. With the increase in the corrosion rate, the bearing capacity and ductility of the beam are decreased, and its brittleness is increased.


Introduction
Among all kinds of civil engineering structures, the reinforced concrete structure is one of the most widely used structural forms.With the increase of service life, the materials of the reinforced concrete structure will degrade under the direct or indirect influence of corrosive media from the outside world, resulting in durable damage (such as surface cracks, carbonization, spalling, and corrosion of rebars).Accidents causing damage to concrete structures often occur in the case of damage resulting in durability damage, and the resulting loss is even more difficult to measure.e sum of the annual maintenance, repair, and reinforcement costs for corrosion of reinforced concrete structures in the world has exceeded 100 billion dollars [1]. According to statistics, among the many factors leading to the damage of the durability of a reinforced concrete structure, rebar corrosion is the most common.Once the rebar in the concrete is corroded, the corrosion rate of the steel is further increased.Corrosion will lead to the degradation of the geometric parameters of steel and mechanical properties, mainly reflected in the following [2][3][4]: (1) loss of mass and effective cross-sectional area; (2) the decline of nominal yield stress and ultimate stress; (3) the decrease of the nominal elongation rate and ultimate strain; and (4) the degradation of bond-slip.Degradation of these properties will inevitably weaken the static bearing capacity of the structure to a certain extent and increase its brittleness.
erefore, it is of great theoretical significance and practical engineering value to study the degradation characteristics of the bearing capacity of concrete structures or members caused by corrosion of rebars.
e previous research on the corrosion direction of RC members mainly focuses on two aspects, which are anticorrosion capability of the members and the mechanical properties of the corroded members [2][3][4].Most of the anticorrosive researches of RC members are to add some additional materials (such as polypropylene fiber [5][6][7][8][9], silica fume, and nanomaterial [10]), or some mechanical treatment schemes for members are adopted to reduce the e ect of corrosion.
At present, noncorroded RC frame structure has been widely studied in engineering [11][12][13][14].However, there are relatively few studies on corroded RC structures or members a ected by corrosion.e existing research focuses on the macroscopic mechanical properties of steel corrosion members.In most natural environments, the corrosion and degradation of reinforced concrete members will take a long time.If it is simulated under natural conditions, it will take a lot of time and cost [13].erefore, the methods of changing the temperature, concentration, and composition of the corrosive medium are used to accelerate the corrosion of the components during the test [14,15], and some scholars also accelerate the simulation process by the method of external current.In order to get enough test data, a lot of samples are needed for the accelerated corrosion test.In addition, the preparation of test specimens and related corrosion media will also cost a lot of time and cost [16][17][18], and there is few relevant impact parameters that can be used for reference during the course of experiment.rough a long period of study, it is found that the methods mentioned above are less similar to the actual corrosion in the natural environment.Based on these restrictions, the performance of corroded reinforced concrete members has not been fully studied.
With the rapid development of computer technology, numerical simulation has been widely used in engineering [19].Moreover, numerical simulation presents numerous advantages such as short computation time and low cost, and it considers the in uences of various parameters.In this paper, the characteristics of the degradation of bearing capacity of reinforced concrete simply supported beams before and after corrosion were studied using ABAQUS software to provide a reference for the performance evaluation and maintenance and reinforcement of rust structures or members.

Constitutive Relationship of Corroded Reinforced Concrete
In this paper, the concrete-damaged plasticity model (CDP model) in ABAQUS is used to simulate the mechanical behaviour of concrete materials [20].e elastoplastic model of the double inclined line is used to simulate the stressstrain relationship of rebar material.e nonlinear spring (SPRING2) is used to simulate the bond-slip between the rebar and the concrete.

Constitutive Relationship of Concrete.
e CDP model is a constitutive relation model for concrete in the implicit algorithm of ABAQUS software.e model was rst proposed by Lubliner et al. [21] and improved by Lee and Fenves [22].e model can automatically consider the biaxial or three-axis stress state of the material unit and calculate the degradation in the materials during stress under the condition of only inputting uniaxial stress-strain parameters and damage parameters of the material.erefore, the model can re ect the mechanical response of the material more truthfully in the process of monotonic loading and repeated loading/unloading [4].

Stress-Strain Relationship of Concrete Material.
Because the mechanical properties of concrete materials have strong randomness characteristics [23,24], this paper adopts the Sargin model [25] with adjustable descent velocity to de ne the compressive stress-strain relationship (Formula (1) and Figure 1).e Arab model [26] is used to de ne the tensile stress-strain relationship of concrete (Formula (2) and Figure 2).e compressive model of Sargin: where E c is the initial elastic modulus of concrete, with data selected from [27]; E g is the secant modulus of concrete; and ε c is the strain corresponding to the peak stress of concrete.
In this paper, a strain value of 0.0020 is used; f c is the peak stress of concrete, and D is the adjustment parameter of the resist compression softening section of concrete; 0 < D < 1. e tensile model of Arab:

Advances in Civil Engineering
where f t is the tensile strength of concrete, with data selected from [28]; ε t is the tensile peak strain of concrete, which can be approximated to ε t � f t /E c ; and β is the adjustment parameter of the resist tensile softening section of concrete; β > 1.
To fit the test results [29], the damage characteristics of concrete materials (3) are defined by the Lee and Fenves damage parameter model [30] with ABAQUS, which makes the calculation results closer to the experimental results: where d is the tensile/pressure damage factor of concrete; σ true is the true stress; ε pl is the plastic strain in the concrete; ε true is the real strain of concrete; ε in is the inelastic strain in the concrete; and b c/t is the scale factor between plastic strain and inelastic strain.In this paper, the tension b c/t was 0.3, and the pressure b c/t was 0.7 during calculation.Corrosion will result in the degradation of the material property of the concrete cover.In a service environment, chloride ions penetrate the surface of the rebar through the primary fractures or cracks of the concrete cover, form a conductive medium on the surface of the rebar, and participate in primary cell reflection of the corroded rebar [31,32].In the rebar corrosion state, the cracks of the concrete cover produce the stress concentration due to the extrusion force of the corroded rebar, which further increases the crack width, causing cracking and spalling.e degree and location of the cracking and exfoliation of the concrete cover have a strong randomness.To date, this process is still difficult to describe with a precise mechanical model.In this paper, the simplified formula recommended in [33,34] is used for calculation as follows: where f c−cor is the compressive peak strength of corroded concrete; c is the correlation coefficient between the surface shape of the rebar and its diameter, assuming 0.1; ε t−cor is the cracked strain of the generalized concrete; b is the primitive width of the member; b cor is the width of the cross section of the corroded member; n is the number of compressive longitudinal tendons of corrosion damage; ω cor is the total width of the corrosion crack; υ cor is the proportion coefficient of the oxidation product of corroded rebar and the volume before the corrosion of the rebar, 2.0; u icor is the crack width of the corrosion member of number i; Χ is the corrosion depth of the rebar; ρ s is the loss rate of corrosion of the rebar cross section; r is the radius before the corrosion of the rebar; and the weighted mean value is used when the diameter of the rebar is different.

Constitutive Relationship of Rebar and Degradation of Corrosion Mechanical Properties.
In this paper, the constitutive relationship of the rebar material is simulated by using a double-fold linear elastoplastic strengthening model (Figure 3).e model can fully reflect the reinforcement behaviour after the rebar yield, so it is widely used in the static analysis of the structure.
ere are two main forms of rebar corrosion [35][36][37]: (1) uniform corrosion and (2) pit corrosion.e characteristic of uniform corrosion is that the section size and shape of the rebar do not change along the longitudinal direction of the rebar, and this state is a perfect state.However, the pit corrosion is a common phenomenon in the test room or in the state of natural corrosion due to the randomness of the size and distribution of primary cracks and voids in the corrosion process of the rebar, which will cause local concentrations of chloride ions, electrons, and random pitting corrosion on the surface of the corroded rebar because the real stress, strain, and other parameters of the steel bar are Advances in Civil Engineering not changed under corrosion conditions.erefore, for the uniform corrosion, the corrosion state of the rebar can be simulated directly by the direct reduction of the section area of the rebar.However, the state of pit corrosion, considering the randomness of pit corrosion distribution, is very difficult to simulate directly using finite element software [38].erefore, we used the equivalent uniform corrosion state to simulate the pit corrosion rebar.
e Wang and Liu model [39] can be used in the nominal yield strength and nominal elastic modulus of the corroded rebar after pit corrosion: where f y is the nominal yield strength of the rebar before corrosion; f yc is the nominal yield strength of the rebar after corrosion; E s is the nominal elastic modulus of the rebar before corrosion; E sc is the nominal elastic modulus of the rebar after corrosion; and δ is the mass loss rate of the corroded rebar, and the transformation relationship between mass loss rate of the corroded rebar and the crosssectional loss rate can be seen in the following equation [39]: For the ultimate strength and ultimate strain of the pit corrosion rebar, we use the fuzzy formula of Wu and Yuan [2] to calculate as follows: where f uc is the nominal ultimate strength of the rebar after corrosion and ε uc is the nominal ultimate strain of the rebar before corrosion.

Constitutive Relationship of Bond-Slip.
e reason why the two parts of the rebar and concrete materials can better display their respective advantages is a good bond-slip relationship between rebar and concrete, which can fully coordinate the deformation of the rebar and concrete to achieve the purpose of coordination of deformation [40].
ABAQUS finite element software uses the single-axis tension and pressure function of the nonlinear spring (SPRING2) to simulate the bond-slip behaviour of the rebar and concrete.e model concentrates the slip shear stress on the surface of the rebar to the nonlinear spring element, and the relative dislocation between the rebar and the concrete is completed by pulling the rebar through the deformation of the spring [41].e mechanics calculation formula of the nonlinear spring element is where F i is the axial force of a single spring unit; A i is interaction area of a single spring unit; τ is the shear stress on the surface of the rebar; R is the radius of a single rebar at the linkage unit; and l is the spacing of the adjacent spring element.rough the research and analysis of the uniaxial tensile test curve of the corroded rebar in the literature [42], the shear stress displacement curve of the corroded rebar is found to be similar to the shear stress displacement curve defined by the European Standard CEB-FIB [43].erefore, the shear stress parameter of CEB-FIB definition ( 8) is adopted in this paper, and the bond-slip model of the corroded rebar is obtained by modifying the key point parameters of the curve (formula (9) and Figure 4): where τ u is the maximum shear stress; s is the value of the displacement; τ f is the minimum shear stress; and α is the exponent of upward section of the curve.e values of each parameter in ( 9) are shown in Table 1.e reduction coefficient of the shear stress peak value of the corroded rebar [36] can be calculated by the following equation: Because the stirrup has the effect of pulling, the slip shear stress will not be reduced to the numerical point with the increase of the slip displacement.After the corrosion, the shear stress parameter τ f will also change accordingly.Almusallam [44], Auyeung et al. [45], and other scholars believe that the final value of the longitudinal rebar corrosion τ f is 0.15 times τ u .erefore, in this paper, the research conclusions of Almusallam [44] and Auyeung et al. [45] are cited to correct corroded τ f .

The Establishment and Reliability Verification of the Finite Element Model of an Uncorroded RC Simply Supported Beam
To establish a more credible finite element analysis model, we first take a simply supported RC beam without corrosion as an example.e finite element analysis is carried out and compared with the test results to verify the reliability of the finite element model.On this basis, the bearing capacity of the simply supported RC beam with different corrosion degrees is analysed.

Introduction of Test Beam Model.
e model of the test beam is from the literature [29], and the beam number is LA6.e total length of the test beam is 2400 mm, and the spacing of the support is 2100 mm. e distribution beam applies the load to the two points at the point of trisection of the top of the beam.
e section size of the beam is 200 mm × 300 mm, and the depth of the concrete cover is 25 mm.e bottom of the beam is configured with two diameters of the 20 mm HRB335 and one diameter of the 12 mm HPB235 longitudinal rebars.e top of the beam is configured with two diameters of the 12 mm HPB235   2, and the parameters of rebar properties are shown in Table 3. e diagrams of geometric parameters and rebar parameters are shown in Figures 5 and 6, respectively.

e Establishment of a Finite Element Model.
e C3D8R element is used to simulate concrete (unit size is 25 mm × 25 mm × 100 mm).e T3D2 element is used to simulate the rebar (unit length is 100 mm).e SPRING2 element is used as the bonding element (the unit is a nonlength unit).Because the rebar is only comparatively remarkable in the longitudinal slip, the SPRING2 element nonlinearity is defined only in the longitudinal direction of the beam in this paper.e slip of the rebar is not considered in the lateral and high directions of the beam, and the spring stiffness in this direction is defined as 2 × 10 12 N/mm.
According to the weakening rule of concrete cover defined by (4), the block modelling technology is used to peel off the concrete cover from the core concrete (Figure 7), but it does not change the force and displacement transfer between nodes.To facilitate the convergence of models, two reference points (RP1 and RP2) are set up at the loading point of the beam.e distribution of the reference points and the top surface of the beam are constrained by the distributed coupling mode, and then the symmetrical displacement load is applied at the loading point.To prevent the stress concentration in the loading process, two rigid cushion plates are set up at the support, and the bottom plate and the bottom of the beam are restrained by the binding (tie) mode.e finite element model is shown in Figure 8.

Reliability Verification of the Finite Element Model.
Figure 9 shows the distribution of cracks after loading of the beams that are not corroded.e distribution of the cracks produced by the results of the finite element analysis is basically consistent with the distribution of the cracks in the test.On the whole, the cracks are symmetrical in the cross section and are distributed in a double arch.

6
Advances in Civil Engineering e load-deformation curves obtained by the finite element calculation and the test are shown in Figure 10 and Table 4.As shown in Figure 10, the finite element analysis result is very close to the test curve.e main difference is that the test model has greater stiffness, while its yield displacement and bearing capacity are relatively small.From the comparison in Table 4, the errors of all parameters are less than 10%, and all are within the acceptable range of error.
e above results indicate that the finite element model established in this paper is reliable.us, the bearing capacity analysis of the corroded members can be further carried out based on this model.

Numerical Analysis of Bearing Performance of the Corroded Reinforced Concrete Simply Supported Beam
Considering that the bottom of the beam is more easily cracked under normal conditions, the chloride ion can corrode the bottom longitudinal reinforcement through a crack in the bottom.In this paper, the bearing capacity of the beam is mainly considered when the longitudinal reinforcement is corroded at the bottom of the beam.e following 4 types of corrosion conditions are considered in the analysis: the corrosion rate distribution of the corresponding rebar section is 0% (noncorrosion), 5% (slight corrosion), 10% (medium corrosion), and 15% (severe corrosion).
e parameters used in modelling of noncorrosion corrosions are exactly the same as the parameters of the previous finite element model.e calculated parameters of the corroded member are shown in Tables 5-7 (the corrosion rates in the table are all the corrosion rates of the cross section).
In the following analysis, the degradation of bond-slip properties, the degradation of mechanical properties of the rebar, and the degradation law of macromechanical properties of the members caused by their coupling are considered.e load-displacement curves of the members with different corrosion rates under the above three conditions are shown in Figures 11-13.
Figure 11 shows that there is no significant effect on the degradation of the rebar-concrete bond-slip behaviour under slight corrosion.Under the medium corrosion condition, the degradation of the rebar-concrete bond-slip properties mainly affects the ductility of the member (deferring yield displacement and limiting the development of ultimate displacement).Under the severe corrosion condition, the degradation of bond-slip properties seriously restricts the exertion of the mechanical properties of the rebar and reduces the ductility and bearing capacity of members, resulting in brittle failure mechanism.e degradation of the Advances in Civil Engineering bond-slip properties of the rebar caused by corrosion above the middle level has a great influence on the ductility development and the mechanical properties of the rebar.e analysis of Figure 12 also shows that the bearing capacity and ductility (ultimate displacement/yield displacement) of members are gradually weakened with the increase of the corrosion rate, while the ductile failure mechanism still exists.e degradation of the mechanical properties of the rebar caused by corrosion directly affects the bearing capacity and ductility of the members, which is equivalent to decreasing the reinforcement ratio of the member section or decreasing the strength level of the reinforcement.
In Figure 13, the effects of the degradation of the bondslip properties and the degradation of mechanical properties of the rebar are taken into consideration.e bearing capacity and characteristic displacement of the beam decrease gradually with the increase in the rebar corrosion rate, and the failure mode of the beams gradually develops to brittleness.Compared with Figure 12, the bearing capacity and ductility degradation speed of members are accelerated, indicating that the mechanical property degradation and the bond-slip property degradation caused by corrosion have a strong coupling effect.Under this coupling effect, the bearing capacity and ductility of the members decrease significantly with the increase of the corrosion rate of the rebar.
To further understand the change characteristics of the bearing behaviour of the members under the coupling action of the degradation of the mechanical properties and the degradation of the bond-slip performance, the key points of load-displacement curves in Figure 13 are listed in Table 8.Advances in Civil Engineering According to the analysis of Table 8, the yield load/ displacement, ultimate load/displacement, and ductility ratio of members are all decreased to some extent after rebar corrosion, and the decrease of yield load/displacement and ultimate load is obvious.erefore, the relationship between the reduced ratio coefficient of the above three parameters (Table 9) and the corrosion rate of the rebar cross section can be fitted.e fitting results can be seen in the following equations: where I yc is the reduction percentage of the yield load; Δ yc is the reduction percentage of the yield displacement; I uc is the reduction percentage of the ultimate load; and ρ s is the corrosion rate of the cross section.

Conclusions
In this paper, the bearing capacity and ductility degradation characteristics of a reinforced concrete simply supported beam with different corrosion rates are simulated and analysed.e main conclusions are as follows: (1) Under slight corrosion conditions, the degradation of the bond-slip properties between rebar and concrete has no significant influence on the bearing capacity of the member, while the degradation of the mechanical properties of the rebar caused by corrosion has a great influence on the bearing capacity and ductility of the member.(2) With the increase in the corrosion rate of the rebar, the bearing capacity and ductility of the member decreased significantly.Under moderate and severe corrosion conditions, the bearing capacity and ductility degradation of the member caused by the degradation of the bond-slip properties are dominant.(3) On the whole, the bearing capacity of the member decreases with the increase of the rebar corrosion, and the failure mode develops from initial ductile failure to brittle failure, which is a disadvantage for the structure hit by an earthquake.

Figure 10 :
Figure 10: Comparison of load and displacement curves between test model and finite element model.

Figure 12 :
Figure 11: e load-displacement curve due to bonding degradation.

Figure 13 :
Figure 13: e load-displacement curve of beams due to corrosion.

Table 1 :
e parameters of the CEB-FIP bond-slip model.

Table 3 :
e material parameters of rebars.

Table 4 :
Comparison of key point mechanical parameters.

Table 5 :
e strength of the concrete cover of the corroded member.

Table 7 :
e degradation of bonding stress.

Table 8 :
e key point value of the load-displacement curve under coupling factors.

Table 9 :
Comparison of parameter lower rate and corrosion rate.Note: reduction percentage � (no corrosion index − corrosion index)/no corrosion index × 100%.