1. Introduction
The evaluation of the cracking time and degradation of reinforced concrete cover in accelerated corrosion tests has some specific features that differ from analogous issues related to natural reinforcement corrosion [
1,
2]. These differences are mainly related to very high values of the intensity and density of the corrosion current [
1,
3]. The nature and kinetics of the process affect the composition of corrosion products and make their composition dependent on the intensity of electrical phenomena [
4]. In accelerated corrosion tests, the value of the electrochemical equivalent of iron is also different, and in the case of natural processes, it can be assumed to be analogous to that for pure iron
, given a charge number of the reaction equal to two.
From the point of view of mechanics, the process of damage evolution in the cover as a result of an accelerated corrosion process is similar to natural corrosion. As a result of the current flow, iron ions are carried into the solution, and corrosion products are formed in the first phase, tightly filling the pore spaces in the transition layer with increased porosity (ITZ interface transition zone) [
5]. In the next phase, microcracks develop within adjacent pore spaces—the so-called corrosion product accumulation area (CAR) [
6]—which influences the intensity of the impact of corrosion products on concrete cover. There are a number of different models describing the process of propagation of damage in the cover as a result of corrosion of reinforcement. A broad overview of analytical models can be found in [
7]. For computer models, the effects of corrosion products are accounted for by introducing internal volumetric strains. This type of approach can be found, among others, in papers relating to the uniform distribution of products on the circumference of the rebar [
8], non-uniform distribution [
9], modeling by imposing displacement fields on contact elements [
10], or approaches introducing a substitute tensor of volumetric strains [
11,
12] (in coupled description).
The difficulty of evaluating the propagation of damage caused by corrosion is the uncertainty of the parameters that characterize the course of the corrosion phenomena. Methods based on the probabilistic approach are not easy to implement due to the difficulty in determining the probability density function. In the group of papers on the probabilistic approach, one can distinguish items concerning the following: the impact of the probability density function on the failure probability of an element, as well as the reliability description of the safety of degrading elements under reinforcement corrosion over time [
13]; in the broader context in the cycle [
14,
15,
16], probabilistic assessment of damage propagation in the cover of reinforced concrete elements in relation to elements subjected to chloride corrosion [
17]; microbial, climatic, and fatigue effects [
18,
19]; and the application of Bayesian methods to account for uncertainties in various nonlinear material models by correlating the results of FEM calculation with experimental results [
20].
In the case of complex nonlinear calculations based on FEM in which uncertainties refer to various model parameters, such as material parameters, the only effective probabilistic computational method is a simulation approach such as the Monte Carlo method. A very in-depth analysis of this type of approach to a broad class of problems of nonlinear structural mechanics can be found in [
21]. A serious disadvantage of this type of approach is the very high computational cost associated with solving multiple problems for various randomly selected parameters.
Alternative approaches to capture uncertainty include methods based on fuzzy numbers [
22], interval arithmetic [
23], and affine numbers, which are free from the problem of overestimation of the solution [
24]. Such approaches also capture uncertainties that are not probabilistic. The application of interval methods in the mechanics of materials and structures can be found, among others, in papers concerning the use of a combined interval and probabilistic approach to the analysis of structural reliability [
25]; the formulation of a method for solving the affine systems of equations, along with applications for FEM analysis of structures, Ref. [
26]; or the formulation of a novel finite element method to limit conservatism affecting classical interval analysis in relation to linear problems of mechanics [
27]. An application of interval methods in relation to the description of volumetric strains caused by the growth of corrosion products is presented in [
28].
In this paper, the proposed approach [
28] to the evaluation of damage propagations in elements subjected to accelerated corrosion tests, taking into account affine numbers, was modified and validated. The interval algorithm for determining increments of volumetric strains was adapted to the problems of evaluating accelerated corrosion tests and further extended. The interval approach was applied to the calculation of complex, nonlinear local and nonlocal models describing the degradation of concrete test specimens. For the purpose of verifying the calculations, experimental tests were carried out on reinforced concrete elements, which were subjected to accelerated corrosion of reinforcement [
29]. The test results were used to verify the approach to estimating the width of crack opening, taking into account the uncertainty of the model parameters. The load was a time-varying increase in the interval tensor of volumetric strains caused by the impact of corrosion products on the concrete of the cover. The limits of corrosive increments of volumetric strains were calculated as follows: lower (
) and upper (
) values were determined using code written in the MATLAB language using calculation libraries that use the INTLAB interval and affine numbers [
30]. The calculations were verified using the Monte Carlo method (MC) [
31], assuming that the model parameters are random variables with a uniform distribution.
2. Purpose and Scope of the Research
As a part of laboratory tests, samples (1) with dimensions of 150 × 150 × 130 mm were made of concrete of class C50/60 (
Figure 1), in which a smooth bar (2) with a diameter of
and length of 150 mm made of St3SX steel (mark of the production period) was placed symmetrically at one edge with a cover thickness of
(samples marked with successive indices P8-P11), [
29].
The symmetrical location of the rod guaranteed the appearance of a dominant crack in the direction of the external surface of the cover, which was beneficial due to the stability of calculations using FEM.
Strong environmental aggression was forced under laboratory conditions using the so-called accelerated corrosion test based on the electrolysis of the reinforcement. The process was regulated by applying an external voltage through a voltage-current stabilizer (6) with a recorder of operating parameters. In the electrolyzer system, the anode was a bar (2) placed inside a concrete sample (1) placed in a vessel with tap water (4), while the cathode was a perforated sheet made of weathering steel (3). The concrete pore liquid served as an electrolyte. The bar and the sheet surrounding the sample were connected with copper wires (5) with a cross-section of 2.5 with positive and negative poles of the external voltage source (6). The front surfaces of the bar (protruding from the concrete) were protected against the uncontrolled flow of electric current with a cover of polyester resin.
A diagram of the experimental research carried, as well as a general outline of the theoretical research (discussed in
Section 4 and
Section 5) are presented in
Figure 2.
During the test, the stabilizer maintained a constant voltage of 20 V in the system. Through the automatic recording of current (I) and resistance (R) with a frequency of 0.0167 Hz, the history of changes in electrical parameters during the test was monitored for each analyzed sample. The electrolysis process continued until cracks with a width of 1 mm appeared, which was accepted as sufficient to stop the process of accelerated corrosion.
The course of changes in the function of the electric current and resistance of the system is shown in
Figure 3.
The adopted research and measurement system, in which the test elements were completely immersed in water, compensated for deviations resulting from differences in the moisture content of the concrete structure [
11]. This assumption made it possible to adopt a uniform distribution of corrosion products on the side surface of the reinforcing bar (2) subjected to the action of an external current source. Steel bars were subjected to gravimetric analysis. The loss of reinforcement mass was compared with the results obtained from theoretical analysis according to Faraday’s law. Before placing the bars in the molds, all bars were weighed (mass,
). After preparing the specimens and performing the accelerated corrosion test, the cracked concrete samples were split, and the bars were removed. After derusting in phosphoric acid and mechanical cleaning, the bars were weighed again (mass,
), separately obtaining the mass loss (
) for each bar. The result of the above measurements was the determination of the
parameter, which captures the difference between the rod mass loss (
) calculated analytically using Faraday’s law and the mass loss obtained from the gravimetric analysis (
) [
11]:
where
is the bar mass loss calculated using Faraday’s law, and
is the bar mass loss determined directly from the gravimetric analysis.
The effective electrochemical equivalent of the reinforcing steel determined individually for each bar placed in the test element (
) was expressed as a function depending on the electrochemical equivalent (
=
), assuming a charge number of the reaction equal to two [
11]:
where
is the electrochemical equivalent of iron.
The obtained test results are summarized in
Table 1. The final theoretical loss of bar mass adopted for further calculations (
) was determined for subsequent bars based on Faraday’s law and the average value of the effective electrochemical equivalent of the reinforcing steel (
). The obtained average differences (
) between the theoretically determined bar mass loss and the loss determined by gravimetric analysis did not exceed 1.5%.
Additionally, during tests, as a result of the increase in pressure with which the corrosion products act on the concrete cover, some of the corrosion products were pushed out of the steel–concrete contact area through the crack that formed in the cover (marked with an arrow), (
Figure 4a) [
29]. Significant losses in reinforcing steel caused by the corrosion process occurring in the place of application of the polyester resin on the front surface of the reinforcing bar immersed in water were also observed (
Figure 4b).
The influence of the corrosion products on the values of the test results was quantitatively determined by analyzing the content of corrosion products in vessels containing tap water (electrolyte) (see
Figure 1). The amount of iron compounds released during the tests was estimated based on the gravimetric analysis of the liquid in which the test samples were immersed. After the accelerated corrosion test, a volume equal to
was taken from the electrolyte solution in each vessel; then, the remaining volume of the solution (
) was measured with an accuracy of 10 mL. The selected volume (
) was filtered and dried at a temperature of more than 100 °C. The difference in the mass of the clean filter (
) before and after the test, together with the dried sludge (
) after the test, allowed us to determine the mass amount of corrosion products (
) in a volume of
of the solution.
The mass amount of iron compounds in the entire volume of the solution (
) taken beyond the volume of the test element (
) is proportional to the mass of products contained in the volume (
). In addition, mass losses of the front surfaces of the reinforcing bars (
) corresponding to the mass of corrosion products produced outside the steel–concrete contact area were estimated. The percentage share of corrosion products that mechanically affect the concrete cover (
) was determined based on the following relationship:
where
is the actual mass of all corrosion products determined gravimetrically,
is the mass of corrosion products transferred to the electrolyte solution, and
is the mass of products transferred to the solution from the front surfaces of the rod not directly affecting the concrete cover. A summary of the mass content of corrosion products transferred to the solution in relation to the samples analyzed in this paper is presented in
Table 2.
Analysis of the test results listed in
Table 2 allows one to estimate the amount of corrosion products that, as a result of the increase in the pressure of the impact of corrosion products on the cover, were pushed out of the concrete and had an impact on the propagation of a damage in the cover at the level of 50–60%. This allows for the estimation of the percentage of corrosion products affecting the cover of concrete during the electrolysis process at the level of about 40-50%. To avoid calculation problems, the average amount of corrosion products acting mechanically on concrete was assumed, which was included in the impact coefficient:
(
).
During the corrosion test, manual measurement of the length of the side of the scratched edge of each sample
was carried out cyclically along the two extreme edges (at a distance of 3 mm from the face of the sample) with an accuracy of 0.1 mm-(
Figure 5; we assume that the change in edge length is approximately equal to the average crack width). To verify this assumption, manual measurement of the width of the cracks in the element was also performed. The first (basic) measurement was performed before connecting to the electrolyzer system (
), with the next measurements conducted in the following hours of the test (
43, 116, 209, 259, and 307
.)The second measurement (
) was performed in time close to the first scratch on the surface of the sample, and the last measurement was conducted at
. The electrolyzer was disconnected from the power supply system after
The elongation of the cracked edge of each sample at successive points of time (t) was calculated as the average value (
) determined from two extreme measurements according to the following relation:
where
and
are the extensions of the sample edges in two extreme positions, while
and
denote the lengths of these edges at time t.
The history of changes in the crack width (
) and the elongation of the sample edge (
) as a function of the duration of the accelerated corrosion test are shown in
Figure 6a,b, respectively [
29].
As can be seen in the figure, the assumption used to identify the average length increment of the edge and the average crack width can be considered correct. It can also be seen that only the results obtained for the P9 and P11 test elements are similar. In the case of the P10 element, a smaller loss of iron ions from the reinforcement was observed, as shown in
Table 1. This situation may result from the poor quality of the connection between the wire and the reinforcement bar. The value of the electrochemical equivalent of iron does not differ significantly for test items P9 and P11 (
Table 1). In the case of the P8 element, a significant deviation in the displacement results obtained in the tests from the results obtained for the P9 and P11 samples can be noticed, which is caused by the appearance of additional corrosion sources located around the front surfaces of the rod as a result of strong damage of the epoxy resin protecting the front surface of the rebar.
4. Material Parameters and Analysis of Calculation Example
Calculations of the edge displacements of the test elements (ANSYS and ATENA) were performed, in which the uncertainty of the model parameters was assumed. The calculations were performed using the finite element method (FEM). The uncertainty of the model parameters was shown by declaring numerical ranges in accordance with the approach proposed in [
28]. As a way of estimation, the calculations were verified using the MC method (500 draws) for a uniform distribution of random variables [
31]. The loadings were the lower (
) and upper (sup
) limits of the coordinates of the interval tensor of volumetric strains (in the case of the MC method, the limit values of a random variable with linear distribution). The purpose of the calculations was to determine the limit increments of displacements of the sample edges as a function of experimental time. The displacement increments were assumed to be approximately equal to the average crack widths.
It was assumed that the corrosion products formed on the surface of the reinforcement are a mixture of hydroxides (
and
). Average parameters characterizing the composition of the corrosion products (
and
) were adopted in accordance with [
36]. To simplify the interval computational algorithm, it was assumed that the critical time, the value of which is, by definition, lower than the cover cracking time, is a deterministic parameter (
) equal to the cracking time of the test element.
As interval parameters, the quantities describing the chemical composition of corrosion products (
and
), porosity (
), the width of the transition layer (
), the multiplier of the corrosion products impact intensity (
), and the electrochemical equivalent of iron (
) were adopted. It was assumed that the specific density of iron ions is a deterministic parameter (
). Deviations of the model parameters (
) from the mean value were calculated according to the following relationship:
where
is the percentage deviation,
is the range,
is the mean value,
is the deviation from the mean value,
is the upper bound of the interval,
is the lower bound of the interval.
Mean values, as well as upper and lower limits of the range numbers, are presented in
Table 3, assuming the variability of the
parameter at the level of 10% and 20%. In the case of Monte Carlo calculations, the declared numerical ranges corresponded to uniformly distributed random variables. The interval defining the critical time was only considered for the MC analysis. In the case of calculations using affine (interval) numbers, the mean value of the critical time (
) was used, which allowed for the simplification of the algorithm in accordance with the comments presented in
Section 3.2.
In order to perform computer analyses, a finite element mesh was generated. Both models were realized as symmetrical in a plane perpendicular to the element axis. The models and means of support implemented in both ATENA and ANSYS are presented in
Figure 7.
In the ATENA software (
Figure 7a) standard eight-node solid elements were used. In ANSYS software (
Figure 7b), the material was modeled with solid185 and cpt215 elements, while contact interactions were modeled with conta174 and targe170 elements [
38]. Due to the use of the EDM gradient model, the mesh was densified in the area where cracks were expected. In the plane perpendicular to the rod axis, finite elements with a dimension of 1 mm were assumed. In the longitudinal direction, the model was divided into eight finite elements (two elements in ATENA software). In the ANSYS program, this type of mesh was also used in calculations using the elastic–plastic model. The total increment of volumetric strains caused by the deposition of corrosion products on the side surface of the reinforcing bar was carried out in 10 calculation steps. Loads were applied in the form of increments of the tensor of volumetric strain to a substitution ring located around the reinforcing bar in a plane perpendicular to the axis of the reinforcement (
Figure 7c). Due to the interval nature of the tensor of increment of volumetric strain, the evaluation of the impact of uncertainty on the degradation of the concrete cover was analyzed for its lower (
) and upper (
) limit values. The courses of changes of the tensor of increments of volumetric strain for the assumed uncertainties (
,
, and
) are shown in
Figure 8. The calculations were performed using three different material models of concrete cover. The elastic parameters common to the models are listed in
Table 4 and
Table 5.
The inelastic parameters are listed in tables depending on the material model and the software: (a) elastic–plastic model of material with cracking (MWR, CC3DnonLinCementitious2), ATENA program,
Table 6; (b) elastic–plastic material model with strengthening and softening (MW with HSD2), ANSYS program,
Table 7; (c) nonlocal gradient-formulated elastic microplane with degradation (EDM), ANSYS program,
Table 8; [
39,
41,
43]. The interactions of steel and concrete were described by introducing a rigid contact model, which was adopted for both the model built in the ATENA program and the models defined in the ANSYS program (
Table 9).
5. Discussion of Calculation Results
As a result of the calculations, the values of the relative displacements of points A and B (
) in the analyzed test element (average crack width) were determined. Graphical images showing the elongations of the edges were generated, taking into account uncertainties of the model parameters (
) of 0, 10, and 20%. The results obtained with model parameter uncertainty of
are presented in
Figure 9. The results oscillate around the average values of displacement increments obtained for the examined test elements. The percentage deviations of the results of computer tests of the analyzed MWR and MW models are characterized by high compliance with the obtained calculation results, with percentage discrepancy between the test results of approximately 4% and 8%, respectively. The deviation from the mean value obtained for the EDM gradient model was greater (approximately 21%;
Table 10).
A graphical image of the sample edge displacements for the assumed uncertainties at the level of
is shown in
Figure 10. The curves obtained from the limit displacement values include the results of the measurements obtained for test elements P9 and P10. In the case of calculations using the MWR model (ATENA program), the obtained boundary curves also include the results of experimental tests obtained for the P11 element. An assumption of
allowed us to achieve limit curves of displacements that include all the obtained results of experimental tests, as shown in
Figure 11.
It should be noted that the results obtained for the affine approach and the MC method are very similar, with deviations close to 10%. With larger deviations of the model parameters, the results start to differ slightly, especially in relation to those obtained for the upper end of the solution interval (
). In the case of a model based on the MWR model, the deviations of the calculation result curves (
and
) using the affine approach and MC start to differ from each other in a significant way. This is most likely due to problems with the convergence of the solution in the ATENA software for increments of total volumetric strains. Furthermore, based on the course of changes of displacements presented in
Figure 9,
Figure 10 and
Figure 11, the limit values of time (
and
) at which the limit width of the opening of the crack can be reached are presented in
Table 11. The estimated limit value (
) was accepted. Time interval (
) in which these cracks may occur, assuming
, varies in the range of
percent of the total duration of the process. In the case of
, the time interval of reaching the critical crack width change is in the range of
as a percentage of the total duration of the process.
To compare damages, the cracks in the real samples (P8–P11) and the total principal tensile strain values obtained for the FEM simulation after
of the accelerated reinforcement corrosion process were analyzed. The crack patterns obtained after examining the actual test elements are shown in
Figure 12a, while in
Figure 12b shows the courses of changes of crack evolution in subsequent samples marked by red doted lines. To verify the experimental tests, damage maps (cracks and total principal tensile strains) obtained by computer simulations are shown in
Figure 12c–f. The following maps are presented: a map of cracks obtained in the MWR model (
Figure 12c), a map of total principal tensile strains in the MWR model (
Figure 12d), a map of total principal strains in the EDM model (
Figure 12e), and a map of total principal tensile strains in the MW model (
Figure 12f).
The computer-generated maps of total principal tensile strains and cracks are similar to the damage images captured on the photographed test elements, in particular with regard to the calculations obtained for the total principal tensile strains in the MW model (
Figure 12f) and both crack distribution and total principal tensile strains for the MWR model (
Figure 12c,d). In the case of the results obtained for the gradient model, the graphical images of total principal tensile strains differ in shape and position from the results of the experimental tests.
The results of simulation tests, maps of displacements (
) in the X direction, cracks in the test element in the case of ATENA software, and total principal tensile strains at
depending on the material model and software are shown in
Figure 13,
Figure 14 and
Figure 15, with model parameter uncertainty of
and a corrosion interaction parameter of
.
Figure 13a,b show graphical images of displacement fields (
) and principal tensile strain (
), respectively. Graphical images were obtained using ANSYS software for the elastic–plastic model MW with HSD2. The displacements in this model are objective and independent of the FEM mesh. In terms of strains, the model is not objective; strains are localized in a single finite element. As previously mentioned, when analyzing the graphical images obtained in
Figure 12, the maps objectively reproduce the image of cracks in the test elements.
Figure 14 presents the results of calculations obtained with ATENA software for the two surface MWR models. The graphical images present the distribution of displacements (
,
Figure 14a), the distribution of principal tensile strain (
,
Figure 14b), and the distribution of the cracks in the test element (
Figure 14c). The maps of the obtained calculation results presented in
Figure 13 and
Figure 14 are similar in terms of both the obtained displacement and strain values. It is worth emphasizing once again that the two local elastic–plastic models are formulated similarly from a theoretical point of view.
Figure 15 shows graphical images obtained for the last series of model studies. The calculations were performed using a gradient EDM model. The following graphical images are presented: the displacement distribution (
,
Figure 15a) and the distribution of principal tensile strain (
Figure 15b). The obtained graphical images of the total principal tensile strain differ from the image of the crack’s distribution obtained in the experimental study (
Figure 12). However, it should be emphasized that the total principal tensile strain in the EDM model representing the location of cracks is objective in terms of strain value and reflects the actual strain results in the test sample (in the major crack). It can also be noted that the displacement maps shown in
Figure 13,
Figure 14 and
Figure 15 are consistent with the results of elongation (
) presented in
Figure 9; despite the discrepancies in the obtained research results, they should be considered a good approximation of the experiment.
6. Conclusions
Analysis of the calculation results indicates the high efficiency of the proposed approach in predicting the propagation of damage in the reinforced concrete elements analyzed during accelerated corrosion tests. The results obtained using the MC method and the interval approach based on affine numbers and the presented algorithm of the procedure achieve similar results in the analyzed uncertainty range. With larger (above
) model parameter uncertainties, a slight increase in discrepancies can be seen between the results obtained using the affine approach and those obtained using the MC method. The proposed approach to estimating the influence of uncertainty on the time of damage propagation in the cover can be easily extended to problems in which we deal with natural corrosion of the reinforcement. The predicted time of possible damage to the elements at which the limit width of the cracks may or may not be reached increases significantly with increasing uncertainty of the model parameters. However, taking into account that a difference in the time of occurrence of a crack width equal to 0.4 mm (equal to the increase in the length of the element) for the extreme calculation results obtained for the MWR and EDM models (
Figure 9) is about 40 h (which is about 13% of the total process time), the obtained results should be considered a reliable assessment of the time interval in which damage to the element may occur.
Analyzing the results of computer research, it should be noted that with respect to the issues in which real corrosion processes and the assessment of the prospects for damage development in a degrading reinforced concrete element take place, there may be very significant discrepancies between the observed state and that which may take place as the corrosion process develops. Even with relatively small differences in forecasts regarding the conditions under which the reinforcement corrosion process occurs, the obtained results may differ significantly.
The calculations clearly show that the proposed approach effectively reflects the impact of uncertainty on the obtained discrepancies in research results. These results are similar to those obtained using the MC method and assuming a linear distribution of random variables, as shown by tests for uncertainties oscillating within limits of 20%. In addition, it should be stated that the proposed approach seems to be a forward-looking one, especially when attempting to narrow down the obtained results, e.g., by introducing fuzzy numbers and a probabilistic approach.
Further development of the research presented in this paper may also be undertaken within the scope of implementation of the presented model for full-size structural elements operating under real conditions of an aggressive environment.