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

Authors: Yijiang Peng; Yinghua Liu; Jiwei Pu; Lijuan Zhang
Journal: Mathematical Problems in Engineering
 (2013)

---

## 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.


---

## Body

## 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; letxα(α=1,2) denote the Lagrangian coordinate system, where P and Q 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:
(1)Pα=∂P∂xα,Qα=∂Q∂xα.Letu denote the displacement of a point; then
(2)u=Q-P.The gradient of displacementuα can be written as follows:
(3)uα=∂u∂xα=Qα-Pα.Then, the Green strainε can be written as
(4)ε=12(ui⊗Pi+Pi⊗ui).In order to describe the stress state at a pointQ, a parallelogram with the edges dx1Q1, dx2Q2 is shown in Figure 1. Define
(5)Tα=dTαdxα+1,dxα⟶0,
where 3=1 for indexes. Quantities Tα(α=1,2) are called the base forces at point Q in the two-dimensional coordinate system xα.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(6)σ=1AQTα⊗Qα.Further, the base forces are given as follows:(7)Tα=ρAQ∂W∂uα=ρ0AP∂W∂uα,
in which W is the strain energy density and ρ0 is the mass density before deformation.Equation (7) expresses the Tα by strain energy directly. Thus, uα is just the conjugate variable of Tα. It can be seen that the mechanics problem can be completely established by means of Tα and uα.


## 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 boundaryS 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
(8)ε-=1A∫AεdA
in which A is the area of element.Substituting (4) into (8), we have
(9)ε-=12A∫A(uα⊗Pα+Pα⊗uα)dA.Using Green’s theorem, (9) becomes
(10)ε-=12A∫S(u⊗n+n⊗u)ds,
where n is the current normal of boundary S.When the element is small enough, (10) can be written as
(11)ε-=12A∑i=13Li(ui⊗ni+ni⊗ui),
where Li is the length of edge i(i=1,2,3), ni denotes the external normal of edge i(i=1,2,3), and ui is the displacement of geometric center of edge i(i=1,2,3).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 forui:
(12)ui=12(uI+uJ),
where uI and uJ denote the displacements of both ends of edge i(i=1,2,3), respectively.Substituting (12) into (11) yields
(13)ε-=12A(uI⊗mI+mI⊗uI).The summation rule is implied in the above equation, andmI is
(14)mI=12(LIJnIJ+LIKnIK),
where LIJ and LIK are the lengths of edges IJ and IK and nIJ and nIK denote the external normals of edges IJ and IK, respectively.Then, for an isotropic material, the strain energy in the element is reduced to(15)WD=AE2(1+ν)[ν1-2ν(ε-:U)2+ε-:ε-],
in which E is Young’s modulus and ν is Poisson’s ratio.Substituting (13) into (15) we have
(16)WD=E4A(1+ν)[2ν1-2ν(uI·mI)2+(uI·uJ)mIJ0+(uI·mJ)(uJ·mI)2ν1-2ν+(uI·uJ)mIJ(uI·mI)2],
where
(17)mIJ=mI·mJ.From (16), we can obtain the force acting on this element at node I(18)fI=∂WD∂uI=KIJ·uJ,
where
(19)KIJ=E2A(1+ν)[2ν1-2νmI⊗mJ+mIJU+mJ⊗mI].
Here, KIJ is a second-order tensor that is called the stiffness matrix.The characteristics of the stiffness matrixKIJ compared with the traditional FEM are as follows. (1) This expression of stiffness matrix KIJ can easily be extended to apply to arbitrary polygonal elements problem in two dimensions or arbitrary polyhedral element problem in three dimensions. (2) The expression of the stiffness matrix KIJ is a precise expression, and it is not necessary to introduce Gauss’s integral for calculating the stiffness coefficient at a point. (3) This expression of KIJ can be used for calculating the stiffness of various elements with a unified method. (4) This expression of stiffness matrix KIJ can be used in any coordinate system. (5) 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:
(20)Pc(D<D0)=Pk(1.065(D0Dmax)1/2-0.053(D0Dmax)4-0.012(D0Dmax)6-0.0045(D0Dmax)8-0.0025(D0Dmax)10),
where Pk is the volume percentage of aggregate volume among the specimens, in general Pk=0.75, D0 is the diameter of sieve pore, and Dmax 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 asE~=E(1-D) as shown in Figure 4, where the damage factor D can be expressed as follows:
(21)D={0ε<ε0,1-η-λη-1ε0ε+1-λη-1ε0<ε≤εr,1-λε0εεr<ε≤εu,1ε>εu,
where ft is the tensile strength of material, the residual tensile strength is defined as ftr=λft, the residual strength coefficient λ ranges from 0 to 1, the residual strain is εr=ηε0, η is the residual strain coefficient, the ultimate strain is defined as εu=ξε0, 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 Table1.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 Figure5.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 Figure6. 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 Figure7.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 Figure8.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 Figure9.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 Figure6. 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 Figure10.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 Figure11.Figure 11
Propagation process of cracks of the RAC specimen by uniaxial tension.


## 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 Figure5.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 Figure6. 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 Figure7.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 Figure8.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 Figure9.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 Figure6. 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 Figure10.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 Figure11.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 matrixKIJ 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.



---

*Source: 292801-2013-09-25.xml*