Research Article  Open Access
Junqing Liu, Fan Zuo, Chao Liu, "A TransScale Young’ Modulus Calculation Model of ITZ Based on Void Shape Randomness and Calcium Hydroxide Enrichment", Advances in Materials Science and Engineering, vol. 2020, Article ID 9430875, 12 pages, 2020. https://doi.org/10.1155/2020/9430875
A TransScale Young’ Modulus Calculation Model of ITZ Based on Void Shape Randomness and Calcium Hydroxide Enrichment
Abstract
The randomness of void shape and enrichment of calcium hydroxide are significant in interfacial transition zone (ITZ) of concrete; however, current theoretical models of ITZ do not include these features. In this article, ITZ was regarded as a threephase composite material, and the pore morphological parameters were defined according to the characteristics of microscopic pores, and the corresponding random distribution function was constructed. The calcium hydroxide enrichment factor was introduced, and a crossscale ITZ Young’s modulus calculation model was established in combination with the Mori–Tanaka method. The reliability of the proposed model in this paper was verified through comparison to experimental results in a reference.
1. Introduction
As the weak area between aggregate and bulk paste in concrete, the interfacial transition zone (ITZ) influences mechanical properties of concrete, including Young’s modulus, strength, and cracking behavior significantly [1–7]. Concrete was treated as a twophase composite composed of aggregate and bulk paste in early studies, but study on limestone aggregate concrete conducted by Nilsen and Monterio [8] revealed that the tested Young’s modulus values were lower than Hashin–Shtrikman lower limit under twophase assumption. Simeonov and Ahamd [5] observed increasing deviation between test results and Hashin–Shtrikman limit under twophase assumption with the increasing of aggregate volume fraction.
ITZ was began to be taken into consideration [3] to solve the problem of large difference between the modelcalculated results and experimental results. Analysis of Chun’an et al. [9] showed that the properties of transition zone between inclusions and the matrix influence the stress field inside the composite, which then influences damage location and failure mode consequently. Studies of Liu et al. [10, 11] revealed the impact of ITZ on longterm deformation and bond behavior of recycled concrete. Threephase models considering ITZ were established by Lee and Park [1], Zhao and Chen [12], and Nadeau [13] et al., and reliable results were obtained. The failure mechanism of concrete under high strain rate with a threephase model was studied by Zhou and Hao [7], and this study showed that it influences the failure mode, failure time, and tensile strength remarkably. Berger [14] held that the wall effect led to higher porosity of cement paste on aggregate surface, which created a suitable circumstance for migration of water and ions and resulted in enrichment of calcium hydroxide and ettringite on aggregate surface. Simeonov and Ahamd [5] et al. pointed out the unsteadiness of ITZ properties along direction perpendicular to aggregate surface.
Experimental and theoretical methods were used to study many aspects of ITZ. The ITZ nanoindentation tests of gravel concrete and limestone concrete by Mondal [15] revealed that Young’s modulus of ITZ has a gradient of change in the direction perpendicular to the surface of the aggregate, and the test results at different measuring points have large dispersion. Other experimental results on ITZ in selfcompacting concrete [16], rockfilled concrete [17], nanosilica modified concrete [18], and sugarcane industrial ash modified concrete [19] also showed similar characteristics. Xray diffraction (XRD) was used by Hagis [20] to study the microstructure of cement hydration products, and the result showed that the microstructure and pore characteristics of calcium hydroxide, ettringite, and hydrated calcium silicate gels are significantly different. Mercury intrusion test was used by Liu et al. [21] to study the influence of porosity on water and ion migration in ITZ. Scanning electron microscopy (SEM) was used by Barnes et al. [22] to study the evolution of the microstructure of the glasscement mortar interface transition zone. The result showed that calcium hydroxide enrichment contributes greatly to the formation of ITZ high porosity. Sun et al. investigated the effects of aggregate surface by coating slag and silica fume on the ITZ microstructure [23] and studied the transport behavior of sulfate ions in concretecontaining ITZ [24]. Their studies showed that the ettringite formation inside the concrete pores influences the elastic deformation behavior of concrete significantly [25].
The most common theoretical method for calculating ITZ Young’s modulus is the discount coefficient method: a discount coefficient is imposed on Young’s modulus of cement paste, and the discounted result is used as the corresponding parameter of ITZ [7, 12, 13]. Luts and Zimmerman [26] believed that ITZ Young’s modulus is related to aggregate size, distance from aggregate surface, and modulus of bulk paste. They assumed that ITZ Young’s modulus can be calculated by a powerlaw function. Nadeau [27] proposed a theoretical calculation method of ITZ Young’s modulus, which comprehensively considered the influence of watercement ratio, aggregate volume fraction, cement specific gravity, ITZ thickness, and other factors. However, so far, there is still no ITZ Young’s modulus calculation model that can essentially consider the void shape randomness and calcium hydroxide enrichment caused by the characteristics of ITZ components and reflect the variation along direction perpendicular to aggregate surface and the discreteness among results of different measured points.
A transscale calculation model of ITZ Young’s modulus is proposed in this paper. In this model, ITZ is adopted to threephase composite, void shape parameter and calcium hydroxide enrichment parameter are defined, and then ITZ Young’s modulus could be calculated with the Mori–Tanaka method. The calculation model comprehensively considers the randomness of void shape and the degree of calcium hydroxide enrichment caused by the characteristics of ITZ components. It has the advantages of covering the upper and lower limits of Young’s modulus, reflecting the gradient of Young’s modulus change and small calculation. By comparison with the test results in a reference, the validity of the calculation model is verified. This model could calculate ITZ Young’s modulus reliably and quickly and provides a theoretical basis for the calculation of concrete mechanical properties.
2. ThreePhase Composite Model of ITZ
ITZ is a thin layer caused by the wall effect occurred between the aggregate and the bulk paste [14]. The resistance of aggregate on water permeation and cement accumulation lead to higher watercement ratio and porosity compared with cement paste; meanwhile, the enrichment of calcium hydroxide crystal is rather notable. Thus, ITZ could be considered as a layer of nearfield cement paste with high watercement ratio, high porosity, and enriched calcium hydroxide crystal, and the watercement ratio, porosity, and enrichment of calcium hydroxide reduce with the distance to the aggregate increasing.
The major components of ITZ remain consistent with the cement paste [15], which include hydrated calcium silicate gel (CSH, which contains lowdensity CSH and highdensity CSH, contracts to LDCSH and HDCSH, respectively), calcium hydroxide (CH), and ettringite. However, the volume fractions of components in ITZ are different with those in cement paste, considering the enrichment of CH. The microcrystal shape and pore features of ITZ components differ largely through electron microscopy [20, 28, 29], as is shown in Figure 1.
(a)
(b)
(c)
(d)
The macromechanical properties of cement paste in far or nearfield are closely related to micromechanical properties and volume fraction of components, while those of components are closely related to corresponding microstructures and pore features. The micro and macrofeatures of components should be taken into consideration in the calculation of ITZ mechanical properties, and a transscale model could be established.
Cement paste was treated as a fourphase composite in previous studies [30]: micropores, LDCSH, HDCSH, and ettringite. However, it is necessary to take the differences between ITZ and cement paste into account: pores stem from the special crystal structures of CH and ettringite rather than air holes; the volume fraction of micropores is not equal to the absolute porosity of bulk paste, which leads to the difference on mechanical properties and volume fractions between ITZ components and bulk paste components.
In addition, phases of ITZ grow and stagger throughout the hydration of cement. Influence of local pores leads to local randomness of property in each phase, following the local randomness of phases leads to general randomness of composite. More attention on local randomness led by pores is necessary in transscale modeling.
This paper presents a transscale ITZ model used for Young’s modulus calculation, which takes the randomness of pores and phases into account. This model adopted a threephase model, which consists of LDCSH, HDCSH, and CH, voids with shape randomness distributed in three phases. The assumption of continuously uniform media was applied into cement paste, and the relative porosity of ITZ to cement paste was set as the porosity of phases; methods of composite materials were used in mechanical property calculation of phases.
The following 3 steps were needed in transscale modeling:(1)In the microscale level 1st step: it was assumed that all the pores distribute in the phases continuously. Each phase was regarded as an accumulation of representative volume element (RVE) with an independent void in the center, and the porosity of each RVE was consistent with others. The shape of voids corresponded with some distribution, the shape parameter was defined, and corresponding distribution functions were constructed so that the features of void shape and distribution were described. The method of composite material was used to gain the mechanical properties of each RVE. 2nd step: each RVE was treated as a part of a certain phase, and thus, each phase could be regarded as a combination of nRVEs, that is, an nphase composite. The mechanical properties of each phase were obtained through the composite material method.(2)In the macroscale level 3rd step: ITZ was considered as a 3phase composite consisted of LDCSH, HDCSH, and CH, and the macro Young’s modulus of ITZ was obtained.
The modeling of the transscale model is shown in Figure 2.
(a)
(b)
3. TransScale Model of Young’s Modulus in ITZ
3.1. Equivalence of Microvoid Shape and the Shape Parameter
An accurate description on void features of phases in ITZ is necessary for calculation. Voids could be treated as special inclusions with zero modulus. The inclusion with a shape of rotating ellipsoid, in which 2 in 3 axes are of equal length, is a kind of inclusions with excellent universal adaptability. Various kinds of inclusions, such as spherical, ellipsoidal, pennyshaped, and cylindrical inclusions [31], are the generalization of rotating ellipsoidal inclusion. Considering the complexity of void shape in ITZ, this paper assumed that all the voids in ITZ phases are in the shape of rotating ellipsoid.
The phases in cement paste appear significantly different according to the SEM images: the crystal polymer of CSH is dense and uniform, and that of CH appears like a peony, while that of ettringite is needlelike, which is easy to form an extremely loose structure. Based on the shape and shape randomness of voids in phases, it is a reasonable assumption that the RVEs corresponding to CSH contain mainly spherical voids, CH contains mainly rotating ellipsoidal voids, and ettringite are treated as a pure void without any matrix. Considering the staggering of all the phases, it is acceptable to put the “ettringite voids” into other RVEs as part of porosity instead of an independent phase, as is shown in Figure 3.
A typical RVE was set to be a cube with 2L in side length, containing a rotating ellipsoidal void with in 3 axis length, respectively, where . The void shape parameter was defined as
The rotating ellipsoidal void could be classified to 3 types according to the value of α: flat rotating ellipsoidal void, spherical void, and prolate rotating ellipsoidal void, as is shown in Figure 4. Eshelby tensors [31] of the voided RVE corresponding to different α are shown in Appendix A.
Based on the fact that and , , it follows thatand by reforming equation (2), we can get
The domain of could be concluded as
Denoteand then the domain of was reformed to
For the purpose of a more tractable parameter, denote , and then equation (6a) could be reformed further as
The void is spherical when , flat rotating ellipsoidal when , and prolate rotating ellipsoidal when . and are equivalent in the following analysis.
3.2. Distribution Functions of Shape Parameter
The differences among microstructures of each phase in cement paste led to the differences among corresponding void shapes. For the voids in the RVE corresponding to CSH, the most probable values of shape parameters were , while those to CH were , or . In accordance with the most probable values, the possible probability density functions of and are shown in Figure 5. Five kinds of distribution functions were constructed consequently to describe the distribution of and .
3.2.1. Distribution of Shape Parameter for
In this case, , ; that is, , . A reasonable assumption is that follows the lognormal distribution [32]:and the distribution function is
3.2.2. Distribution of Shape Parameter for
In this case, . The lognormalPareto combined distribution was constructed [33] to describe the distribution of shape parameters, denoted as , and the distribution function iswhere
In the extreme case, , and . follows Pareto distribution, denoted as . The distribution function is
3.2.3. Distribution of Shape Parameter for
In this case, . The domain of is completely conjugated with that in Section 3.2.2. Therefore, transforming to a conjugated parameter and then the conjugated distribution function could be reconstructed. The method is as follows.
Denote as the conjugated parameter of , as the conjugated distribution function of . Obviously,
The monotonicity of and is in contrary, then and should be monotone decreasing consequently, and numerical range of which should be . Guaranteeing the monotonicity and numerical range, could be constructed as
The reconstructed distribution is named conjugated Paretolognormal combined distribution, denoted as . The distribution function is
The calculation of and is the same with equation (11).
In the extreme case, and . follows the conjugated Pareto distribution, denoted as , and the distribution function is
As the distribution functions are obtained, could be obtained with a known , and other parameters in distribution functions could be calculated with a corresponding quantile.
3.3. Equivalent Modulus of RVEs Corresponding to Phases
The porosity of the ith phase was denoted as , the Poisson ratio as , and Young’s modulus as . This phase is composed of nRVEs, each RVE contains a void with a shape parameter of , and the porosities of all the RVEs are the same . Massive values could be gain through the Monte Carlo sampling method, with a certain distribution function correspondingly.
The average mechanical properties of RVEs are worthy to be adopted in consideration of the randomness of void axis directions in the RVE accumulation. The average volume modulus and the average shear modulus of an RVE with an valued shape parameter and random axis direction could be obtained by Mori–Tanaka method [34]:respectively, where and are decided by the Eshelby tensor [34].
The macroporosity, Poisson ratio, and Young’s modulus of each phase would be determined in the next part then the macro Young’s modulus of ITZ, consequently.
3.4. Related Parameters of Cement Paste
As the distance between the cement paste and the aggregate surface increases, the influence of the side wall effect gradually weakens, the local watercement ratio decreases, the porosity decreases, the degree of calcium hydroxide enrichment weakens, and the porosity and volume fraction of the constituent phases decrease. The porosity and volume fraction of the constituent phases also change.
ITZ is a layer of highporosity cement paste. Taking the porosity of the cement paste and the compositional phase volume fraction as a reference, the local porosity of ITZ was determined based on the local watercement ratio, and the relative porosity of ITZ was obtained. The degree of calcium hydroxide enrichment of ITZ was determined, based on which the volume fraction of the cement mortar composition phase was adjusted, and the volume fractions of the ITZ composition phases were obtained. After ITZ’s relative porosity and composition phase volume fraction were obtained, the ITZ modulus calculation can be performed.
The volume fractions of components in cement paste were denoted as and and CH as . The relative volume fractions among LDCSH, HDCSH, and CH are almost unchanged under the watercement ratio of 0.5 and 0.65, respectively [15]. Without test result under higher watercement ratio, it is reasonable to assume that, for the cement paste with a watercement ratio , the relative volume fractions of the 3 phases keep as the average value of that under 0.5 and 0.65 watercement ratio, as is shown in Table 1.

Young’s modulus of phases in cement paste is shown in Table 2, in which that of micropores(MC) was set to be zero in this paper. The Poisson ratios of all the phases were set to be 0.2.
 
The absolute porosity of cement paste, φ bulkabs, was set to be 30.7% when _{0} = 0.5 [35].
It is worth noting that Young’s modulus of phases in cement paste provided by [15] is affected by the micropores in the phases themselves, while the micropores are contained in the absolute porosity of cement paste. There would be a higher estimation on actual porosity if the test results in [13] are adopted directly. To avoid this problem, it is supposed to treat the volume fraction of micropores, , as the relative porosity of cement paste, . Then, the porosity of phases could be calculated as follows:in which and when [35], and thus, .
3.5. Macro Young’s Modulus of ITZ considering the Enrichment of CH
3.5.1. Local WaterCement Ratio and Relative Porosity of ITZ
As is discussed in 3.4, in accordance with the influence of bulk paste porosity on tested modulus of cement paste phases, it is necessary to adopt the relative porosity of ITZ to bulk paste in calculation of ITZ modulus, to avoid an over estimation on porosity.
The local watercement ratio, , could be calculated as [27]wherein which is the thickness of ITZ, is the radius of aggregate, is the volume fraction of aggregate, and is the watercement ratio of bulk paste, as is shown in Figure 6.
The relation between the distance from aggregate surface and porosity is [35]where is the absolute porosity when . , and [35] when .
With a determined , the relative porosity, , could be obtained as follows:
3.5.2. Enrichment Parameter of CH
The enrichment degree of CH is proportional to the local relative porosity. An enrichment parameter, , was introduced to estimate the enrichment degree of CH as follows:
The form of equation (23) depends on test results, which were determined in the 4th part of this paper. In accordance with an assumption of a constant volume fraction ratio between LDCSH and HDCSH, the volume fraction ratio of ITZ with consideration of CH enrichment could be obtained as follows:
3.5.3. Macro Young’s Modulus of ITZ
The volume and shear modulus of isotropic multiphase composite could be obtained by the Mori–Tanaka method [36]:where
The subscript 0 represents the bulk matrix, r is the rth phase, and is the volume fraction ratio of rth phase.
4. Result and Verification
Nanoindentation tests on gravelcement ITZ [13] were chosen for verification. Samples were made using Type I Portland cement, river sand, and gravel. Watercement ratio was fixed at 0.5, river sand in the size of 1.18 mm to 2.36 mm, and gravels with an approximate size of 5 mm were used. All the samples were designed to be 1 inch × 1 inch × 0.5 inch, with a single gravel in the center, and cured under water at 25°C.
Seven groups of test results located at 2.5 μm to 27.5 μm from the aggregate surface were provided, 2 of which were used for parameter fitting, and the thickness of ITZ, δ, was determined to be 20 μm. With these data, according to least square principle, the form of enrichment parameter was determined as follows:and fitted parameters in which are shown in Table 3.

Considering the randomness of void shape, the calculation results were not equal to test results directly, so it is reasonable to adopt all the test results for verification. The calculated local watercement ratio and local porosity and are shown in Table 4, and the parameters in distribution are shown in Appendix B.

Assumed the void shape of CH follows LN, P, con_P, PLN, and con_PLN distribution, respectively, 5 curves of ITZ Young’s modulus along the distance from aggregate surface were obtained, and comparisons with test results were executed, as is shown in Figures 7 and 8.
(a)
(b)
Among the tested points within ITZ (0 μm < d < 20 μm), the nearer located to the surface of aggregate, the more significant discreteness of the test results, and the lager difference among all the calculated results under different distributions. Based on the analysis in Part 3, the local porosity increases with the distance from the surface of aggregate decreases. The calculated results showed that high porosity has a greater influence on Young’s modulus of ITZ than low porosity, which was proved by the discreteness of test results.
Among the tested points within ITZ (0 μm < d < 20 μm), the further located to the surface of aggregate, the lower porosity of ITZ, and the smaller difference among all the calculated results under different distributions. This phenomenon showed that the low porosity has a lesser influence on Young’s modulus of ITZ, and which keeps consistent with the discreteness of test results.
One possible reason for the large discreteness of test results within bulk paste (d > 20 μm) is the thickness of ITZ was set 20 μm, but the actual thickness of ITZ among the tested points varied. For the 2nd and 4th group of test results, the peak value of Young’s modulus appeared at d = 22.5 μm, which indicated a lager thickness of ITZ than 22.5 μm; similarly, the 6th group of test results indicated 27.5 μm. Another possible reason is the randomness of Young’s modulus in bulk paste, which results in the difference of corresponding ITZ Young’s modulus. The deviation between test and calculated results would significantly decrease when ignoring the 2nd, 4th, and 6th group.
The calculated results were suitable for the estimation on upper limit of Young’s modulus when the void in CH was assumed to follow P or PLN distribution, the lower limit when con_P or con_PLN distribution, and the average value when LN distribution. The estimation on upper and lower limit is conservative relatively, while that on average value showed small deviation. Considering the discreteness of test results themselves, there is no sense to analyse the discreteness of calculated results. Deviations about all the estimation are shown in Tables 5–7.



5. Conclusion
A macro Young’s modulus calculation model of ITZ was proposed in this paper. With the introduction of void shape parameter and calcium hydroxide enrichment parameter, and the combination of local watercement ratio, aggregate volume fraction, and local porosity, Young’s modulus of ITZ could be obtained.
This model, with an advantage of small computation, could cover the upper and lower limit and average value of Young’s modulus of ITZ and reflect the gradient variation along direction perpendicular to aggregate surface. The availability of the model when = 0.5 was verified through comparison with test results in a reference paper.
For the sake of limited research scope of gravelcement paste transition zone and limited test data for verification, the proposed model was mainly adopted to ITZ in natural aggregate concrete at present. In addition, this suggested model could provide both theoretical support for meticulous modeling and failure analysis on concrete, and theoretical tool for improvement of ITZ properties in concrete strength research.
The suggested formula of calcium hydroxide enrichment parameter was empirical, while the coupling effect of water enrichment and ion migration was the basic cause of calcium hydroxide enrichment. This coupling effect would be the theme of further study.
Appendix
A. Eshelby Tensor Corresponding to Different Void Inclusion Shape
Remark : the expressions were derived from [31], unwritten components were calculated according to Voigt symmetry, and those inconsistent with Voigt symmetry were 0 (Tables 8–10).



B. Distribution Parameter Values Corresponding to Different Porosity (Tables 11–15)





Data Availability
Previously reported ITZ Young’s modulus data were used to support this study and are available at https://www.mendeley.com/catalogue/65b2e568a37e3f4c8fbc677ba0dbc068/. These prior studies (and datasets) are cited at relevant places within the text as references [15].
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The authors would like to acknowledge the financial support provided to this study by the National Natural Science Foundation of China (grant no. 51878546), the Innovative Talent Promotion Plan of Shaanxi Province (grant no. 2018KJXX056), the Key Research and Development Projects of Shaanxi Province (grant no. 2018ZDCXLSF030302), the Science and Technology Innovation Base of Shaanxi Province (grant no. 2017KTPT19), and the Outstanding Youth Science Foundation of Shaanxi Province (grant no.2020JC46).
References
 K. M. Lee and J. H. Park, “A numerical model for elastic modulus of concrete considering interfacial transition zone,” Cement and Concrete Research, vol. 38, no. 3, pp. 396–402, 2008. View at: Publisher Site  Google Scholar
 M. P. Lutz, P. J. M. Monteiro, and R. W. Zimmerman, “Inhomogeneous interfacial transition zone model for the bulk modulus of mortar,” Cement and Concrete Research, vol. 27, no. 7, pp. 1113–1122, 1997. View at: Publisher Site  Google Scholar
 S. Mindess, “Bonding in cementitious composites: how important is it?” MRS Online Proceedings Library Archive, vol. 114, 1987. View at: Publisher Site  Google Scholar
 K. L. Scrivener, A. K. Crumbie, and P. Laugesen, “The interfacial transition zone (ITZ) between cement paste and aggregate in concrete,” Interface Science, vol. 12, no. 4, pp. 411–421, 2004. View at: Publisher Site  Google Scholar
 P. Simeonov and S. Ahmad, “Effect of transition zone on the elastic behavior of cementbased composites,” Cement and Concrete Research, vol. 25, no. 1, pp. 165–176, 1995. View at: Publisher Site  Google Scholar
 K. Z. J. Wu, “The influence of the matrixaggregate bond on the strength and brittleness of concrete,” MRS Proceedings, vol. 114, 1987. View at: Publisher Site  Google Scholar
 X. Q. Zhou and H. Hao, “Mesoscale modelling of concrete tensile failure mechanism at high strain rates,” Computers & Structures, vol. 86, no. 2122, pp. 2013–2026, 2008. View at: Publisher Site  Google Scholar
 A. U. Nilsen and P. J. M. Monteiro, “Concrete: a three phase material,” Cement & Concrete Research, vol. 23, no. 1, pp. 147–151, 1993. View at: Publisher Site  Google Scholar
 T. Chun’an, F. Yufang, and Z. ancheng, “Numerical approach to effect of interface properties on failure modes in particle filled composite,” Acta Materiae Compositae Sinica, vol. 16, pp. 112–120, 1999, in Chinese. View at: Google Scholar
 C. Liu, Z. Fan, X. Chen, C. Zhu, H. Wang, and G. Bai, “Experimental study on bond behavior between section steel and RAC under full replacement ratio,” KSCE Journal of Civil Engineering, vol. 23, no. 3, pp. 1159–1170, 2019. View at: Publisher Site  Google Scholar
 C. Liu, Z. Lv, C. Zhu, G. Bai, and Y. Zhang, “Study on calculation method of long term deformation of RAC beam based on creep adjustment coefficient,” KSCE Journal of Civil Engineering, vol. 23, no. 1, pp. 260–267, 2019. View at: Publisher Site  Google Scholar
 X. H. Zhao and W. F. Chen, “Effective elastic moduli of concrete with interface layer,” Computers & Structures, vol. 66, no. 23, pp. 275–288, 1998. View at: Publisher Site  Google Scholar
 J. C. Nadeau, “A multiscale model for effective moduli of concrete incorporating ITZ watercement ratio gradients, aggregate size distributions, and entrapped voids,” Cement and Concrete Research, vol. 33, no. 1, pp. 103–113, 2003. View at: Publisher Site  Google Scholar
 R. L. Berger, D. S. Cahn, and J. D. Mcgregor, “Calcium hydroxide as a binder in portland cement paste,” Journal of the American Ceramic Society, vol. 53, 1970. View at: Publisher Site  Google Scholar
 P. Mondal, “Nanomechanical properties of cementitious materials,” Northwestern University, Evanston, IL, USA, 2008, Ph.D. thesis. View at: Google Scholar
 W. Zhu, “Microstructure and properties of interfacial transition zone in SCC,” in Proceedings of the SCC’2005China1st International Symposium on Design, Performance and Use of SelfConsolidating Concrete, Changsha, China, 2005. View at: Google Scholar
 Y. Xie, D. J. Corr, F. Jin, H. Zhou, and S. P. Shah, “Experimental study of the interfacial transition zone (ITZ) of model rockfilled concrete (RFC),” Cement & Concrete Composites, vol. 55, pp. 223–231, 2015. View at: Publisher Site  Google Scholar
 J. Xu, B. Wang, and J. Zuo, “Modification effects of nanosilica on the interfacial transition zone in concrete: a multiscale approach,” Cement & Concrete Composites, vol. 81, pp. 1–10, 2017. View at: Publisher Site  Google Scholar
 J. A. Rossignolo, M. S. Rodrigues, M. Frias, S. F. Santos, and H. S. Junior, “Improved interfacial transition zone between aggregatecementitious matrix by addition sugarcane industrial ash,” Cement & Concrete Composites, vol. 80, pp. 157–167, 2017. View at: Publisher Site  Google Scholar
 C. W. Hargis, Advances in Sustainable Cements, 2013, https://escholarship.org/uc/item/4n67v8d6.
 C. Liu, H. Liu, C. Zhu, and G. Bai, “On the mechanism of internal temperature and humidity response of recycled aggregate concrete based on the recycled aggregate porous interface,” Cement and Concrete Composites, vol. 103, pp. 22–35, 2019. View at: Publisher Site  Google Scholar
 B. D. Barnes, S. Diamond, and W. L. Dolch, “The contact zone between portland cement paste and glass “aggregate” surfaces,” Cement and Concrete Research, vol. 8, no. 2, pp. 233–243, 1978. View at: Publisher Site  Google Scholar
 D. Sun, H. Shi, K. Wu, S. Miramini, B. Li, and L. Zhang, “Influence of aggregate surface treatment on corrosion resistance of cement composite under chloride attack,” Construction and Building Materials, vol. 248, Article ID 118636, 2020. View at: Publisher Site  Google Scholar
 D. Sun, K. Wu, H. Shi, L. Zhang, and L. Zhang, “Effect of interfacial transition zone on the transport of sulfate ions in concrete,” Construction and Building Materials, vol. 192, pp. 28–37, 2018. View at: Publisher Site  Google Scholar
 D. Sun, K. Wu, H. Shi, S. Miramini, and L. Zhang, “Deformation behaviour of concrete materials under the sulfate attack,” Construction and Building Materials, vol. 210, pp. 232–241, 2019. View at: Publisher Site  Google Scholar
 M. P. Lutz and R. W. Zimmerman, “Effect of the interphase zone on the bulk modulus of a particulate composite,” Journal of Applied Mechanics, vol. 63, no. 4, pp. 855–861, 1996. View at: Publisher Site  Google Scholar
 J. C. Nadeau, “Watercement ratio gradients in mortars and corresponding effective elastic properties,” Cement and Concrete Research, vol. 32, no. 3, pp. 481–490, 2002. View at: Publisher Site  Google Scholar
 J. Xiao, W. Li, Z. Sun, D. A. Lange, and S. P. Shah, “Properties of interfacial transition zones in recycled aggregate concrete tested by nanoindentation,” Cement and Concrete Composites, vol. 37, pp. 276–292, 2013. View at: Publisher Site  Google Scholar
 A. C. Gray, Crystal Morphology and Surface Reactivity Studies of Calcium Hydroxide, Brunel University, Uxbridge, UK, 1990.
 Z. Weiling, S. Wei, C. Cuicui, and M. Changwen, “Characterization for micro mechanical properties of cementitious materials by nanoindentation technique,” Journal of Southeast University (Natural Science Edition), vol. 41, pp. 160–165, 2011, in Chinese. View at: Google Scholar
 X. Jin, D. Lyu, X. Zhang, Q. Zhou, Q. Wang, and L. M. Keer, “Explicit analytical solutions for a complete set of the Eshelby tensors of an ellipsoidal inclusion,” Journal of Applied Mechanics, vol. 83, no. 12, 2016. View at: Publisher Site  Google Scholar
 D. G. Raheja and M. Allocco, Assurance Technologies Principles and Practices: A Product, Process, and System Safety Perspective, John Wiley & Sons, Hoboken, NJ, USA, 2006.
 D. P. M. Scollnik and P. M. David, “On composite lognormalPareto models,” Scandinavian Actuarial Journal, vol. 2007, no. 1, pp. 20–33, 2007. View at: Publisher Site  Google Scholar
 G. P. Tandon and G. J. Weng, “Average stress in the matrix and effective moduli of randomly oriented composites,” Composites Science and Technology, vol. 27, no. 2, pp. 111–132, 1986. View at: Publisher Site  Google Scholar
 Y. Gao, G. De Schutter, G. Ye, Z. Tan, and K. Wu, “The ITZ microstructure, thickness and porosity in blended cementitious composite: effects of curing age, water to binder ratio and aggregate content,” Composites Part B Engineering, vol. 60, pp. 1–13, 2014. View at: Publisher Site  Google Scholar
 G. J. Weng, “Some elastic properties of reinforced solids, with special reference to isotropic ones containing spherical inclusions,” International Journal of Engineering Science, vol. 22, no. 7, pp. 845–856, 1984. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Junqing Liu 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.