1. INTRODUCTION
Among the factors for the sustainable development and economic growth of modern society are the reliability and durability of structures and infrastructure facilities, especially reinforced concrete structures. However, such structural systems are vulnerable to deterioration processes resulting from chemical degradation and physical damage, which, over time, can lead to unsatisfying structural performance under service loads or accidental actions.
In the last years, there have been significant advances in modeling, analysis and design areas relate to the deterioration in structures, besides as new approaches have been proposed to evaluate the service life of structures (Ellingwood, 2016; Biondini et al., 2017, Andrade et al., 2017).
The corrosion of reinforcement is the pathological manifestation with the highest occurrence rate in reinforced concrete structures (KARI et al., 2014). In Brazil, for example, this rate varies from 14 to 64% depending on analysis zone (Carmona et al., 1988; Dal Molin, 1988; Andrade, 1992).
A precise and computationally efficient structural modeling of corrosion deterioration is an essential for structural reliability and service life analysis in order to reduce the maintenance costs of structures; especially reinforced concrete structures (Vu et al., 2000, Rao et al., 2017).
Corrosion of steel reinforcement in concrete is an electrochemical process caused by differences in the concentrations of dissolved ions, so that part of the metal becomes cathodic and another anodic, resulting in the loss of material volume and the formation of corrosion products, which is a secondary material with a volume of 3 to 10 times greater than the initial one (Mehta et al., 2014; Geiker et al., 2016).
Concrete damage modeling due to uniform corrosion has been used phenomenological processes that segment and synthesize steel corrosion, immersed in reinforced concrete, in two stages, initiation and propagation.
The initiation phase corresponds to the period that the transport of aggressive agents (CO2) occurs in the porous matrix of concrete, resulting in reducing the pH (from approximately 12.5 to 8.5) and in the depassivation of steel reinforcement. Conversely, the propagation phase is characterized by the loss of steel mass and the formation of corrosion products, which causes cracking of concrete cover or, in advanced stages, the concrete spalling (Tuutti, 1982; Bakker, 1988; Rao et al. al., 2017).
The projects costs of reinforced concrete structures that consider only the initiation period for corrosion are not the most economical, precisely because they do not consider the maintenance costs caused during the corrosion propagation, an important consideration in the service life analysis of reinforced concrete structures (Yanaka et al., 2016).
Given the above, this work analyzes the service life of reinforced concrete structures through the coupling of two models, one responsible for the time estimation at which occurs the steel depassivation (design life - DL) and another referring to the time at which the concrete element reaches the serviceability limit state (service life - SL).
2. SERVICE LIFE OF CONCRETE STRUCTURES UNDER CORROSION
The service life, safety, reliability and risk of civil infrastructure systems have become emerging problems in recent years due to natural and human disasters, sustainability issues and global warming.
The durability management of civil infrastructure involves significant expenses and, in an era of limited public resources, requires difficult decisions to establish maintenance, rehabilitation, and replacement priorities. In this regard, the definition of service life stands out as an important concept, which serves as the basis for a holistic design approach. Structures should be designed for structural safety and maintenance for a specified period, which includes a design for durability and sustainability. With the aim of design a structure with a low-maintenance during its service life, measures must be taken in the early stages of design, and it is necessary to carry out the control during the structure service life (Ellingwood et al., 2016).
Corrosion is one of the main causes of reduced service life of reinforced concrete structures, since this involves the material loss from steel surface as a result of a chemical action. The material loss leads to an effective area reduction of cross section, and, consequently, decreased load-bearing capacity.
However, corrosion can be delayed by adopting a medium or high durability concrete (depending on the exposure environment) or by considering a suitable thickness for the concrete cover. Broomfield (2007) and Dyer (2015) report that the lower the concrete cover and quality, the greater the corrosion possibility and, consequently, its degradative effects, e.g., the cracking in concrete.
The concrete alkalinity is due to the high concentrations of soluble calcium, sodium and potassium oxides present in the concrete microscopic pores. These oxides form hydroxides, alkaline in the presence of water, creating an optimal pH condition (between 12 and 13). In this way, the concrete protects the steel from corrosion both physically by forming a protective layer for the reinforcement, and chemically through of the alkaline condition that induces the formation of a passive film on the steel surface, very dense, thin layer of oxide, and leading to a very slow rate of corrosion (Broomfield, 2007; Köliö et al., 2017).
The uniform corrosion process, carbonation corrosion, can be segmented into two phases, initiation and propagation. The diffusion process of CO2 characterizes the initiation phase. Carbon dioxide enters the concrete porous matrix reacting with the calcium hydroxides (Ca(OH)2) present in the cement paste, leading to the formation of calcium carbonate (CaCO3) (Figure 1). In the literature, this process is called carbonation, and is responsible for some alterations in carbonated concrete, e.g., permeability reduction (Neville, 1997). In addition, this process reduces the pH of concrete (approximately 12.5 to 8.5), resulting in the destruction of a chemical layer that protects steel from corrosion mechanisms (Chang et al., 2006) (Figure 1).
The carbonation reaction begins at the concrete surface progressing internally over time. As the carbonation front reaches the steel depth, the corrosion progression phase begins, with no actual damage to the structure to this point (Tuutti, 1982; Possan, 2010; Köliö et al., 2017).
In the propagation stage, for the corrosion products formation, it is necessary, firstly, the transformation of ferrous hydroxide into ferric hydroxide (1), and then, the transformation into hydrated ferric oxide, also called the corrosion product (2).
Propagation is determined by the corrosion rate (governed by oxygen availability, relative humidity and temperature) and by the capacity of concrete cover to withstand internal stresses. The unhydrated ferric oxide (Fe2O3) has a volume of about twice that of steel it replaces when fully dense. Ferric oxide upon hydration expands further, becoming porous, increasing the volume at the steel-concrete interface about six to ten times, and causing loss of effective steel area.
The concrete expansion due to formation of corrosion products results in the cracking of concrete cover, this occurs when the stresses induced by the increasing layer of corrosion products exceed the tensile strength capacity of concrete, especially in structures with low concrete covers. The cracking in concrete facilitates and accelerates the diffusion of external agents, and may cause greater damages to the concrete, e.g., the concrete spalling (Figure 2), where a whole segment of concrete cover detaches from the structure surface. (Broomfield, 2007; Köliö et al., 2017).
3. CONCRETE STRUCTURES SERVICE LIFE
3.1 Description of the implemented model for the service life analysis
In order to determine the service life of reinforced concrete structures subjected to uniform corrosion, a FORTRAN code was developed. The code possibilities the geometric non-linear structure analysis of two-dimensional fibers composites solid based on the Positional Finite Element Method (PFEM) described in Coda (2003). To represent and simulate the effects and damages from the uniform corrosion propagation phase, an adaptation was made in the code, where the useful steel area of the reinforcement is actualized in function of the corrosion time progression as given in (3) and (4).
where d(t p ) is the service steel reinforcement (mm) in function of the corrosion progression time t p (years), d 0 is the initial bar diameter (mm), η is the rate corrosion (µA/cm²), w/c is the water/cement ratio and cob is the concrete cover (mm).
To determine when occur the steel reinforcement depassivation (referent to the DL), a concrete carbonation depth predictive model based on Artificial Neural Network Model (ANN) and developed by Felix (2016), was coupled to the PFEM code. The Figure 3 present the model inputs and the ANN topology (the neural architecture) utilized to developed the concrete carbonation model.
The time corresponding to the service life (SL) was determined in this work as the instant of time in which the characteristic crack width (5) is greater than the value described on NBR 6118 (ABNT, 2014) or the maximum vertical beam displacement is greater than 250/beam length, also described on NBR 6118 (ABNT, 2014).
where w is the characteristic crack width (mm), ϕ i is the bar diameter (mm), λ i is the surface conformation coefficient of the bar (1.0 to steel plain bar, 1.4 to steel notched bar e 2.5 to steel ribbed bar), σ si is the tensile stress (kN/cm²), E si is the steel elasticity modulus (kN/cm²), fctm is the concrete mean tensile (kN/cm²) e ρri is the reinforcement ratio in relation to the surrounding concrete area.
The Figure 4 show the flowchart referent to the model code implemented to determine the service life of reinforced concrete structures subjected to the uniform.
3.2 Validation of the concrete carbonation depth predicting model
The Figure 5(a-d) presents the concrete carbonation depth in function of the time for concrete structures in four different environment conditions (sceneries). The environment conditions are detail in the Table 1. With the purpose to demonstrate the model applicability and accuracy, the results obtained by the model were compared to different analytical models (Smolczyk, 1976; Vesikari, 1988; Bob et al., 1993; EHE, 2008; Possan, 2010) and real data Mehta (2014). More details about the analytical and the numerical model can be obtained in Felix (2016).
Scenery | CO2
concentration (%) |
Relative
humidity (%) |
Rain exposure | Cement
type |
Compressive strength (MPa) |
---|---|---|---|---|---|
I | 0.04 | 70.00 | Protected | CP II - E | 30.00 |
II | 0.04 | 70.00 | Protected | CP III | 40.00 |
III | 0.04 | 65.00 | Unprotected | CP IV | 40.00 |
IV | 0.04 | 65.00 | Unprotected | CP V | 40.00 |
OBS.: In all scenarios the addition content is zero and the analysis time is 60 years.
The results show the model applicability and that it presents as an efficient numerical tool to predict the carbonation concrete depth.
3.3 Validation of the reinforced concrete corrosion model
In order to validate the numerical model implemented and proposed in this work, a reinforced concrete beam subjected to the uniform corrosion was modeled, the results are compared with the experimental results obtained by Graeff (2007). The modeled structure consist in a reinforced concrete beam with a rectangular section. The beam dimensions is specified in the Figure 6.
The concrete beam was discretized with two different mesh, where to represent the concrete matrix we used 134 triangular finite elements, and the reinforcements was represented with 120 simple bars elements.
Regarding the materials properties, the elasticity modulus of the concrete is 2600 kN/cm², the compressive strength is 2.5 kN/cm², the tensile strength adopted was 0.179 kN/cm² and the Poisson coefficient is 0.2. About the reinforcements, the elasticity modulus is 21000 kN/cm² and the tensile strength is 50 kN/cm². With the purpose to represent and modeling elastic-linear materials we used the Saint-Venant-Kirchhoff constitutive law.
A comparison between the experimental, realized by Graeff (2007), and numerical (obtained with the model proposed) results of the vertical beam displacements is shown in the Figure 7. It is possible verified that the numerical model proposed represent efficiently the beam displacement field in the elastic linear regime of the material.
With the purpose to verify the applicability of the model to the determination of the effects caused by corrosion to reinforced concrete structures, the rate of increase of the vertical displacement of the center of the beam was compared as the reinforcement suffered degradation (loss of area), with the one obtained numerically by Graeff (2007).
It is observed in Figure 8, that the difference between the responses of the two models is increasing with the evolution of the deformation of the armature, this is due to Graeff (2007) adopting in its model a nonlinear constitutive law. However, the model proposed in this work obtains displacements equivalent to the Graeff (2007), presenting 3.20% of average deviation.
3.4 Description of the analyzed structure
The structure analyzed in this work consists of a reinforced concrete beam dimensioned according to the procedures of NBR 6118 (ABNT, 2014). The Figure 9 present the values of the loads utilized to design the reinforced concrete beam and analyze the service life. In the Figure 9 the construction and geometric details have also shown. In order to analyze the concrete structure durability exposed to a moderately aggressiveness environment, the structure was dimensioned with three different concrete cover (i.e., 20, 25 and 30 mm).
The concrete beam was discretized with two different mesh, where to represent the concrete matrix we used 486 triangular finite elements, and the reinforcements was represented with 900 simple bars elements. With the purpose to represent and modeling elastic-linear materials we used the Saint-Venant-Kirchhoff constitutive law.
The data relative to the properties of the materials are presented in Table 2(a), while in Table 2(b) the concrete and conditions exposure data are presented.
Material | Property (kN/cm²) | Value |
---|---|---|
Concrete | Elasticity modulus | 2607.16 |
Poisson coefficient | 0.20 | |
Compressive strength | 3.00 | |
Tensile strength | 0.21 | |
Steel | Elasticity modulus | 21000.00 |
Tensile strength | 50.00 |
4. RESULTS
Initially, the concrete carbonation depth advance and the design life (DL) of the beams with different cover thicknesses (20, 25 e 30 mm) are presented in Figure 10. It is immediately noticeable the importance of using an adequate concrete cover due to the aggressiveness of the medium to which the concrete structure is exposed, since adequate covering increases the DL and consequently the time required for the first human interventions and structural repairs.
In Figure 11 it is possible to analyze the advance of the degradation of the reinforcement due to the action of corrosion (the graph shows the steel area of the longitudinal reinforcement). It is observed that during the corrosion initiation phase (DL) the steel area of the reinforcement is not reduced, and the section loss is more pronounced in the early years of the corrosion propagation period, in agreement with the study of Vu et al. (2000). The authors report that there is a tendency of greater loss of steel area in the initial period of corrosion because the corrosion rate of the reinforcement is higher during this period, decreasing over time with an exponential behavior.
Figure 12 shows the DL of the reinforced concrete beams analyzed in this work. The crack width limit adopted was w lim = 0.3 mm, as prescribed by NBR 6118 (ABNT, 2014) for concrete structures with class II of environmental aggression (moderate).
It is observed that the slope of the crack opening curve decreases with the increase of the concrete cover, increasing the time required for the reinforced concrete structure to reach the cracking-opening limit state (relative to the SL).
Finally, Figures 13, 14 and 15 show the longitudinal tensions and vertical displacements of the reinforced concrete beam with 20, 25 and 30 mm of concrete cover, at the moment the cracks reach their limit value (DL).
Thought the color maps presented in Figures 13, 14 and 15 it is possible to observe that the displacements of the beam decrease with the increase of the concrete cover because of the fact that there is a greater layer of concrete for the diffusion of CO2, delaying the corrosion beginning and, consequently, the degradation of the reinforcement, responsible for giving greater rigidity to the concrete beam. This result demonstrates the need to know correctly the environment to which the structure will be exposed, since the concrete elements will be designed with the required durability, guaranteeing the design life of the structures and the safety of their users.
5. CONCLUSÕES
From the analysis of the presented results regarding the validation of the developed model, it can be concluded that the formulation for the determination of the Design Life (DL) through the use of Artificial Neural Networks (ANN) and the formulation for the mechanical analysis of reinforced concrete structures under uniform corrosion based on the Finite Element Method (FEM) are consistently implemented and present consistent results when compared to literature data.
From the simulations performed with the model generated by the coupling of the MEF with the ANN, it was possible to conclude that:
An adequate concrete cover (cover which takes into account the degree of aggressiveness to which the structure is exposed) guarantees not only a longer period to lead to the depassivation of the reinforcement, but also to reach the crack-opening limit recommended by NBR 6118 (ABNT, 2014);
Consideration of the loss of steel mass as a consequence of corrosion results in the development of cracks larger than the recommendation of the NBR 6118 (ABNT, 2014), when the level of corrosion reaches 60 %;
When the corrosion rate of the reinforcement reaches 50%, there is an increase of 30% in the beam vertical displacement, regardless of the concrete cover;
By increasing the concrete cover, the maximum displacements and the deformation of the beam are reduced at the end of the DL, due to the lower degree of corrosion to which the reinforcement is subjected;
Although all the beams at the end of the DL are with the concrete in the same state of cracking, the mechanical capacity of the beam with 30 mm of concrete cover is higher than the others beams, since the reinforcement has a lower degree of corrosion and a greater area of steel, demonstrating that a concrete cover coherent with the aggressiveness of the environment to which the structure is exposed provides durability and security for its users.
The results obtained indicate the applicability of the developed code as an efficient and alternative tool for the analysis of the design life of reinforced concrete structures subjected to uniform corrosion.
In addition, the developed program demonstrates the efficiency of the use of numerical tools in the development of models for the study and simulation of complex mechanisms of the degradation of reinforced concrete structures, corroborating with the studies of the pathology of constructions.