Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 562410, 7 pages

http://dx.doi.org/10.1155/2013/562410

## An Elastoplastic Damage Constitutive Model for Cementitious Materials under Wet-Dry Cyclic Sulfate Attack

College of Harbour, Coastal, and Offshore Engineering, Hohai University, Nanjing 210098, China

Received 7 February 2013; Revised 17 April 2013; Accepted 22 April 2013

Academic Editor: Irene Penesis

Copyright © 2013 Da Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The mechanical properties of cement mortars subjected to wet-dry cyclic sulfate attack were studied by the compression strength test. The results showed that the ultimate compressive strength increased with number of cycles at the initial stage. However, after a certain time, it started to decrease with further increases in the number of cycles. Moreover, the concentration of the sodium sulfate solution proved to be an important factor affecting the ultimate compressive strength. Based on continuum damage mechanics theory, an elastoplastic damage constitutive model is presented to describe the mechanical behavior of cementitious materials under compressive stress. The results obtained agree well with the experimentally observed elastic, plastic, and damage characteristics of cement mortars under compressive stress.

#### 1. Introduction

Sulfate attack is one of the most aggressive forms of environmental deterioration affecting the long-term durability of concrete structures [1]. As concrete suffers sulfate attack, the sulfate ions react with calcium hydroxide and calcium aluminate hydrate and form more voluminous phases [2] (gypsum and ettringite). The formation of gypsum and ettringite leads to expansion, cracking, deterioration, and disruption of concrete structures. Studying the mechanical and mineralogical properties of cement mortar with different pozzolanic compositions in a sulfate medium, Ahmet and Sükrü [1] concluded that compact structures have more effective properties against sulfate attack. In addition, these results indicated that the formation of gypsum causes cracks inside the mortar, and the formation of ettringite causes spalling on the surface of mortar. Irassar et al. [3] used X-ray diffraction to study mortar specimens (containing limestone filler) exposed for two years to sodium sulfate solution and found that the sulfate attack is characterized by an inward movement of the reaction. The first reaction in mortar was the formation of ettringite; later, gypsum deposition occurred, followed by thaumasite formation. Lee et al. [4] studied the effectiveness of silica fumes in controlling the damage arising from sulfate attack. The results showed that the presence of silica fumes had a beneficial effect on strength for mortars exposed to a sodium sulfate environment. However, mortars containing silica fumes were severely damaged in a magnesium sulfate environment. In addition, the compressive strength results indicated that the water-to-cementitious materials ratio is the most critical parameter influencing the resistance of concrete to sulfate attack. Investigating the effect of sulfate attack on the dynamic strength and deformation characteristics of concrete using the compression test, Liu et al. [5] reported that the compressive strength, elastic modulus, peak strain, and energy absorption capacity of concrete suffering from sulfate attack tend to decrease compared with control samples (concrete cured in a room environment without sulfate attack). Assuming the sulfate attack process to be a function of time incorporating diffusion, chemical reaction, and structural damage, Sarkar et al. [6] proposed a continuum damage mechanics model to assess the degradation of cementitious materials under sulfate attack and reported that the model could be used to determine ion and mineral profiles. Idiart et al. [7] presented a chemomechanical model to assess the degradation of concrete subjected to external sulfate attack. The results showed that the coupling effect (effect of cracks on the reactive transport process) is very important in the case of external sulfate attack. Cracks represent preferential penetration channels for the diffusion of sulfates into the sample. In turn, the acceleration of the diffusion-reaction process leads to a higher degree of mechanical degradation.

Concrete structures exposed to marine environments have always suffered from sulfate attack. In particular, the deterioration of concrete members placed in the tidal zone is more pronounced because of wet-dry cycles. To simulate sulfate attack on offshore concrete structures, a wet-dry cyclic sulfate attack process was examined here. The evolution of the mechanical properties of mortars suffering wet-dry cyclic sodium sulfate attack was studied. In addition, an elastoplastic damage model was proposed to describe the mechanical behavior of these mortars under compression. The constitutive model presented here could be used to assess the degradation of the mechanical properties of concrete exposed to wet-dry sulfate attack. For example, it could be used to estimate the degradation of concrete structures exposed to the tidal zone in offshore engineering.

#### 2. Experiments

##### 2.1. Materials and Experiments

Cubic samples, of size 70.7 mm and prepared with ordinary Portland cement (P.O 42.5), were used in this study. The sand-cementitious material ratio and the water-cement ratio (W/C) of the mortars were kept at 3 and 0.5, respectively. The mortar samples were removed from their molds after 24 h. Cement mortar samples were maintained for 28 days in a curing chamber, where the temperature and humidity were kept at °C and %, respectively. After this, the samples were dried in a drying oven at °C for 48 h. The dried specimens were then prepared for the wet-dry cyclic sulfate attack experiment. For each cycle, after 16 h immersed in Na_{2}SO_{4} solution, the samples were dried in air for 1 h. Then, the samples were dried at °C for 6 h and cooled at room temperature for 1 h. Afterwards, the samples were reimmersed in the sulfate solution and started a new wet-dry cycle. Solutions containing 0.5%, 5%, and 10% sodium sulfate were chosen for the wet-dry cycle tests. After different wet-dry cycle times, compression tests were conducted with a universal testing machine at a rate of 0.02 mm/min.

##### 2.2. Results

According to the results shown in Figure 1, the stress-strain curves of the mortars started with a linear rampup but then showed loss of linearity. Liu et al. [5] had reported the nonlinear zone (before reaching ultimate strength) and attributed it to the formation of interior cracks in the mortars. After reaching peak stress, the compressive stress decreased rapidly, and visible cracks formed. As the strain further increased, the mortar fractured.

As shown by the results for ultimate compressive strength (Figure 2), the strength increased with number of cycles in the early stages. However, after a certain time, it began to decrease with further increases in the number of cycles. Moreover, the transition time to ultimate compressive strength decreased with increasing sulfate concentration. For example, for mortars exposed to 0.5% sodium sulfate solution, the transition occurred after 90 wet-dry cycles. However, the transition occurred after only 60 cycles for mortars exposed to 10% sodium sulfate solution. Moreover, variations in compression strength were more significant for the samples exposed to the higher concentration sulfate solution.

#### 3. Numerical Modeling of Sulfate Attack

##### 3.1. Mechanical Damage

In the present study, an elastoplastic damage model is proposed. For simplicity in applying the model in the engineering domain, the damage state of materials is assumed to be isotropic. The bulk modulus and the shear modulus can fully characterize the elastic behavior of isotropic materials. Cementitious materials in structures suffer mainly from compressive stress. Under compressive stress, the failure mechanism of materials presents mainly in the form of shear failure. Compared with the variance in shear modulus, the change of bulk modulus under compressive stress is small. Therefore, similar to Lemaitre’s model [8], it is assumed here that only the parameter is affected by material damage. The effect of the mechanical damage variable on can be expressed as follows: where is the shear modulus of undamaged materials. As a result, the thermodynamic potential of a damaged material can be expressed as where is the elastic volumetric strain and is the components of an elastic deviation strain tensor. Therefore, the stress-strain relation of the damaged material can be expressed as follows: where is the effective stiffness tensor. Consequently, the stress increment can be expressed as In addition, the conjugate force associated with the damage variable is defined by

Note that the theory expressed above is based on the assumption that the materials maintain an elastic state. In fact, plastic deformation occurs in mortar during the compressive process. Therefore, a heuristic approach has been proposed to express the effect of plastic strain on damage evolution. Moreover, damage evolution was thought to depend separately on elastic and plastic strains. To express this relationship, the driving force was introduced to explain damage evolution and can be defined as follows: where is the components of the plastic deviation strain tensor. Including the plastic deformation effects, the stress increment can be expressed as follows: To assess the mechanism of material damage, an exponential function [9] was introduced to describe damage evolution: The parameter defines a critical value of damage corresponding to the residual strength of damaged material. controls the kinetics of damage evolution, and defines the threshold of damage force. Therefore, the mechanical damage variable can be expressed as follows:

##### 3.2. Chemical Damage

Under sulfate attack, ettringite and gypsum minerals form in mortar [1]. As the formed sulfate minerals exceed a certain amount, they can cause partial breakage of the mortar. Furthermore, the breakage can be defined by the degradation of compressive strength: where is the chemical damage variable, and are the ultimate compressive strengths of the undamaged samples, and the sulfate-attacked samples respectively, and defines a critical value of damage corresponding to the residual strength of the chemically damaged material.

The compressive strengths of the samples are closely related to the concentration of the sulfate solution and the number of wet-dry cycles . Therefore, the compressive strength can be expressed as follows: where is a function of the sulfate concentration and the number of wet-dry cycles . Then, combining (10) and (11), the following expression for can be written:

Xu [10] has reported that there is a power function relationship between compressive strength and sulfate attack time. In addition, a linear relationship between compressive strength and sulfate concentration has also been determined. Based on experimental data, compressive strength can be expressed as a function of sulfate concentration and sulfate attack time as follows:

It has been assumed in this study that chemical damage to materials can be described by thermodynamic theory. The scalar damage variable can be physically identified as a microcrack density [9]. For concrete materials suffering from sulfate attack (chemical damage), the formed gypsum and ettringite lead to the formation of internal cracks in the concrete. Therefore, the mechanical properties of concrete are affected by chemical damage. In fact, the intrinsic characteristics of chemical damage are the same as those of mechanical damage. Both mechanical and chemical damage appear incrementally during the damages process. Therefore, the total damage to materials, which includes mechanical damage and chemical damage , can be expressed as follows:

##### 3.3. Plasticity Characteristics

Small strains and isothermal conditions were assumed in the present study. The total strain tensor is composed of an elastic part and a plastic part as follows: During the microcrack extension process, local plastic deformation of mortar induced stress redistribution. To account for this redistribution of stresses, an effective stress tensor for damaged material has been introduced and is defined as follows [10]: where is the total damage to materials and and represent the nominal mean stress and the components of the nominal deviation stress tensor, respectively. Compared with the classical model used in continuum damage mechanics [10], the proposed model presented in (16) is a more general form.

The classic Mohr-Coulomb criterion [11] based on linear rules was used in the present study to describe the failure surface of the mortar: where is the mean stress, is the equivalent shear stress, and ,, , and , where is the damage deviation stress tensor. The parameter is the inverse of the slope of the failure surface, and is related to the ultimate tensile strength of the material.

Generally, the failure surface is considered as the ultimate state of a yield surface. According to the plastic hardening law, the yield surface can be defined as follows: where . The scalar-valued function defines the plastic hardening law of materials, which can be expressed as: where ,,, and the parameter controls the kinetics of plastic hardening. Note that the failure surface is created when .

When subjected to compressive stress, the volume of the mortars shrank in the initial stage. According to the generation and propagation of cracks in mortars, there would normally be a transition from compressibility to dilatancy as the applied stress further increased. To describe this transition, the nonassociated plastic flow rule proposed by Pietruszczak et al. [12] has been used here: where , defines the intersection of the plastic potential surface with the axis, and is the slope of the transition line between compressibility and dilatancy.

The plastic constitutive model of concrete can be expressed as follows: where is the plastic factor. When , according to the plastic compatibility conditions, . Assuming that damage evolution can be ignored , can be expressed as follows: where , , and is the effective elastic stiffness tensor.

##### 3.4. Computational Procedure

The iterative method consisting of elastic prediction and plastic correction was chosen for numerical computation in the proposed elastoplastic damage model. In the numerical computation process, the total damage to materials is made up of mechanical damage and chemical damage. The chemical damage variable , which is related simply to the concentration of the sodium sulfate solution and the number of wet-dry cycles, can be calculated as . The elastic materials parameters ( and ) before mechanical damage can also be calculated from the chemical damage variable . The mechanical damage is simply related to the total strain, which is known. As an example, the th computational step can be described as follows.(1) After the -st step, the parameters, including the stress , the total strain , the plastic strain , the plastic hardening variable , and the damage variable , have been obtained. Then, the strain increment is input to the th step.(2) The total strain is calculated as . Based on the proposed damage criterion, the mechanical damage variables are computed, and the total damage is obtained using .(3) Based on the assumption that the strain increment is fully elastic, the stress increment and the total stress are calculated as and respectively.(4) The function is used to evaluate whether plastic yielding has occurred.

If , plastic yielding has occurred, and the plastic correction must be carried out. An iterative calculation is used to determine the corresponding plastic strain increment . Then, the increment of the plastic hardening parameter, , is also calculated.

If the function is less than zero, the plastic yielding has not occurred. Update the corresponding variables, and go to the th step.

#### 4. Numerical Simulations

Using sample subjected to different sulfate solution concentrations and numbers of wet-dry cycles, the parameters involved in the constitutive model were determined by compression tests; the results are presented in Table 1.

Figure 3 presents the simulated compressive strength of mortars in 5% sulfate solution based on the proposed elastoplastic damage model. As the results show, similary to the test data presented in Figure 1(b), the simulated ultimate compressive strength and elastic modulus increased with the number of wet-dry cycles in the initial stage. Then, these values began to decrease as the number of wet-dry cycles further increased. Similarly, in axial compressive experiments on concrete specimens, Liu et al. [5] reported that sulfate attack led to decreases in elastic modulus and peak strain. On the other hand, as shown by the simulated results presented in Figure 4, the concentration of the sodium sulfate solution appears to have a significant effect on the mechanical properties of mortars, which is consistent with the test results. Moreover, as shown in Figure 5, the simulated results fit well with the stress-strain curve of the measured values.

Figure 6 presents a comparison between the simulated ultimate compressive strength values and test data for different sulfate solution concentrations and numbers of wet-dry cycles. As the results show, the proposed elastoplastic damage model fitted the test data very well. This indicates that the model is suitable for simulating the mechanical properties of sulfate-attacked mortar subjected to compressive stress.

#### 5. Conclusions

In this study, the mechanical properties of cement mortar in a wet-dry cyclic sulfate attack environment have been investigated. The results showed that the ultimate compressive strength of mortars increased with the number of wet-dry cycles in the early stage. Then, the ultimate strength began to decrease with further increases in the number of wet-dry cycles. In addition, the transition to ultimate strength occurred earlier for samples exposed to higher sulfate solution concentrations.

An elastoplastic damage model was proposed to simulate the mechanical behavior of sulfate-attacked cement mortar subjected to compression stress. The results showed generally good agreement between the simulated values and test data. Besides accurately simulating the effects of sulfate concentration and numbers of wet-dry cycles, the proposed model could also capture the main features of observed mechanical behavior. The mechanical behavior of cementitious materials is closely related to the composition of the materials, the water-cement ratio, and other factors. To apply the proposed constitutive model in an engineering environment, the model parameters must be determined from compressive experiments. Once the corresponding parameters have been brought into the constitutive model, the damage progression of the mechanical properties of concrete exposed to a certain sulfate concentration can also be evaluated.

#### Acknowledgments

The authors are grateful to the National Natural Science Foundation of China (Grant no. 51009061) and the Key Project of the National Natural Science Foundation of Jiangsu Province (BK2011026) for their support of this work.

#### References

- C. Ahmet and Y. Sükrü, “Investigation of mechanical and mineralogical properties of mortars subjected to sulfate,”
*Construction and Building Materials*, vol. 24, no. 11, pp. 2231–2242, 2010. View at Publisher · View at Google Scholar - R. D. Hooton, “Influence of silica fume replacement of cement on physical properties and resistance to sulfate attack, freezing and thawing, and alkali-silica reactivity,”
*ACI Materials Journal*, vol. 90, no. 2, pp. 143–151, 1993. View at Google Scholar · View at Scopus - E. F. Irassar, V. L. Bonavetti, and M. González, “Microstructural study of sulfate attack on ordinary and limestone Portland cements at ambient temperature,”
*Cement and Concrete Research*, vol. 33, no. 1, pp. 31–41, 2003. View at Publisher · View at Google Scholar · View at Scopus - S. T. Lee, H. Y. Moon, and R. N. Swamy, “Sulfate attack and role of silica fume in resisting strength loss,”
*Cement and Concrete Composites*, vol. 27, no. 1, pp. 65–76, 2005. View at Publisher · View at Google Scholar · View at Scopus - T. J. Liu, D. J. Zou, J. Teng, and G. L. Yan, “The influence of sulfate attack on the dynamic properties of concrete column,”
*Construction and Building Materials*, vol. 28, no. 1, pp. 201–207, 2012. View at Publisher · View at Google Scholar - S. Sarkar, S. Mahadevan, J. C. L. Meeussen, H. van der Sloot, and D. S. Kosson, “Numerical simulation of cementitious materials degradation under external sulfate attack,”
*Cement and Concrete Composites*, vol. 32, no. 3, pp. 241–252, 2010. View at Publisher · View at Google Scholar · View at Scopus - A. E. Idiart, C. M. López, and I. Carol, “Chemo-mechanical analysis of concrete cracking and degradation due to external sulfate attack: a meso-scale model,”
*Cement and Concrete Composites*, vol. 33, no. 3, pp. 411–423, 2011. View at Publisher · View at Google Scholar · View at Scopus - J. Lemaitre,
*A Course on Damage Mechanics*, Springer, Berlin, Germany, 2nd edition, 1992. - D. Chen, I. Yurtdas, N. Burlion, and J. F. Shao, “Elastoplastic damage behavior of a mortar subjected to compression and desiccation,”
*Journal of Engineering Mechanics*, vol. 133, no. 4, pp. 464–472, 2007. View at Publisher · View at Google Scholar · View at Scopus - H. Xu,
*Research on damage behaviour of concrete under sulfate corrosion [Doctoral thesis]*, China University of Mining and Technology, 2012. - D. Sfer, I. Carol, R. Gettu, and G. Etse, “Study of the behavior of concrete under triaxial compression,”
*Journal of Engineering Mechanics*, vol. 128, no. 2, pp. 156–163, 2002. View at Publisher · View at Google Scholar · View at Scopus - S. Pietruszczak, J. Jiang, and F. A. Mirza, “An elastoplastic constitutive model for concrete,”
*International Journal of Solids and Structures*, vol. 24, no. 7, pp. 705–722, 1988. View at Google Scholar · View at Zentralblatt MATH · View at Scopus