Advances in Condensed Matter Physics

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

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

## Study on the Macromesoscopic Yield Criterion of Concrete Material

^{1}State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Hohai University, Nanjing 210098, China^{2}Department of Engineering Mechanics, Hohai University, Nanjing 210098, China^{3}Communication and Electron College, Jiangxi Science & Technology Normal University, Nanchang 330013, China^{4}Department of Civil Engineering, Nanjing Institute of Technology, Nanjing 211167, China

Received 6 September 2013; Revised 23 November 2013; Accepted 23 November 2013

Academic Editor: Zaoyang Guo

Copyright © 2013 Xia XiaoZhou 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

In mesoscopy, the concrete material is a kind of composite consisting of mortar matrix, aggregate, and the interface between them. And the shape of aggregate is assumed to be spherical and the mortar matrix is supposed to satisfy the Drucker-Prager yield criterion. The energy density support function is introduced to reflect the yield surface of mortar matrix. In order to solve the nonderivability on the yield boundary, the function approximation series has been constructed to substitute for the energy density support function. Finally, based on the function series, a macromesoscopic yield criterion of concrete material has been derived by the nonlinear homogenized technique. Through the macromesoscopic yield criterion established, the influence law of the bonding status of interface and the fraction of aggregate on the macrofriction coefficient is explored.

#### 1. Introduction

In mesoscopy, the concrete can be taken as a multiphase composite material which is made of aggregate, mortar matrix, and the interface between them. And the macroscopic mechanical properties of concrete depend on the mechanical characteristic of the mesocomponents and their connection status. Compared with the well-rounded estimation theory of the effective modulus, the estimation of composite strength which is very sensitive to the mesostructures is still not mature enough. For all that, there are still many scholars developing this research. Huaxiang et al. [1] used the homogenization method of mesomechanics to estimate the ultimate loads applied and the macroyield criterion of composites is further obtained through the corresponding relationship between the macromechanical property and the meso structures. Barthelemy and Dormieux [2] studied the relationship of macro friction coefficient with the matrix friction coefficient and aggregates volume occupancy under the perfect interface case and the non-friction sliding interface case by using the meso-inclusion theory and obtained the macroconstitutive model which can reflect the mesoparameters. For the composite strength of cement matrix, by many components experiment, Maalej et al. [3] obtained the conclusion that the cement matrix occupy contributes the cohesive strength as proportional pattern and has almost no effect on the macrofriction coefficient and made explanation of mechanical mechanism for this experimental results. For the estimation of effective property of composite material considering the nonlinear characteristic of matrix, Willis [4, 5] studied the effective property estimation problem of composite material in which the nonlinear matrix includes the rigid inclusion and void inclusion by using Hashin-Shtrikman variational structure method [6, 7]. Ponte Castañeda [7–11] studied the estimation of effective property of composite material in which the matrix is nonlinear as power form and the bonds between matrix and inclusions are supposed perfect by using his proposed variational structure method. Zhang and Xia [12] studied the nonlinear composite properties of concrete and explored the influence of interface connect status on the concrete effective properties considering the nonlinear soft characteristic of cement matrix by variational structure method proposed by Ponte Castañeda.

In this paper, the aggregate component is supposed to never yield and its shape is assumed to be spherical. For the mortar matrix as friction type of material, the Drucker-Prager yield criterion is applied to describe its yield behavior. The macroscopic yield criterion can be established by introducing the potential function method and combining the estimation of the effective modulus. And the influence law of aggregate content and bond characteristic to the inner friction coefficient of concrete can be explored by the established macroscopic yield criterion. The potential function will use a kind of support function, the energy density support function.

The average Eshelby tensor which can reflect the connecting status of interface between the inclusion and the matrix including imperfect connecting status [13, 14] is introduced to explore the influence law of interface connecting status on macroyield characteristic of concrete.

The energy density support function [14] has been used to reflect the yield surface of the mortar matrix successfully. And through this function, the influence of the air-entrained water reduce agent on concrete inner friction coefficient was explored subsequently [15]. This paper will still adopt the energy density support function and its function approximation series to describe the yield surface of concrete, and the macromesoscopic yield criterion is expected to be established.

#### 2. The Construction of the Function Approximation Series

The concrete sample can be taken as a representative volume element which is a multiphase composite material made of cement mortar (matrix), aggregate, and the interface between them. For the mortar and interface phase, the Drucker-Prager yield criterion is applied as follows: where , , , and . Among them, , and are the cohesion and the friction angle of the material, respectively.

Energy density is defined as follows: where , , and . Since is the convex domain, the energy density support function can be taken to the supremum. The support function form can be written as . For Drucker-Prager criterion, the expression is as follows:

When applying a constant strain on the representative volume element boundary, the basic equation is as follows: where is the mortar domain, and is the aggregate domain, is the energy density support function of the matrix. Let , and consider the nonderivability at of function; a derivable convex function series should be constructed. And when parameter tends to be zero, the function series constructed should approach to function. For Drucker-Prager yield function form, it can be written as follows: The solution of (4) can be obtained as follows: As long as when the parameter tends to be zero, the function reaches positive infinity as presented in Figure 1, and then the function series also reaches positive infinity . So, the constructed derivative convex function series meet the requirement of approaching to function. Thus, the third subformula of (4) becomes as follows:

#### 3. Nonlinear Homogenization Technique

The stress-strain relationship can be derived as follows by the potential function (6): where , , and , which are the bulk modulus, two times of shear modulus, and the spherical prestressed, respectively. Since is negative and except for the stress status , , are positive and depend on and . In addition, if , the equal formula of can be obtained from the definition of the function . Since the strain fields both in the matrix domain and in the aggregates domain of the question (7) are nonuniform, , take on nonuniform distribution. Based on the classical homogenization method [16] which is only to deal with the heterogeneous material consisted of homogeneous phases, the real strain field is substituted by an equivalent alternative strain , and the non-uniform problem (7) can use the classical homogenization method to deal with as follows: where, Among them, . Then (7) becomes as follows:

In order to ensure that the effective strain is not rougher compared with true strain fields, the first-order form is adopted to estimate the bulk strain and the second-order form to estimate deviator strain invariants as follows: where is an averaged operator. The relationship between the average strain of inclusion and the applied strain at far field can be deduced by the following energy equivalent principle:

The partial derivatives of to and can be gained respectively, by formula (12) as follows: It can be further obtained as follows by the assumption of material isotropic characteristic: The Mori-Tanaka method [17] is adopted to estimate the average strain of the inclusions. This method can consider the interaction between inclusions by changing the far field strain to the average strain of matrix which is equivalent to the action of average strain of matrix applied to far field when each inclusion is put embedded into the infinity matrix. So, where is the strain localization tensor, and . Here, is the average Eshelby tensor. and are softness and stiffness tensor of the matrix, respectively. Because of the application of the average Eshelby tensor in the strain localization tensor, the inclusion average strain includes the influence of relative displacement generated by the imperfect interface. So the formula (15) can be further trimmed as follows: And subsequently, the bulk average strain of inclusion is obtained as The mechanical property of all the aggregate is assumed to be the same, and the stiffness ratios are defined as and . The average Eshelby tensor expressions can be given as follows by [14]: where , .

Among , (for spherical inclusion). The determination of the coefficients , , , , , , and are related with the interface parameters. Assuming that the interface is of elastic and isotropic characteristic, the interface constitutive relation equation can be decomposed into two parts presented as follows when the interface embedding effect and the couple effect of constitutive relationship between the normal direction and tangential direction is not considered: where , denotes the external normal direction of the interface, the displacement at interface, and the interface stress component. It can be seen from (18) that and denote the interface characterization behaviors. For example, is corresponding to the perfect interface, is corresponding to the slip in tangential direction and perfect in normal direction, and is corresponding to the free sliding status of interface. The linear interface condition is essentially equivalent to using a nonthickness linear spring layer to simulate the connection behavior of the imperfect interface. The details can be found in [14]. The estimation formula of the effective bulk modulus and shear modulus for the isotropic case [17, 18] can be easily derived by the formula (18). Consider where is the volume occupancy of the th phase aggregate.

#### 4. The Macroscopic Yield Criterion

Since the constitutive relation of the matrix and aggregates in (11) is linear, it is decomposed into two basic problems, one is that the representative volume element of concrete is subjected by and , respectively, the other is that the representative volume element of concrete is subjected by and , respectively.

Here, we use and to denote the macroscopic stress invariant. Among them, .

Obviously the macroscopic stress generated by subproblem one is ; For subproblem two, the macroscopic stress generated depended on the effective modulus of the representative volume element of concrete. Thus, the macrostress-strain relationship can be gained as follows after superposition: Here, and are given by formula (17) and and are given by formula (9b). Since , the following formula can be gained: It can be further gained as follows: Now, we will derive the relationship between and . It can be analyzed from expression (14) that the solution of expression and is the key of gaining the relationship between and . Since the matrix modulus involved in and is constant, the two partial derivative functions can be solved as follows from formula (14). Consider where Substituting the formulas (16) and (24) into the second formula of the formula (13) and further simplification, the relation between and can be gained as follows: Substituting the formulas (17), (20), and (27) into formula (21), the macrobulk stress and macrodeviator stress can be gained as follows: Substituting the formulas (17) and (9b) into the formulas (27) and (28) respectively, and further simplification, the macrobulk stress and macrodeviator stress become as follows: The macro-micro constitutive relation of the concrete can be obtained as the follows by eliminating of formula (30) and formula (31): The above formula is further abbreviated as follows: where Among them, , .

It is worth to explain that the macroscopic friction coefficient could be of singularity because of the possibility that its denominator is zero. So it has to regulate the following expression to avoid the occurrence of singularity: where is a small tolerant threshold.

#### 5. Results Analysis

The relationship between the macroinner friction coefficient and microstructure can be explored by the macro-microconstitutive relations (33). And the laws of macroscopic inner friction coefficient with aggregate occupancy, stiffness ratio of mortar matrix to aggregate, and softness of interface for different matrix inner friction coefficient are presented as Figures 2–6. It can be seen from Figures 2 and 3 that the interface connection status is a key factor obviously influencing the overall inner friction coefficient , and the effective inner friction coefficient even decreased with the increment of the imperfect interface occupancy. What is more, the effective inner friction coefficient increases more rapidly with the increment of aggregate occupancy for infinite aggregate stiffness case, as presented Figure 3 than that of common aggregate stiffness case, as shown Figure 2. When all of the interface connection statuses are imperfect, the overall inner friction coefficient decreases with the increment of aggregate occupancy, as presented Figures 2–4. And it can be seen from Figure 4 that the effective inner friction coefficient decreases with the increment of stiffness ratio of mortar matrix to aggregate. Compared with Figure 5 reflecting the influence of tangential softness of interfaces on the macroinner friction coefficient, the influence of the normal softness of interface on the macroinner friction coefficient is more violent, as presented Figure 6. All the changing trends presented in the Figures 2–6 agree well with the experience intuition.

#### 6. Summary

In this paper, the concrete material is taken as a composite material consisting of aggregate, mortar matrix, and interface between them in the mesoscopic. Based on the assumption that the shape of aggregate is spherical, the aggregate phase is never yield and the mortar matrix satisfies Drucker-Prager yield criterion, and the macromesoscopic yield criterion of concrete material is derived by introducing the energy density supporting function to reflect the yield surface of the mortar matrix and the nonlinear homogenization technique. An approximation function series can be constructed to solve the nonderivable problem at the yield boundary. Finally, the influence laws of the interface connection status, aggregate content, and the rigidity ratio of matrix to aggregate on macroscopic inner friction coefficient are explored by the macromesoscopic yield criterion derived.

The results show that the influence laws gained agree well with the real case. Compared with the [2], which can only consider the estimation of the macroinner friction coefficient under the nonfriction interface case, this paper can not only deal with that of arbitrary interface connection status but also is not limited by the assumption of the aggregate rigid infinity. So this theory proposed takes a big step forward, but there is still some unestimated or estimation failure problem because of the singularity.

#### Conflict of Interests

The authors declared that they have no conflict of interests to this work.

#### Acknowledgments

This work is partially supported by the National Natural Science Foundation of China (nos. 51179064, 11372099, 11132003, and 51108231) and the state Key Laboratory Open Foundation of Hydrology-Water Resources and Hydraulic Engineering at Hohai University (no. 2011490911). Their financial support is gratefully acknowledged.

#### References

- L. Huaxiang, L. Yinghua, F. Xiqiao, and C. Zhangzhi, “A micromechanical method for determining the macroscopic yield criterion of composites,”
*Acta Mechanical Solida Sinica*, vol. 23, no. 2, pp. 133–139, 2002. View at Google Scholar - J. F. Barthelemy and L. Dormieux, “A micromechanical approach to the strength criterion of Drucker-Prager materials reinforced by rigid inclusions,”
*International Journal for Numerical and Analytical Methods in Geomechanics*, vol. 28, pp. 565–582, 2004. View at Google Scholar - Y. Maalej, L. Dormieux, J. Canou, and J. C. Dupla, “Strength of a granular medium reinforced by cement grouting,”
*Comptes Rendus Mécanique*, vol. 335, no. 2, pp. 87–92, 2007. View at Publisher · View at Google Scholar · View at Scopus - J. R. Willis, “The overall behaviour of a nonlinear matrix reinforced by rigid inclusions or weakened by cavities,” in
*Random Media and Composites*, 1984. View at Google Scholar - J. R. Willis, “Variational and related methods for the overall properties of composites,” in
*Advances in Applied Mechanics*, C.-S. Yih, Ed., vol. 21, pp. 1–78, Academic Press, New York, NY, USA, 1981. View at Google Scholar - Z. Hashin and S. Shtrikman, “A variational approach to the theory of the elastic behaviour of polycrystals,”
*Journal of the Mechanics and Physics of Solids*, vol. 10, no. 4, pp. 343–352, 1962. View at Google Scholar · View at Scopus - Z. Hashin and S. Shtrikman, “A variational approach to the theory of the elastic behaviour of multiphase materials,”
*Journal of the Mechanics and Physics of Solids*, vol. 11, no. 2, pp. 127–140, 1963. View at Google Scholar · View at Scopus - P. P. Castañeda, “The effective mechanical properties of nonlinear isotropic composites,”
*Journal of the Mechanics and Physics of Solids*, vol. 39, no. 1, pp. 45–71, 1991. View at Google Scholar · View at Scopus - P. P. Castañeda, “Exact second-order estimates for the effective mechanical properties of nonlinear composite materials,”
*Journal of the Mechanics and Physics of Solids*, vol. 44, no. 6, pp. 827–862, 1996. View at Publisher · View at Google Scholar · View at Scopus - P. Ponte Castañeda and P. Suquet, “Nonlinear composites,”
*Advances in Applied Mechanics*, no. 34, pp. 171–302, 1998. View at Google Scholar - P. Ponte Castañeda, “Second-order homogenization estimates for nonlinear composites incorporating field fluctuations—I—theory,”
*Journal of the Mechanics and Physics of Solids*, vol. 50, no. 4, pp. 737–757, 2002. View at Google Scholar · View at Scopus - Q. Zhang and X. Xia, “A macro-meso constitutive law for concrete having imperfect interface and nonlinear matrix,”
*Science in China E*, vol. 52, no. 3, pp. 634–640, 2009. View at Publisher · View at Google Scholar · View at Scopus - Z. Zhong and S. A. Meguid, “On the imperfectly bonded spherical inclusion problem,”
*Journal of Applied Mechanics, Transactions ASME*, vol. 66, no. 4, pp. 839–844, 1999. View at Google Scholar · View at Scopus - X. Xiao-Zhou,
*The meso numerical simulation and the macro-meso mechanics research for concrete material [Doctor dissertation]*, Hohai University, 2007. - J. Qun, X. Xiao-zhou, C. Ai-jiu, and Z. Qing, “Investigation on the influence of the air-entrained water reduce agent to the mechanical properties of recycle concrete,”
*Journal of SanXia University*, vol. 31, no. 3, pp. 46–51, 2009. View at Google Scholar - R. J. Asaro and D. M. Barnett, “The non-uniform transformation strain problem for an anisotropic ellipsoidal inclusion,”
*Journal of the Mechanics and Physics of Solids*, vol. 23, no. 1, pp. 77–83, 1975. View at Google Scholar · View at Scopus - T. M. Tanaka, “Average stress in matrix and average elastic energy of materials with misfitting inclusions,”
*Acta Metallurgica*, vol. 21, pp. 571–574, 1973. View at Google Scholar - T. Mura,
*Micromechanics of Defects in Solids*, Martinus Nijhoff Publishers, 1982.