Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 292801, 8 pages

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

Research Article

Application of Base Force Element Method to Mesomechanics Analysis for Recycled Aggregate Concrete

Yijiang Peng,¹ Yinghua Liu,² Jiwei Pu,¹ and Lijuan Zhang¹

¹Key Lab of Urban Security & Disaster Engineering, Department of Civil Engineering, Beijing University of Technology, Beijing 100124, China

²AML, Department of Engineering Mechanics, Tsinghua University, Beijing 100084, China

Received 18 July 2013; Accepted 23 August 2013

Academic Editor: Song Cen

Copyright © 2013 Yijiang Peng 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 base force element method (BFEM) on potential energy principle is used to analyze recycled aggregate concrete (RAC) on mesolevel. The model of BFEM with triangular element is derived. The recycled aggregate concrete is taken as five-phase composites consisting of natural coarse aggregate, new mortar, new interfacial transition zone (ITZ), old mortar, and old ITZ on meso-level. The random aggregate model is used to simulate the mesostructure of recycled aggregate concrete. The mechanics properties of uniaxial compression and tension tests for RAC are simulated using the BFEM, respectively. The simulation results agree with the test results. This research method is a new way for investigating fracture mechanism and numerical simulation of mechanics properties for recycled aggregate concrete.

1. Introduction

“Concrete” is considered as heterogeneous composites whose mechanical performance is much related to the microstructure of material. The mechanical performances are usually obtained using experimental methods. However, testing usually consumes a large amount of manpower and material resources, and test results are usually more discrete. In order to overcome this defect, the concept of numerical concrete was presented by Wittmann et al. [1] based on micromechanics. Subsequently, some scholars did some creative works in this field and made a number of models. Among them, the two important models are the lattice model and the random aggregate model. For example, Schlangen et al. [2, 3] applied the lattice model to simulate the failure mechanism of concrete. Liu and Wang [4] adopted the random aggregate model to simulate cracking process of concrete using FEM. Peng et al. [5] adopted the random aggregate model to simulate the mechanics properties of rolled compacted concrete on meso-level using FEM. Du et al. [6] simulated the failure mechanism of beam under impact loading and triangular cyclic loading by using displacement-controlled FEM, stress-strain curves, and dynamic bending strengths of specimens.

“Recycled aggregate concrete” which is used as a green building material has attracted more and more researchers with the shortage of resources and an increasing number of construction wastes. They have carried out a series of experiments and some conclusions have been obtained. An overview of study on recycled aggregate concrete has been given by Xiao et al. [7]. However, because of the complexity of recycled coarse aggregates, conclusions made by different researchers are usually not very accordant, even opposite sometimes. To remove effects of experimental conditions, some numerical researches on meso-level were considered. For example, numerical simulation on stress-strain curve of recycled aggregate concrete was made by Xiao et al. [8] with lattice model under uniaxial compression. A method on mesomechanics analysis was proposed by Peng et al. [9, 16] and Zhou et al. [10] using FEM for recycled aggregate concrete based on random aggregate model. However, the numerical researches on the damage mechanism for recycled aggregate concrete material have just begun.

In recent years, a new type of finite element method, the base force element method (BFEM), has been developed by Peng et al. [9, 11–17] based on the concept of the base forces by Gao [18]. In this paper, the BFEM on potential energy principle is used to analyze recycled aggregate concrete (RAC) on meso-level. The mechanics properties of recycled aggregate concrete in uniaxial compression and tension tests are simulated using the BFEM, respectively. The simulation results agree with the test results. This research method is a new way for investigating fracture mechanism and numerical simulation of mechanics properties of recycled aggregate concrete.

2. Basic Equation

Consider a two-dimensional domain of solid medium; let denote the Lagrangian coordinate system, where and are the position vectors of a material point before and after deformation, respectively. Two triads for original and current configurations can be defined as follows:

Let denote the displacement of a point; then

The gradient of displacement can be written as follows:

Then, the Green strain can be written as

In order to describe the stress state at a point , a parallelogram with the edges , is shown in Figure 1. Define

where for indexes. Quantities are called the base forces at point in the two-dimensional coordinate system .

Figure 1: The base forces.

According to the definitions of various stress tensors, the relation between the base forces and various stress tensors can be given. The Cauchy stress is

Further, the base forces are given as follows:

in which is the strain energy density and is the mass density before deformation.

Equation (7) expresses the by strain energy directly. Thus, is just the conjugate variable of . It can be seen that the mechanics problem can be completely established by means of and .

3. Model of BFEM with Triangular Element

We will derive explicit expressions for stiffness matrices of a triangular element now, based on the concept of “base forces.” Consider a triangular element with boundary as shown in Figure 2.

Figure 2: A triangular element.

For the small displacement case, the real strain can be replaced by . We can obtain the average strain in element as

in which is the area of element.

Substituting (4) into (8), we have

Using Green’s theorem, (9) becomes

where is the current normal of boundary .

When the element is small enough, (10) can be written as

where is the length of edge , denotes the external normal of edge , and is the displacement of geometric center of edge .

Further, we assume that any edge of the triangular element in the deformation process keeps its edges straight lines. Then, we can obtain the following expression for :

where and denote the displacements of both ends of edge , respectively.

Substituting (12) into (11) yields

The summation rule is implied in the above equation, and is

where and are the lengths of edges and and and denote the external normals of edges and , respectively.

Then, for an isotropic material, the strain energy in the element is reduced to

in which is Young’s modulus and is Poisson’s ratio.

Substituting (13) into (15) we have

where

From (16), we can obtain the force acting on this element at node

where

Here, is a second-order tensor that is called the stiffness matrix.

The characteristics of the stiffness matrix compared with the traditional FEM are as follows. This expression of stiffness matrix can easily be extended to apply to arbitrary polygonal elements problem in two dimensions or arbitrary polyhedral element problem in three dimensions. The expression of the stiffness matrix is a precise expression, and it is not necessary to introduce Gauss’s integral for calculating the stiffness coefficient at a point. This expression of can be used for calculating the stiffness of various elements with a unified method. This expression of stiffness matrix can be used in any coordinate system. The method of constructing the stiffness matrix does not regulate the introduction of interpolation.

The model of the base force element method will be used to analyze the damage problem for recycled aggregate concrete and to analyze the relationships between mesostructure and macroscopic mechanical performance of recycled aggregate concrete.

4. Random Aggregate Model for RAC

Based on the Fuller grading curve, Walraven and Reinhardt [19] put the three-dimensional grading curve into the probability of any point located in the sectional plane of specimens, and its expression is as follows:

where is the volume percentage of aggregate volume among the specimens, in general , is the diameter of sieve pore, and is the maximum aggregate size.

According to (20), the numbers of coarse aggregate particles with various sizes can be obtained. By Monte Carlo method, random to create the centroid coordinates of all kinds of coarse aggregate particles, namely, to generate random aggregate model as Figure 3. According to the projection method, we dissect the specimens of RAC with different phases of materials. Then, the phase of recycled coarse aggregate, the phase of new hardened cement, the phase of old hardened cement, and the phase of new and old interfacial transition zone (ITZ) can be judged by a computer code.

Figure 3: Attribute recognition figure.

5. Damage Model of Materials

Components of RAC such as recycled coarse aggregate, new mortar, old mortar, new interfacial transition zone (new ITZ), and old interfacial transition zone (old ITZ) are basically quasi-brittle material, whose failure patterns are mainly brittle failure. In this paper, according to the characteristics of recycled aggregate concrete on meso-structure, the damage degradation of recycled aggregate concrete is described by the bilinear damage model, and the failure principal is the criterion of maximum tensile strain. Damage constitutive model is defined as as shown in Figure 4, where the damage factor can be expressed as follows:

where is the tensile strength of material, the residual tensile strength is defined as , the residual strength coefficient λ ranges from 0 to 1, the residual strain is , is the residual strain coefficient, the ultimate strain is defined as , where is ultimate strain coefficient, and is principal tensile strain of element.

Figure 4: Bilinear damage model.

6. Numerical Examples

According to the test results achieved from the experiments, material parameters of recycled aggregate concrete are selected for the numerical simulations. Material parameters of numerical simulations are shown in Table 1.

Table 1: Material parameters of numerical simulations.

Materials
Elastic modulus/GPa
Poisson’s ratio
Tensile strength/MPa

Natural coarse aggregate
70
0.16
10
0.1
5
10

Old ITZ
13
0.2
2
0.1
3
10

Old cement mortar
25
0.22
2.5
0.1
4
10

New ITZ
15
0.2
2
0.1
3
10

New cement mortar
30
0.22
3
0.1
4
10

6.1. Simulation on Uniaxial Compressive Strength of RAC

Based on the uniaxial compression experiments, uniaxial compressive strengths for RAC specimens are investigated using the BFEM as shown in Figure 5.

Figure 5: Numerical analysis of uniaxial compressive strength for RAC specimen.

The dimension of compression specimen is 100 mm × 100 mm × 100 mm. The numbers of large aggregates (representative size is 17.5 mm), middle aggregates (representative size is 12.5 mm), and small aggregates (representative size is 7.5 mm) are 3, 9, and 37, respectively. In order to simplify the numerical calculation, the three-dimensional structure is schematized as a plane stress problem and the geometry of the analytical specimen is 100 mm × 100 mm. The random aggregate models of a RAC specimen are shown in Figure 6. Three specimens with different random aggregate distribution are simulated and the mean compressive strength of these specimens is obtained.

Figure 6: Random aggregate models of recycled aggregate concrete with a gradation.

The uniaxial compressive stress-strain curve of recycled aggregate concrete is as shown in Figure 7.

Figure 7: Compressive stress-strain curve of recycled aggregate concrete.

The compressive strengths of the three specimens are 20.36 MPa, 20.25 MPa, and 20.05 MPa. The average uniaxial compressive strength of the specimen group is 20.22 MPa. The result of BFEM on mesodamage analysis for RAC is consistent with the test results [20].

The propagation process of cracks of the RAC specimen by uniaxial compression is shown in Figure 8.

Figure 8: Propagation process of cracks of the RAC specimen by uniaxial compression.

6.2. Simulation on Uniaxial Tensile Strength of RAC

Based on the uniaxial tension experiments, uniaxial tensile strengths for RAC specimens are investigated using the BFEM as shown in Figure 9.

Figure 9: Numerical analysis of uniaxial tensile strength for RAC specimen.

The dimension of tension specimen is 100 mm × 100 mm × 100 mm. The numbers of large aggregates (representative size is 17.5 mm), middle aggregates (representative size is 12.5 mm), and small aggregates (representative size is 7.5 mm) are 3, 9, and 37, respectively. In order to simplify the numerical calculation, the three-dimensional structure is schematized as a plane stress problem and the geometry of the analytical specimen is 100 mm × 100 mm. The random aggregate model of a RAC specimen is shown in Figure 6. Three specimens with different random aggregate distribution are simulated and the mean tensile strength of these specimens is obtained.

The uniaxial tensile stress-strain curve of recycled aggregate concrete is as shown in Figure 10.

Figure 10: Tensile stress-strain curve of recycled aggregate concrete.

The tensile strengths of the three specimens are 2.64 MPa, 2.44 MPa, and 2.63 MPa. The average uniaxial tensile strength of the specimen group is 2.57 MPa. The result of BFEM on meso-damage analysis for RAC is consistent with the test results [21].

The propagation process of cracks of the RAC specimen by uniaxial tension is shown in Figure 11.

Figure 11: Propagation process of cracks of the RAC specimen by uniaxial tension.

7. Conclusions

(1)In this paper, a model of the base force element method (BFEM) is proposed for the damage analysis problem and is used to simulate the relations of mesostructure and macro-strength of recycled aggregate concrete (RAC). The characteristics of the BFEM are that the expression of the stiffness matrix is a precise expression and it is not necessary to introduce the interpolation function and Gauss’s integral for calculating the stiffness coefficient at a point. The numerical results show that this method can be used for damage analysis of the RAC.(2)In order to simulate the meso-structure of recycled aggregate concrete material, recycled aggregate concrete is taken as five-phase composites consisting of natural coarse aggregate, new mortar, new interfacial transition zone (ITZ), old mortar, and old ITZ on meso-level in this paper. The random aggregate models are used for the numerical simulations of uniaxial compressive and tensile performances of recycled aggregate concrete, respectively. The results of the BFEM on meso-damage analysis for RAC are consistent with the test results.(3)The numerical simulation provides a new way for research on mechanical properties of recycled aggregate concrete.

Conflict of Interests

The authors declare that there is no conflict of interests.

Acknowledgment

This work is supported by the National Science Foundation of China, nos. 10972015, 11172015 and 91016026.

References

F. H. Wittmann, P. E. Roelfstra, and H. Sadouki, “Simulation and analysis of composite structures,” Materials Science and Engineering, vol. 68, no. 2, pp. 239–248, 1985.

E. Schlangen and J. G. M. van Mier, “Simple lattice model for numerical simulation of fracture of concrete materials and structures,” Materials and Structures, vol. 25, no. 9, pp. 534–542, 1992.

E. Schlangen and E. J. Garboczi, “Fracture simulations of concrete using lattice models: computational aspects,” Engineering Fracture Mechanics, vol. 57, no. 2-3, pp. 319–332, 1997.

G. T. Liu and Z. M. Wang, “Numerical simulation of fracture of concrete materials using random aggregate model,” Journal of Tsinghua University, vol. 36, no. 1, pp. 84–89, 1996 (Chinese).

Y. J. Peng, B. K. Li, and B. Liu, “Numerical simulation of meso-level mechanical properties of roller compacted concrete,” Journal of Hydraulic Engineering, vol. 32, no. 6, pp. 19–22, 2001 (Chinese).

X. L. Du, R. J. Tian, Y.-J. Peng, and Y. D. Tian, “Numerical simulation on the three-graded concrete beam under dynamic loading,” World Information on Earthquake Engineering, vol. 24, no. 1, pp. 1–5, 2008 (Chinese).

J. Z. Xiao, W. G. Li, Y. H. Fan, and X. Huang, “An overview of study on recycled aggregate concrete in China (1996–2011),” Construction and Building Materials, vol. 31, pp. 364–383, 2012.

J. Z. Xiao, J. T. Du, and Q. Liu, “Numerical simulation on stress-strain curve of recycled concrete under uniaxial compression with lattice model,” Journal of Building Materials, vol. 12, no. 5, pp. 511–514, 2009 (Chinese).

Y. J. Peng, N. N. Dang, and J. Cheng, “A method on meso-mechanics analysis for recycled aggregate concrete based on random aggregate model,” in Proceedings of the Chinese Congress of Theoretical and Applied Mechanics (CCTAM '11), pp. 1–6, 2011, (Chinese).

H. P. Zhou, Y. J. Peng, N. N. Dang, and J. W. Pu, “Numerical simulation of uniaxial compression performance for recycled concrete using micromechanics,” Applied Mechanics and Materials, vol. 253, no. 1, pp. 550–554, 2013.

Y. J. Peng and M. Jin, “New complementary finite-element method based on base forces,” Chinese Journal of Applied Mechanics, vol. 23, no. 4, pp. 649–652, 2006 (Chinese).

Y. J. Peng and M. Jin, “Application of the base forces concept in finite element method on potential energy principle,” Journal of Beijing Jiaotong University, vol. 31, no. 4, pp. 1–4, 2007 (Chinese).

Y. J. Peng and M. Jin, “A new finite element method on potential energy principle by base forces,” Journal of Beijing University of Technology, vol. 33, no. 7, pp. 687–692, 2007 (Chinese).

Y. J. Peng and M. Jin, “Finite element method for arbitrary meshes based on complementary energy principle using base forces,” Engineering Mechanics, vol. 24, no. 10, pp. 41–45, 2007 (Chinese).

Y. J. Peng and Y. H. Liu, “Base force element method of complementary energy principle for large rotation problems,” Acta Mechanica Sinica, vol. 25, no. 4, pp. 507–515, 2009.

Y. J. Peng, Z. L. Dong, B. Peng, and Y. H. Liu, “Base force element method (BFEM) on potential energy principle for elasticity problems,” International Journal of Mechanics and Materials in Design, vol. 7, no. 3, pp. 245–251, 2011.

Y. J. Peng, Z. L. Dong, B. Peng, and N. N. Zong, “The application of 2D base force element method (BFEM) to geometrically non-linear analysis,” International Journal of Non-Linear Mechanics, vol. 47, no. 3, pp. 153–161, 2012.

Y. C. Gao, “A new description of the stress state at a point with applications,” Archive of Applied Mechanics, vol. 73, no. 3-4, pp. 171–183, 2003.

J. C. Walraven and H. W. Reinhardt, “Theory and experiments on the mechanical behaviour of cracks in plain and reinforced concrete subjected to shear loading,” Heron, vol. 26, no. 1, pp. 1–68, 1981.

J. Z. Xiao, “Experimental investigation on complete stress-strain curve of recycled concrete under uniaxial loading,” Journal of Tongji University, vol. 35, no. 11, pp. 1445–1449, 2007 (Chinese).

J.-Z. Xiao and Y. Lan, “Investigation on the tensile behavior of recycled aggregate concrete,” Journal of Building Materials, vol. 9, no. 2, pp. 154–158, 2006 (Chinese).