Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 528590, 11 pages

http://dx.doi.org/10.1155/2014/528590

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

The Key Laboratory of Urban Security and Disaster Engineering, Ministry of Education, Beijing University of Technology, Beijing 100124, China

Received 22 June 2014; Accepted 2 July 2014; Published 4 August 2014

Academic Editor: Song Cen

Copyright © 2014 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 damage mechanics is used to analyze the compressive strength, the size effects of compressive strength, and fracture process of concrete at mesolevel. The concrete is taken as three-phase composites consisting of coarse aggregate, hardened cement mortar, and interfacial transition zone (ITZ) on mesolevel. The random aggregate model is used to simulate the mesostructure of concrete. The mechanical properties and fracture process of concrete under uniaxial compression loading are simulated using the BFEM on damage mechanics. The simulation results agree with the test results. This analysis method is the new way for investigating fracture mechanism and numerical simulation of mechanical properties for concrete.

#### 1. Introduction

Concrete is considered as heterogeneous composites whose mechanical performance is much related with the microstructure of material. The composite behavior of concrete is exceedingly complex and up to now many details such as strain softening, microcrack propagation, failure mechanisms, and size effects are still far from being fully understood. Since it is difficult to look inside concrete to observe the actual crack propagation or to experimentally determine the microscopic stress field, it has become obvious that further progress based exclusively on experimental studies will be limited [1].

In order to overcome this defect, the concept of numerical concrete was presented by Wittmann et al. [2] 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. [3, 4] applied the lattice model to simulate the failure mechanism of concrete. Liu and Wang [5] adopted the random aggregate model to simulate cracking process of concrete using FEM. Peng et al. [6] adopted the random aggregate model to simulate the mechanical properties of rolled compacted concrete on mesolevel using FEM. Du et al. [7–11] researched a mesoelement equivalent method for the simulation of macromechanical properties of concrete and a mesoscale analysis method for the simulation of nonlinear damage and failure behavior of reinforced concrete members.

The finite element method (FEM) is one of the most important numerical methods developed from 1950s, and it has been the most popular and widely used numerical analysis tool for problems in engineering and science. Over the past 50 years, numerous efforts techniques have been proposed for developing finite element models [12–15] and some other improvement and alternative methods have been proposed and developed, such as boundary element methods [16, 17] and meshless methods [18, 19]. In recent years, a new type of finite element method, the base force element method (BFEM), has been developed by Peng et al. [20–25] based on the concept of the base forces by Gao [26]. Further, the base force element method (BFEM) on potential energy principle was used to analyze recycled aggregate concrete on mesolevel [27].

In this paper, the base force element method (BFEM) on damage mechanics is used to analyze the size effect on compressive strength for concrete at mesolevel.

#### 2. Basic Formula of Base Force Theory

##### 2.1. Base Forces

Consider a two-dimensional material domain. denote the Lagrangian coordinates. Let denote the position vector of a point after deformation. The tangent basis vectors can be written in the form

The position vector of a point can be written as

In order to describe the stress state at a point , a parallelogram with the edges is shown in Figure 1. Let denote the force acting on the edge. We calculate the limit where we promise for indexes. Quantities are called the base forces at point in the two-dimensional coordinate system .

In order to further explain the meaning of , let us compare with the stress vector which represent the forces per unit area in the deformed body. That is, where is called base area .

The base forces can also be understood as stress flux.

According to the definitions of various stress tensors, the relation between the base forces and various stress tensors can be given. For example, the Cauchy stress is where is the dyadic symbol and the summation rule is implied.

When the body force per unit mass is zero, the equilibrium equation can be written as

##### 2.2. Conjugate Variable of Base Forces

Let denote the position vector of a point before deformation; the displacement of a point is

The gradient of displacement can be written as

Further, the elastic law can be given as follows: in which is the strain energy density and is the mass density after deformation.

Equation (9) 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 .

For the small deformation case, the Green strain can be written as where is the conjugate of .

#### 3. Model of the BFEM

##### 3.1. Stiffness Matrix

Consider a triangular element as shown in Figure 2; the stiffness matrix of a base force element can be obtained [27] as in which is Young’s modulus, is Poisson’s ratio, is the area of an element, is the unit tensor , , and can be calculated from where and are the lengths of edges and of an element and and denote the external normals of edges and , respectively.

Further, the stiffness matrix of a base force element can be written as where the summation rule is implied.

For the Cartesian coordinate system, we introduce the unit basic vectors and have or

Substitute (15) into (13) and consider , ; we can obtain or

For the plane stress problem, it is necessary to replace by and by in (16) and (17).

Further, the stiffness matrix of a base force element for the plane stress problem can be written as in which and can be calculated from or

##### 3.2. Strain Tensor of an Element

When the element is small enough, the real strain can be replaced by the average strain . We can obtain the average strain in the element as in which is the area of an element.

Substituting (10) into (21), we can easily derive the following equation: where is the local code number for the nodes of a triangular element and the summation rule is implied.

Substitute and into (22); we can obtain the explicit expression of the tensor formula as or

##### 3.3. Stress Tensor of an Element

When the element is small enough, the real stress can be replaced by the average stress . According to the generalized Hooke law, the stress component expressions of an element can be obtained for the plane stress problem as follows:

For the plane strain problem, it is necessary to replace by and by in (25).

#### 4. Random Aggregate Model

Based on the Fuller grading curve, Walraven and Reinhardt [28] put the three-dimensional grading curve into the probability of any point which 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 (26), 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.

According to the projection method, we dissect the specimens of concrete with different phases of materials. Then, the phase of coarse aggregate, the phase of hardened cement, and the phase of interfacial transition zone (ITZ) can be judged by a computer code as in Figure 3.

#### 5. Damage Model of Materials

Components of concrete such as coarse aggregate, cement mortar, and interfacial transition zone (ITZ) are basically quasibrittle material, whose failure patterns are mainly brittle failure.

In this paper, according to the characteristics of concrete on mesostructure, the damage degradation of concrete is described by the bilinear damage model, and the failure principal is the criterion of maximum tensile stress. Damage constitutive model is defined as , shown in Figure 4, where is the initial Young modulus; 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 peak strain, is the residual strain coefficient, the ultimate strain is defined as , where is ultimate strain coefficient, and is principal tensile strain of element.

#### 6. Numerical Example

According to the test results obtained from the experiment, material parameters of concrete are selected, which is completely coherent with model material used in the experiment. Material parameters of numerical simulation are shown in Table 1. The concrete specimens were loaded by displacement steps.

##### 6.1. Compressive Strength of Concrete Specimen

For the size of compression specimen, the numbers of coarse aggregate particles can be obtained according to (26). The random aggregate model is generated as in Figure 5.

The uniaxial compressive stress-strain curve of concrete is got as shown in Figure 6. The compressive strengths of the three specimens were 24.28 MPa, 24.40 MPa, and 24.41 MPa. The uniaxial compressive strength average of the specimen group is 24.36 MPa. The result BFEM on mesodamage analysis for concrete is consistent with the test results [29].

##### 6.2. Size Effects of Compressive Strength

In order to study the size effects of compressive strength for concrete specimen, we used three different sizes of specimens in the numerical simulations, which are , , and .

###### 6.2.1. Compressive Strength of 100 mm × 100 mm × 100 mm Specimens

For the size of compression specimen, the numbers of coarse aggregate particles can be obtained according to (26). The random aggregate model is as in Figure 7.

The uniaxial compressive stress-strain curve of concrete is got as shown in Figure 8. The compressive strengths of the four specimens were 25.85 MPa, 25.33 MPa, and 25.67 MPa. The uniaxial compressive strength average of the specimen group is 25.62 MPa.

###### 6.2.2. Compressive Strength of 300 mm × 300 mm × 300 mm Specimens

For the size of compression specimen, the numbers of coarse aggregate particles can be obtained according to (26). The generated random aggregate model is as in Figure 9.

The uniaxial compressive stress-strain curve of concrete is got as shown in Figure 10. The compressive strengths of the four specimens were 22.59 MPa, 22.72 MPa, and 22.70 MPa. The uniaxial compressive strength average of the specimen group is 22.67 MPa.

The results of specimens are shown in Figure 6. The size effects of mechanical properties of concrete under uniaxial compression loading are shown in Table 2.

##### 6.3. Fracture Process of Concrete

In the calculation, we use the incremental method and the displacement loading steps. According to the first strength theory, the element causes damage when the first principal stress exceeds the allowable value of the element. All the damaged elements were represented by the black color.

The propagation process of cracks of the specimen of concrete with by uniaxial compression is shown in Figure 11.

The propagation process of cracks of the specimen with by uniaxial compression is shown in Figure 12.

The propagation process of cracks of the specimen with by uniaxial compression is shown in Figure 13.

#### 7. Conclusions

In this paper, the base force element method (BFEM) on damage mechanics is used to analyze the size effect on compressive strength for concrete at mesolevel. The characteristics of the BFEM compared with the traditional FEM are as follows. The expression of the stiffness matrix is a precise expression, and it is not necessary to introduce Gauss’ integral for calculating the stiffness coefficient at a point. This expression of stiffness matrix 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 was used to analyze the damage problem for concrete and was used to analyze the relationships of mesostructure and macroscopic mechanical performance of concrete in this paper. The concrete is taken as three-phase composites consisting of coarse aggregate, hardened cement mortar, and interfacial transition zone (ITZ) on mesolevel. The random aggregate model is used for the numerical simulation of uniaxial compressive performance of concrete. The results by the BFEM show that the uniaxial compressive strengths of specimens are approximately coincident with the experiment results and the size effect of specimens is in agreement with the common rule. The following conclusions can be drawn.(1)It shows that the BFEM with the mesodamage model is feasible and effective to study failure process and mechanical parameter of concrete.(2)When specimens are under the same conditions, the compressive strength of concrete decreases with increasing specimen size.(3)Macrofailure of concrete apparently lags behind the growth of microdamage; that is, macrofailure can be considered as the accumulation of microdamage.(4)The interface is the weakest part of concrete, and the selection of mechanical parameters of interface plays an important role on the numerical results of uniaxial compression for specimens. More future studies should be done on researching the mechanical property of the interface.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This work is supported by the National Science Foundation of China nos. 10972015 and 11172015 and the preexploration project of Key Laboratory of Urban Security and Disaster Engineering, Ministry of Education, Beijing University of Technology, no. USDE201404.

#### References

- Z. M. Wang, A. K. H. Kwan, and H. C. Chan, “Mesoscopic study of concrete I: generation of random aggregate structure and finite element mesh,”
*Computers and Structures*, vol. 70, no. 5, pp. 533–544, 1999. View at Google Scholar · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Scopus - G. T. Liu and Z. M. Wang, “Numerical simulation study of fracture of concrete materials using random aggregate model,”
*Journal of Tsinghua University*, vol. 36, no. 1, pp. 84–89, 1996, (Chinese). View at Google Scholar · View at Scopus - 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). View at Google Scholar · View at Scopus - X. L. Du, L. Jin, and G. W. Ma, “Meso-element equivalent method for the simulation of macro mechanical properties of concrete,”
*International Journal of Damage Mechanics*, vol. 22, no. 5, pp. 617–642, 2013. View at Publisher · View at Google Scholar · View at Scopus - X. L. Du, L. Jin, and G. W. Ma, “Macroscopic effective mechanical properties of porous dry concrete,”
*Cement and Concrete Research*, vol. 44, pp. 87–96, 2013. View at Publisher · View at Google Scholar · View at Scopus - L. Jin, X. L. Du, and G. W. Ma, “Macroscopic effective moduli and tensile strength of saturated concrete,”
*Cement and Concrete Research*, vol. 42, no. 12, pp. 1590–1600, 2012. View at Publisher · View at Google Scholar · View at Scopus - X. L. Du and L. Jin, “Meso-scale numerical investigation on cracking of cover concrete induced by corrosion of reinforcing steel,”
*Engineering Failure Analysis*, vol. 39, pp. 21–33, 2014. View at Google Scholar - X. L. Du, L. Jin, and G. W. Ma, “A meso-scale analysis method for the simulation of nonlinear damage and failure behavior of reinforced concrete members,”
*International Journal of Damage Mechanics*, vol. 22, no. 6, pp. 878–904, 2013. View at Publisher · View at Google Scholar · View at Scopus - S. Cen, X. R. Fu, and M. J. Zhou, “8- and 12-node plane hybrid stress-function elements immune to severely distorted mesh containing elements with concave shapes,”
*Computer Methods in Applied Mechanics and Engineering*, vol. 200, no. 29–32, pp. 2321–2336, 2011. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - S. Cen, G. Zhou, and X. R. Fu, “A shape-free 8-node plane element unsymmetric analytical trial function method,”
*International Journal for Numerical Methods in Engineering*, vol. 91, no. 2, pp. 158–185, 2012. View at Publisher · View at Google Scholar · View at Scopus - H. A. F. A. Santos, “Complementary-energy methods for geometrically non-linear structural models: an overview and recent developments in the analysis of frames,”
*Archives of Computational Methods in Engineering*, vol. 18, no. 4, pp. 405–440, 2011. View at Publisher · View at Google Scholar · View at Scopus - H. A. F. A. Santos and C. I. Almeida Paulo, “On a pure complementary energy principle and a force-based finite element formulation for non-linear elastic cables,”
*International Journal of Non-Linear Mechanics*, vol. 46, no. 2, pp. 395–406, 2011. View at Publisher · View at Google Scholar · View at Scopus - C. Y. Dong and G. L. Zhang, “Boundary element analysis of three dimensional nanoscale inhomogeneities,”
*International Journal of Solids and Structures*, vol. 50, no. 1, pp. 201–208, 2013. View at Publisher · View at Google Scholar · View at Scopus - C. Y. Dong and E. Pan, “Boundary element analysis of nanoinhomogeneities of arbitrary shapes with surface and interface effects,”
*Engineering Analysis with Boundary Elements*, vol. 35, no. 8, pp. 996–1002, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at Scopus - S. S. Chen, Q. H. Li, Y. H. Liu, and Z. Q. Xue, “A meshless local natural neighbour interpolation method for analysis of two-dimensional piezoelectric structures,”
*Engineering Analysis with Boundary Elements*, vol. 37, no. 2, pp. 273–279, 2013. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - S. S. Chen, Y. H. Liu, J. Li, and Z. Z. Cen, “Performance of the MLPG method for static shakedown analysis for bounded kinematic hardening structures,”
*European Journal of Mechanics A: Solids*, vol. 30, no. 2, pp. 183–194, 2011. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Scopus - 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. View at Publisher · View at Google Scholar · View at Scopus - Y. J. Peng, J. W. Pu, B. Peng, and L. J. Zhang, “Two-dimensional model of base force element method (BFEM) on complementary energy principle for geometrically nonlinear problems,”
*Finite Elements in Analysis and Design*, vol. 75, pp. 78–84, 2013. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus - Y. J. Peng, L. J. Zhang, J. W. Pu, and Q. Guo, “A two-dimensional base force element method using concave polygonal mesh,”
*Engineering Analysis with Boundary Elements*, vol. 42, pp. 45–50, 2014. View at Publisher · View at Google Scholar · View at MathSciNet - Y. J. Peng, N. N. Zong, L. J. Zhang, and J. W. Pu, “Application of 2D base force element method with complementary energy principle for arbitrary meshes,”
*Engineering Computations*, vol. 31, no. 4, pp. 1–15, 2014. View at Google Scholar - 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. View at Publisher · View at Google Scholar · View at Scopus - Y. J. Peng, Y. H. Liu, J. W. Pu, and L. J. Zhang, “Application of base force element method to mesomechanics analysis for recycled aggregate concrete,”
*Mathematical Problems in Engineering*, vol. 2013, Article ID 292801, 8 pages, 2013. View at Publisher · View at Google Scholar - J. C. Walraven and H. W. Reinhardt, “Theory and experiments on the mechanical behavior of cracks in plain and reinforced concrete subjected to shear loading,”
*Heron*, vol. 26, no. 1, pp. 1–68, 1981. View at Google Scholar · View at Scopus - Z. L. Wang, F. Lin, and X. L. Gu, “Numerical simulation of failure process of concrete under compression based on mesoscopic discrete element model,”
*Tsinghua Science and Technology*, vol. 13, supplement 1, pp. 19–25, 2008. View at Publisher · View at Google Scholar · View at Scopus