Research Article  Open Access
Application of Base Force Element Method to Mesomechanics Analysis for Concrete
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 threephase 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 twodimensional 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 twodimensional 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 threedimensional 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 stressstrain 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 stressstrain 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 stressstrain 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.
(a) Specimen 1 with
(b) Specimen 2 with
(c) Specimen 3 with
The propagation process of cracks of the specimen with by uniaxial compression is shown in Figure 12.
(a) Specimen 1 with
(b) Specimen 2 with
(c) Specimen 3 with
The propagation process of cracks of the specimen with by uniaxial compression is shown in Figure 13.
(a) Specimen 1 with
(b) Specimen 2 with
(c) Specimen 3 with
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 threephase 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
 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 Site  Google Scholar
 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 Site  Google Scholar
 E. Schlangen and E. J. Garboczi, “Fracture simulations of concrete using lattice models: computational aspects,” Engineering Fracture Mechanics, vol. 57, no. 23, pp. 319–332, 1997. View at: Publisher Site  Google Scholar
 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
 Y. J. Peng, B. K. Li, and B. Liu, “Numerical simulation of mesolevel mechanical properties of roller compacted concrete,” Journal of Hydraulic Engineering, vol. 32, no. 6, pp. 19–22, 2001 (Chinese). View at: Google Scholar
 X. L. Du, L. Jin, and G. W. Ma, “Mesoelement 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 Site  Google Scholar
 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 Site  Google Scholar
 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 Site  Google Scholar
 X. L. Du and L. Jin, “Mesoscale 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 mesoscale 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 Site  Google Scholar
 S. Cen, X. R. Fu, and M. J. Zhou, “8 and 12node plane hybrid stressfunction 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 Site  Google Scholar  MathSciNet
 S. Cen, G. Zhou, and X. R. Fu, “A shapefree 8node 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 Site  Google Scholar
 H. A. F. A. Santos, “Complementaryenergy methods for geometrically nonlinear 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 Site  Google Scholar
 H. A. F. A. Santos and C. I. Almeida Paulo, “On a pure complementary energy principle and a forcebased finite element formulation for nonlinear elastic cables,” International Journal of NonLinear Mechanics, vol. 46, no. 2, pp. 395–406, 2011. View at: Publisher Site  Google Scholar
 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 Site  Google Scholar
 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 Site  Google Scholar  Zentralblatt MATH
 S. S. Chen, Q. H. Li, Y. H. Liu, and Z. Q. Xue, “A meshless local natural neighbour interpolation method for analysis of twodimensional piezoelectric structures,” Engineering Analysis with Boundary Elements, vol. 37, no. 2, pp. 273–279, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 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 Site  Google Scholar  MathSciNet
 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 Site  Google Scholar
 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 Site  Google Scholar
 Y. J. Peng, Z. L. Dong, B. Peng, and N. N. Zong, “The application of 2D base force element method (BFEM) to geometrically nonlinear analysis,” International Journal of NonLinear Mechanics, vol. 47, no. 3, pp. 153–161, 2012. View at: Publisher Site  Google Scholar
 Y. J. Peng, J. W. Pu, B. Peng, and L. J. Zhang, “Twodimensional 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 Site  Google Scholar  MathSciNet
 Y. J. Peng, L. J. Zhang, J. W. Pu, and Q. Guo, “A twodimensional base force element method using concave polygonal mesh,” Engineering Analysis with Boundary Elements, vol. 42, pp. 45–50, 2014. View at: Publisher Site  Google Scholar  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. 34, pp. 171–183, 2003. View at: Publisher Site  Google Scholar
 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 Site  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
 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 Site  Google Scholar
Copyright
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.