Abstract
The RC beam usually works under the combined action of external loads and corrosion expansion pressures, and the RC beams without corrosion of reinforcements have good mechanical performance. In this study, the mesoscopic models of the three-point bending RC beams with different aggregates, different aggregate contents, and different corrosion expansion pressures were built. The plastic strain and the deflection are analyzed. The results show that the cracks of the RC beam with polygon aggregates are more dispersed than circle aggregates. As the aggregate content increases, the maximum plastic strain and maximum deflection of the RC beam decrease. The aggregate contents increased from 25% to 35%, the maximum plastic strain decreased from 3.4 × 10−3 to 2.1 × 10−3, and the maximum deflection of the lower boundary of the beam decreased from 0.005 m to 0.004 m. Furthermore, the corrosion expansion pressure makes the plastic strain increase in the concrete cover.
1. Introductions
Corrosion-induced crack is the primary cause of the bearing-capacity degradation of RC beams [1, 2]. Reinforcement corrosion not only reduces the bonding performance between reinforcement and concrete [3–8], but also causes concrete cracking and even spalling of concrete cover [9–11]. In addition, RC beams are often subjected to external loads. Under the combined action of reinforcement corrosion and external loads, cracks in concrete will grow faster, which will result in structure performance degradation [12–14]. Therefore, it is very significant to analyze this problem.
In the past two decades, a tremendous amount of work has been carried out to investigate the cracking behavior of RC beams under loads and reinforcement corrosion by using experimental methods and numerical methods. In experimental studies, the galvanic method and artificial climate exposure method are used to study the crack characteristics and mechanical performance for the corroded RC beams under sustained loading [5, 15–21]. The crack-propagation behaviors of the RC beams under different loading levels and corrosion levels were described [15–18], and the deflection variation is compared in the test [5, 19–21]. The work of concrete cracking under the coupled effect of sustained loading and corrosion is quite extensive.
Comparing with the experimental study, the numerical calculation is an effective tool to analyze the internal cracking mechanism of RC beams. Chernin and Val [22] combined numerical calculation and experimental verification and discussed the initial cracking time of RC beams. A simple nonlinear finite element (FE) is then formulated and validated for predicting the time to corrosion-induced cover cracking. Zhu and Zi [23] developed a two-dimensional model that coupled with ingress of chloride ions, carbonation, electrochemical reaction, and corrosion-induced mechanical damage and studied the expansive pressure and the growth of corrosion-induced cracks of RC beams. Using the model, they can predict the nonuniform distribution of the corrosion product, the expansive pressure, and the growth of corrosion-induced cracks. Ožbolt et al. [24] built a three-dimensional model that coupled with chemo-hygro-thermo-mechanical for predicting crack patterns due to corrosion-induced damage. The comparison between numerical results and experimental evidence shows that the model is able to realistically predict experimentally observed crack pattern. However, all these efforts have been carried out based on the assumption that concrete is a homogeneous medium and the effect of external loads is not considered.
In reality, the mechanical behaviors of concrete, for instance, the macroscopic mechanical properties and the failure pattern, are closely associated with the microstructure of concrete. In the numerical simulation, the heterogeneity of concrete should be considered [25–31]. Therefore, Pacheco et al. [32] used the mesoscale lattice model to study the radial expansion behavior of concrete. Du and Jin [33] assumed that the concrete is a three-phase composite composed of circular aggregate, mortar matrix, and the interfacial transition zone (ITZ) and studied the crack pattern of microconcrete components in the direction of cross section. The present work is just an attempt to study the mesoscale cracking behavior in concrete cover. There are still some efforts that need to be further discussed [33].
Therefore, the mesoscale characteristic of the corroded RC beam due to external load will be studied in this paper and the plastic strain crack and the deflection of the RC beam will be analyzed. This paper is organized as follows. Section 2 will introduce the method of building the mesoscale concrete model. A 2D model of three-point bending RC beam will be built in Section 3. Sections 4, 5, and 6 will analyze the effect of the aggregate shape, the aggregate content, and the corrosion expansion pressure. Conclusions are given in Section 7.
2. Building Method of the Mesoscale Concrete Model
2.1. Aggregate Generation
For real aggregates, the outline shapes are very complex. In this paper, a hand-held laser scanner was used to scan the real aggregate for the 3D outline shapes. For the 2D aggregate model, the original contour points of the reconstructed aggregate will make the meshing and calculation in the finite element model highly complex. Therefore, it is necessary to simplify the original contour points.
The key control points can be extracted from the original contour points [34], such as given in Figure 1. These key control points can describe the basic outline shape of the real aggregate using as few as possible contour points, named the aggregate base of the real aggregate.

(a)

(b)
2.2. Aggregate Placement
In this study, we propose an aggregate placement algorithm based on the intrusion judgment of front polygons-rear circumcircle. The algorithm is proposed as follows. Determine the circumcircles of each aggregate in the real aggregate library and then place the aggregates from large to small according to gradation. When placing a new aggregate, determine if the new aggregate circumcircle and the placed aggregate contours overlap. The algorithm of aggregate placement effectively solves the problem of low discriminating efficiency between polygons with significant numbers of control points while simultaneously satisfying the gradation variation of the aggregates without creating large gaps between the aggregates [34]. The average particle size range of aggregate is 3 mm–16 mm. Aggregate placement is shown in Figure 2.

2.3. Aggregate Overlap Criteria
An aggregate placement process is complete when the circumcircle of the aggregate is placed into the designated area. Therefore, it is necessary to judge whether the circumcircles of the existing aggregates and the new aggregate overlap when placing, and this can be regarded as the intrusion relationship between polygons and circles. The method simplifies the judgment required between traditional polygons. Simultaneously, the principle of placing aggregates from large to small circumvents the problem of low aggregate content because of gaps between the circumscribed circles and the polygons.
When judging the location of the relationship between circles and polygons, it can be divided into three cases: (1) control point inside the circle, (2) polygonal edges intersect with circles, and (3) circle inside the polygon. For the first case, it can be immediately determined if the aggregate control point is inside the circumcircle. For the second case, it is required to judge whether the distance from the circumcircle center to the straight line is smaller than the circumcircle radius and whether the vertical point falls on the line segment formed by the adjacent control point; we call it criterion A. For the third case, calculate each angle between the circumcircle center and two adjacent vertices and then determine if the sum of these angles is equal to 360°; we call it criterion B. The flowchart of the judgment criterion is shown in Figure 3.

3. Mesoscale Model of Three-Point Bending RC Beam
3.1. Sizes of the RC Beam
Considering the calculation time and efficiency of the mesoscale concrete model, the length of the RC beam is 300 mm and the height is 50 mm. The thickness of the concrete cover is 7.5 mm. The diameter of the steel bar is 4 mm. The length of the steel bar is 275 mm. The two ends of the steel bar are 12.5 mm from each end of the RC beam, as shown in Figure 4. The concrete does not include the aggregates, aggregate shrinkage within 0.5 mm is used as the interface between the aggregate and mortar, and the interface constitutive adopts the plastic damage constitutive.

3.2. Constitutive Relation
The mesoscale concrete model will be calculated using ABAQUS. The damage plasticity model [35, 36] is always used to describe the fracture behavior of concrete materials. In this paper, the evolution model of the double-fold damage variable is adopted as the constitutive of the aggregate, the mortar, and the interface, as shown in Figure 5.

The damage variables can be expressed aswhere is the largest primary tensile strain value of the loading history. The relationship between parameters is , , and .
The elastic modulus of the aggregate is 50 GPa [37], tensile strength is 6 MPa, Poisson’s ratio is 0.2, is 0.1, is 10, and is 5. The elastic modulus of the mortar is 26 GPa, tensile strength is 3.2 MPa, Poisson’s ratio is 0.2, is 0.1, is 10, and is 4. The elastic modulus of the interface is 25 GPa, tensile strength is 3 MPa, Poisson’s ratio is 0.2, is 0.1, is 10, and is 4.
3.3. Simulation of Interaction between Steel Bar and Concrete
While formulating the rust expansion model, the following basic assumptions are made [38]: (1) the corrosion process is spatially uniform around the steel reinforcement which results in a uniform radial expansive pressure at the steel-concrete interface; (2) the concrete around the steel reinforcing bar is modeled as a thick-walled cylinder, and the wall thickness equals to the thinnest concrete cover.
In the initial stage of steel corrosion, the rust layer will fill the gap between steel and concrete, the gap layer thickness is , as shown in Figure 6. With the accumulation of rust layer, the rust expansion pressure will cause compressive deformation of the concrete; the deformation thickness of the concrete is . In Figure 6, the diameter of steel bar is indicated by D, the combined diameter of noncorroded steel plus free-expansive corrosion products is indicated by Df, the radial loss of steel is , the diameter of noncorroded steel bar , and the corresponding radial displacement is indicated by .

(a)

(b)

(c)
When R1 = D/2 + and R2 = Df/2, the deformation compatibility equation at steel-concrete interface is given by
The displacements of the concrete and corrosion products are defined as follows [39]:where is the effective elastic modulus of the concrete with , in which is the creep coefficient of the concrete [40], is Poisson’s ratio of the concrete, and and are elastic modulus and Poisson’s ratio of the corrosion products, respectively.
The relationship between corrosion rate and parameter of Dn with Df is given as follows:where n is the ratio of volume expansion of corrosion products, is the total mass of steel, and is corrosion loss mass of steel.
Then, from equation (2) using the relationship given in equations (3)–(7), the relationship between the radial pressure and the corrosion rate can be expressed as
Concrete cracking is affected by the corrosion of internal steel bar [41]. Furthermore, in view of the bonding slip between the steel bar and the concrete, the tangential springs are set per 30 mm along the length direction of the steel bar [42], as shown in Figure 7.

3.4. Boundary and Load Conditions
The left boundary is the vertical constraint and the horizontal constraint. The right boundary is the vertical constraint. The upper middle of the RC beam is loaded, as shown in Figure 4. The pressure applied to the top of the beam is 5 MPa. The number of elements is 111265.
4. RC Beam with Different Aggregates
Two kinds of aggregates are analyzed. The first one is the polygon aggregate based on the laser scan of the real aggregate, as shown in Figure 8(a). The second one is the circular aggregate, as shown in Figure 8(b). The plastic strain of the RC beam and the deflection curve of the RC beam under the external load are analyzed, respectively.

(a)

(b)
Figure 9 shows the plastic strain of the three-point bending RC beam with two different aggregates. Both of the aggregates’ content is 35%. Comparing Figure 9(a) with Figure 9(b), the maximum plastic strain of the RC beam with the circular aggregate is greater than the polygon aggregate. Because of the complex outline shapes of the real aggregates, the plastic strain cracks [42] of the RC beam are more irregular. It makes the plastic strain cracks more dispersed. In contrast, the plastic strain cracks of the RC beam with the circular aggregates are almost arc. It makes the plastic strain cracks more concentrated. In addition, there are some plastic strain cracks in the concrete cover of the RC beam using the polygon aggregates.

(a)

(b)
Figure 10 shows the deflection curve of the lower boundary of the RC beam. The vertical displacement of the lower boundary of the RC beam with the circular aggregates is much larger than the polygon aggregates. It implies that the RC beam with the polygon aggregate has better flexural rigidity.

Therefore, we can infer that the RC beam with the circular aggregate is more likely to form a penetrating large crack in the subsequent force process. In the following section, only the polygon aggregate is analyzed.
5. RC Beam with Different Aggregate Contents
There are three kinds of aggregate contents, 25%, 30%, and 35%, in this section. The plastic strain cracks of the RC beam are analyzed, as shown in Figure 11. The maximum plastic strain of the RC beam with 25% aggregate content is the largest among these three aggregate contents. With the increase in aggregate contents, the maximum plastic strain decreases. In addition, Figure 11 shows that the more aggregate content, the more dispersed the plastic strain cracks. It reflects that more aggregate content will make the RC beam have better flexural rigidity. The output is PE in Figures 11 and 12.

(a)

(b)

(c)

(a)

(b)

(c)
6. RC Beam with Different Corrosion Expansion Pressures
According to equation (8), the ultimate corrosion expansion pressure is calculated as 2.69 MPa. In this section, three different corrosion expansion pressures, 2% ultimate corrosion expansion pressure, 3% ultimate corrosion expansion pressure, and 3.5% ultimate corrosion expansion pressure, are analyzed.
Figure 12 shows the plastic strain cracks of the RC beam with 25% aggregate content under different corrosion expansion pressures. With the increase in the corrosion expansion pressure, the maximum plastic strain increases. The development of plastic strain is mainly in the concrete cover. Comparing Figures 11(a) and 12(a), it is obvious that the corrosion expansion pressure causes the large area of the plastic strain cracks in the concrete cover. The RC beam with 30% aggregate content and the RC beam with 35% aggregate content under different corrosion expansion pressures have the same pattern as Figure 12.
Figure 13(a) shows the deflection of the upper boundary and the lower boundary of the RC beam with no corrosion. From Figure 13(a), the deflection of the upper boundary is a little greater than the lower boundary. This is due to the larger elastic modulus of the steel bar than the concrete materials. In order to satisfy displacement coordination, the concrete deformation is constrained by the steel bar. Therefore, the vertical displacement of the lower boundary is smaller than the upper boundary.

(a)

(b)
Figure 13(b) shows the deflection of the upper boundary and the lower boundary of the RC beam with 2% ultimate corrosion expansion pressure. From Figure 13(b), the vertical displacement of the lower boundary is much greater than the upper boundary. This is different from the results shown in Figure 13(a). The reason for this phenomenon is that corrosion expansion pressure is applied directly on the concrete cover. It makes the vertical displacement of the concrete cover increase rapidly. This result can be explained in mesoscale through Figures 11 and 12.
Figure 14 shows the deflection of the lower boundary of the RC beam with different aggregate contents under different corrosion expansion pressures. It is obvious that, with the increase in the corrosion expansion pressure, the vertical displacement of the lower boundary of the RC beam increases. In addition, the displacement of the lower boundary of the RC beam with 35% aggregate content is the smallest among these three aggregate contents.

(a)

(b)
7. Conclusions
Based on the mesoscale concrete model, the three-point bending RC beams with different aggregate shapes, different aggregate contents, and different corrosion expansion pressures are simulated. The plastic strain crack and the deflection of the RC beam are analyzed. According to the numerical calculation, the conclusions are drawn as follows:(1)The plastic strain cracks of the RC beam with the polygon aggregates are more dispersed than the circle aggregates. It presents that the RC beam with the polygon aggregate has better flexural rigidity.(2)As the aggregate content increases, the maximum plastic strain and maximum deflection of the RC beam decrease. The aggregate contents increased from 25% to 35%, the maximum plastic strain decreased from 3.4 × 10−3 to 2.1 × 10−3, and the maximum deflection of the lower boundary of the beam decreased from 0.005 m to 0.004 m.(3)The corrosion expansion pressure makes the plastic strain increase in the concrete cover. With the increase in the corrosion expansion pressure, the deflection of the lower boundary of the beam will increase, which will lead to the peeling of the concrete cover.
Data Availability
Data will be made available on request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This study was sponsored by the National Nature Science Foundation of China (Grant no. 51908428) and the Natural Science Foundation of Shanghai (Grant no. 19ZR1460400).