Tensile Fracture Property of Concrete Affected by Interfacial Transition Zone

As a weak link between aggregate and mortar in concrete, interfacial transition zone (ITZ) usually plays a key role in concrete fracture. To investigate the tensile fracture property of concrete affected by the mechanical properties of ITZ numerically, the geometrical models of heterogeneous concrete were established with the parameterization modeling. They include three phases, namely, mortar, ITZ, and randomly distributed aggregates with distinct sizes and orientations. The cracking behaviors of mortar and ITZ were characterized by the bilinear cohesive zone constitutive model. Based on the experiments, the mechanical properties of ITZ were mediated by changing the water–cement ratio of mortar, the aggregate surface roughness and the content of silica fume in interfacial agent. A series of numerical simulations were conducted on the concrete models in tension after the numerical modeling method was validated. The macroscopic tensile fracture properties of concrete were quantitatively connected with some microscopic variables, including the water–cement ratio of mortar, the aggregate surface roughness and the silica fume content in interfacial agent. It was found that the tensile fracture properties of concrete have negative linear correlations with the water–cement ratio of mortar, while the effects of the aggregate surface roughness and the silica fume content in interfacial agent are very complex. The tensile fracture mechanical properties of concrete have a bilinear relationship with the aggregate surface roughness and an approximate quadratic parabola relationship with the content of silica fume in the interfacial agent. This study is beneficial to improve the fracture resistance of concrete by some interface handling measures.


Introduction
As the most widely used engineering building material, concrete has many advantages, including low price, easy construction, mechanized mass production, high compressive strength, good durability and wide strength grade. Concrete consists of aggregate, mortar and their interfacial transition zone (ITZ). Compared with mortar matrix, the ITZ usually has larger porosity (Liao et al., 2004;Tschegg et al., 1995;Yang, 1998), higher watercement ratio and permeability (Scrivener et al., 2004), and lower strength, toughness (Wu et al., 2016) and elastic modulus (Bourdette et al., 1995), so it is a weak link between aggregate and mortar in concrete and cracks tend to occur and propagate inside it. Accordingly, the ITZ property generally has a great influence on the fracture property of concrete, and to a certain extent determines crack path selection, as well as adhesion or cohesion failure mode of concrete subjected to external loads (Chen & Qian, 2017;Ding & Chen, 2019;Hong et al., 2014;Zhang et al., 2018aZhang et al., , 2018b. In general, there are two internal mechanisms for the bonding between aggregate and mortar (Guinea et al., 2002;Liao et al., 2004). One is the physical bonding mechanism. The bonding force is mainly related to the mechanical interlocking between the aggregates and the mortar matrix, and depends on the mineral properties and surface texture characteristics of aggregates (Caliskan, 2003;Gu et al., 2013;Hong et al., 2014;Rao & Prasad, 2004;Shen et al., 2019), as well as their shape and size (Caliskan, 2003). The other is the chemical bonding mechanism. When the aggregates are chemically active to the mortar, a chemical bond would be formed between the two (Asbridge et al., 2001;Bui et al., 2005;Duan et al., 2013;Ji et al., 2020;Jiang, 1999;Kuroda et al., 2000;Rossignolo, 2009;Siddique, 2011). Thus, the properties of ITZ can be physically or chemically changed. For example, it is believed that there is a poor bonding between the cement paste and asphalt-coated aggregates. Brand and Roesler () compared the ITZ for mortar with reclaimed asphalt pavement (RAP) aggregates relative to dolomite aggregates through image analysis of backscattered electron micrographs. They found that the ITZ with RAP aggregates was larger and more porous with less calcium silicate hydrate and calcium hydroxide at the interface than the ITZ with dolomite aggregates.
Based on the abovementioned two mechanisms of interfacial interaction, the ITZ properties can be improved from the following three aspects: (a) selecting appropriate aggregate type and aggregate size combination, and improving aggregate surface roughness (Caliskan, 2003;Jiang, 1999;Rao & Prasad, 2004;Shen et al., 2019); (b) adding mineral additives or admixtures with good compactness and chemical activity into the mortar (Asbridge et al., 2001;Duan et al., 2013;Jiang, 1999;Rao & Prasad, 2004); and (c) treating the aggregate surface chemically to produce functional groups that can be coupled to the mortar. Among them, the use of mineral additives or admixtures is a very effective way to improve the ITZ properties. For example, fly ash and silica fume were, respectively, added to concrete as the substitutes for cement to improve the mechanical properties of concrete (Golewski, 2021a(Golewski, , 2021b(Golewski, , 2021cSzostak & Golewski, 2020). It was demonstrated by a variety of experimental means that these additives could significantly increase the strength and fracture toughness of concrete and reduce the micro-crack size in the ITZ of concrete. The effects of fly ash and silica fume on the fracture properties of high-performance concrete were also experimentally investigated (Smarzewski, 2019;Zhang & Li, 2012). In these studies, however, the additives were added to cement paste, so that they affect not only the interface properties but also the properties of the matrix material.
By improving the adhesive effect of ITZ, the mechanical properties of concrete can be enhanced effectively. Due to the complex internal structure of concrete, different shapes and sizes of aggregates, and their random distribution, there are many factors affecting the ITZ bonding effect. Accordingly, it is very difficult to evaluate quantitatively the contribution of the ITZ property improvement to the macroscopic mechanical properties of concrete based on experiments, excluding the mortar property effect.
Many efforts have been made to develop numerical micromechanical models capable of predicting the mechanical behavior of heterogeneous concrete. For example, a concurrent multiscale modeling approach was proposed by Rodrigues et al. (2021;Rodriguesa et al., 2020) to investigate crack propagation in conventional concrete and recycled aggregate concrete. Utilizing the "take-and-place" method, coarse aggregates with regular shapes were generated from a given grading curve and then placed into the mortar matrix randomly. The mesh fragmentation technique (Manzoli et al., 2016;Sánchez et al., 2014) was used to treat the cracking of mortar and ITZ by introducing interface solid finite elements equipped with a tensile damage model. Recently, a modeling procedure named "meshplacement-identification-assignment" was developed by Zheng et al. (2021) to establish the mesoscale concrete model composed of aggregate and mortar. The effects of aggregate shape and aggregate content on the tensile strength of concrete were discussed. In our previous work Yin et al., 2011Yin et al., , 2012Yin et al., , 2015, the parameterization modeling was also developed for heterogeneous asphalt mixture and concrete. In the parameterization modeling, polyhedral or polygonal aggregates with a certain volume content and gradation were generated through some given parameters, and then randomly packed in a specified area to establish a numerical model of concrete including aggregates, matrix and their ITZs. Accordingly, it is potential that the parameterization modeling can be combined with experiments to investigate quantitatively the effect of the ITZ property improvement to the mechanical properties of concrete.
Recently, the effects of water-cement ratio of mortar, aggregate surface roughness and the content of silica fume in interface agent on the tensile fracture properties of ITZ were studied through the three-point bending test (Ji et al., 2020) and the splitting test (Zou et al., 2021) on the granite-mortar hybrid specimens. Based on the fitted experimental curves, the ITZ properties can be mediated by changing the water-cement ratio of mortar, aggregate surface roughness and the content of silica fume in interface agent. With the help of the parameterization modeling, the tensile fracture behaviors of concrete are simulated to investigate the direct improvement of the tensile fracture property of concrete from the improved ITZ properties quantitatively in this paper. This is very beneficial for engineering applications of concrete.

Parameterization Modeling
Taking the creation of two-dimensional concrete model as an example, the main steps of parameterization modeling Yin et al., 2011Yin et al., , 2012Yin et al., , 2015 are introduced as follows. The first step is to determine the total area of all diameter levels of aggregates according to a given aggregate gradation. The aggregate size of 4.75 mm was used as the cutoff size distinguishing coarse and fine aggregates (Li & Metcalf, 2005). The coarse aggregates are dispersed in the concrete as the reinforcing phase, while the fine aggregates are combined with cement and water to form the mortar matrix. To improve the modeling efficiency, circular aggregates would be first generated and randomly packed into a prescribed region, and then inscribed into the corresponding regular octagons, as shown in Fig. 1. In view of the loss of aggregate area caused by the inscription, it is necessary to modify the total areas of the aggregates with different particle sizes by means of the area conversion factor, k= 2 √ 2/π. The second step is to generate and pack circular aggregates. The generated circular aggregates with different sizes are randomly packed into a region, whose shape and size is prescribed according to the real specimen. To avoid conflicts between packed aggregates, their sizes are reduced by a common shrinking factor before packed, and restored to their original sizes after packed.
The third step is to perform a perturbation process until a stable equilibrium is reached for the aggregate system. After all the aggregates are restored to their initial sizes, a number of conflicts would inevitably occur between aggregates and aggregates or the wall, namely, the packing region boundary. It is necessary to redistribute their locations under the great repulsive forces between them through the perturbation process, so as to form a reasonable aggregate distribution.
In the fourth step, the circular aggregates are inscribed to their corresponding regular octagonal ones, so that the effect of aggregate shape can be considered in numerical simulations.
More detailed descriptions on the parameterization modeling can be seen in Ref. .

Cohesive Constitutive Model
It is assumed that the material is composed of simple cells connected by some virtual bonding surfaces, and damage and fracture only occur on these virtual bonding surfaces (Zhang et al., 2016). When the stress on a certain bonding surface exceeds the strength of the material, it would open (Yin et al., , 2015. In general, there is a nonlinear relationship between the stress acting on the bonding surface and its relative displacement. The graphical area enclosed by the stress-displacement curve and the horizontal axis is called the fracture energy.
It is further assumed that the relationship between the stress σ and the displacement δ can be described by a bilinear model (Zhang et al., 2018a(Zhang et al., , 2018b, as shown in Fig. 2. It has two independent parameters: the bond strength t 0 and the critical displacement δ c . When the stress is lower than the bond strength, the material is linear elastic, but after the stress reaches the bond strength, the material enters the linear softening stage and local damage occurs on the bonding surface. Finally, when the displacement reaches the critical displacement, the bonding surfaces would completely separate and a macro-crack would start. This is the bilinear cohesive constitutive model. In this case, the fracture energy G = 1/2δ c t 0 . When simulating the fracture process of concrete, considering that aggregates usually have higher cracking resistance than mortar matrix and ITZ, it is reasonably assumed that a crack can only form and propagate in the mortar matrix and ITZ. To save computing resources, the cohesive elements with zero thickness were inserted only in the mortar matrix and ITZ in finite-element modeling. More details can be found in Ref. (Yin et al., , 2015.

Three-Point Bending Fracture Experiment
To verify the feasibility of the above numerical simulation method, the three-point bending fracture experiments were carried out on six concrete specimens in the same batch. The experimental raw materials included P.I.42.5 type cement, crushed stone as coarse aggregate, medium sand as fine aggregate, and ordinary tap water, as shown in Fig. 3. After the crushed stones were screened by the sieving shaker, the 5-grade aggregates were obtained, as shown in Table 1. The concrete mix ratio is shown in Table 2. The mass ratio between water, cement and medium sand is 1:2:4.37.
The concrete specimen size is 400 mm × 100 mm × 100 mm. After poured in the mold, the concrete was fully tamped in the vibrostand. The specimen was released from the mold after being placed in the room for 24 h, and cured in the standard curing room for 7 days. Then, a 20 mm deep notch was cut in the middle span of the length direction of the specimen with the help of a diamond saw blade under a low machining speed and water-cooling conditions.
The experiment was carried out on the WDW-100E universal testing machine. The loading setup is shown in Fig. 4, and the spacing between the two supports is 300 mm. During the experiment, the loading head slowly descends and contacts with the specimen. To fully reduce the gap between the loading head and the specimen, a load of 0.5 kN was applied in advance, followed by stable loading at a constant rate of 0.1 mm/min until the specimen was completely broken. During the loading process, the load and the displacement of the loading head were measured by the built-in sensors of the testing machine.   The load-displacement curves of the fracture process and the crack propagation paths were recorded, as shown in Figs. 5 and 6. Due to inhomogeneity and randomness of the internal structure of concrete, the load-displacement curves are highly dispersed, but they all have obvious three stages, namely, rising stage, sharp falling stage and gentle falling stage. The average peak fracture load is about 6.545 kN. The crack paths are highlighted in red in Fig. 6 for easy observation. It can be seen that, although the crack paths are very tortuous and have significant differences, they all have a tendency to go straight forward. They originate from the tip of the notch, and put their end points near the loading head.

Validation of Numerical Modeling Method
According to Table 1, after excluding fine aggregates with particle sizes less than 4.75 mm, the total areas of packed coarse aggregates with different particle sizes were calculated, as shown in Table 3. Then, the parameterization modeling Yin et al., 2011Yin et al., , 2012Yin et al., , 2015 was adopted to randomly pack 253 aggregates with different particle sizes into a rectangular area with a length of 400 mm and a width of 100 mm. After the geometric modeling was completed, a two-dimensional three-point bending concrete beam model was obtained and then a 20 mm deep notch located in the middle of its bottom, as shown in Fig. 7.
The geometric model was discretized with the triangular plane stress solid elements and two-dimensional 0-thickness cohesive elements. Fig. 8 exhibits the finiteelement mesh with the highlighted cohesive elements in the mortar for (a) and in the ITZ for (b), respectively. The efficiency and accuracy of numerical simulation depend on the finite-element grid density (Zhang et al., 2018b). In general, the denser the grid is, the more accurate the calculation is, but the higher the time cost of the calculation is. To seek a balance between the accuracy and efficiency of the calculation, the element size was changed from 0.75 to 1.00 mm, 1.50 mm, and 2.00 mm. It was proved that the mesh size of 1.50 mm could balance the efficiency and accuracy. In this case, the computational model totally has 75,447 nodes, 37,779 solid elements, 29,704 mortar matrix cohesive elements and 4664 ITZ cohesive elements.    The solid elements and cohesive elements were characterized with the linear elastic model and bilinear softening model, respectively. Their material parameters are shown in Table 4 (Ji et al., 2020;Zhang et al., 2016Zhang et al., , 2018aZhang et al., , 2018bZou et al., 2021). The tensile strengths and fracture energies of mortar and ITZ were obtained through the splitting tensile test (Zou et al., 2021) and three-point bending test on the granite-mortar hybrid specimens and pure mortar specimens (Ji et al., 2020). The tangential parameters of the cohesive elements were assumed to be the same as the corresponding normal parameters (Zhang et al., 2016(Zhang et al., , 2018a(Zhang et al., , 2018b. The loads and constraints are shown in Fig. 9. A displacement load from 0 to 0.6 mm was applied. It includes 10,000 load steps with a constant load increment. Fig. 10 shows the load-displacement curve. For comparison, the envelope region of the experimental load-displacement curves for the six specimens in the same batch is also given in this figure. It can be seen that the numerical load-displacement curve completely fall within the experimental envelope region. The numerical peak load is 5.688 kN, which is about  13.09% lower than the experimental average peak load of 6.545 kN. Considering the following factors, this difference is acceptable: (1) the size, shape and position of the aggregates inside the concrete are random, and the constitutive model parameters of different component materials are approximate, so the numerical model cannot be completely consistent with the physical truth; (2) There are certain errors in the constraint and loading locations, specimen and notch geometries, and so on; (3) A two-dimensional model with the plane stress assumption is used for the numerical simulation; (4) In the process of making and cutting concrete specimens, it is inevitable to introduce initial defects and damage. Fig. 11 shows the experimental and numerical crack propagation paths. As you can see, they have a very good consistency in shape and trend. This indicates that all the used numerical simulation method, constitutive models and their parameters are reasonable. It should be pointed out that coarse aggregate distribution generally has very slight effect on the prepeak mechanical behaviors but significant effect on the post-peak behavior. Considering that the aggregate distribution effect was carefully evaluated Zhang et al., 2016), it is no longer discussed here.

Results and Analyses
In this section, the numerical simulation method was adopted to study the effects of the water-cement ratio of mortar, the roughness of bonding surface and the content of silica fume in interface agent on the tensile fracture properties of concrete.

Computational Model and Conditions
The parameterization modeling Yin et al., 2011Yin et al., , 2012Yin et al., , 2015 was used to establish a two-dimensional concrete model with the size of 100 mm × 100 mm. In the model, coarse aggregate has the total area content of 44.33%, so the packed circular aggregate has the total area content of 49.24%. According to Table 1, the areas of coarse aggregates with different particle sizes were calculated, as shown in Table 5. In the numerical simulation, the material parameters of ITZ and mortar were changed according to the varying mortar water-cement ratio, bonding surface roughness and silica fume content in the interface agent, while the other material parameters remained unchanged. The lateral and longitudinal displacements of the lower surface of the model were constrained, and the upper surface was subjected to a longitudinal displacement increasing from 0 to 0.16 mm at a constant rate of 1.6 mm/s, as shown in Fig. 12. A series of tensile fracture simulations were carried out.

Effect of Mortar Water-Cement Ratio
The water-cement ratio of mortar is an important factor affecting the mechanical properties of concrete. Increasing the water-cement ratio would increase the porosity of ITZ and weaken the aggregate-mortar bonding. The splitting tensile test on the granite-mortar hybrid specimens (Zou et al., 2021) showed that the tensile strength of ITZ is relative to the water-cement ratio of mortar by The unit of the tensile strength is MPa. From the splitting tensile test on the pure mortar specimens, the tensile strength of mortar experimentally depends on the watercement ratio of mortar by Zou et al. (2021) On the other hand, according to the three-point bending test on the granite-mortar hybrid beams, there is the following connection between the ITZ fracture energy and the mortar water-cement ratio (Ji et al., 2020).
(1) f st = −9.296w/c + 5.872 (2) f st = −5.427w/c + 5.730  Ji et al. Int J Concr Struct Mater (2023) 17:2 The unit of the fracture energy is N/m. From the threepoint bending test on the pure mortar beams, the mortar fracture energy is expressed in terms of the mortar water-cement ratio as (Ji et al., 2020) Considering that the water-cement ratio of mortar has a range from 0.3 to 0.5 in the three-point bending test and splitting tensile test, the tensile fracture behaviors of concrete were investigated when the water-cement ratio of mortar was 0.3, 0.35, 0.4, 0.45 and 0.5, respectively. Fig. 13 shows the load-displacement curves for different mortar water-cement ratios. It can be seen that with the decrease of mortar water-cement ratio, both the peak  24.6%, 38.8%, 50.9% and 62.8%, respectively. The influence of water-cement ratio on the post-peak deflection behavior of concrete does not seem consistent with common perception that a lower water-cement ratio generally leads to more brittle concrete. However, it is also noticed from Fig. 13 that this influence is mainly from the impact of water-cement ratio on the peak load. Therefore, the above inconsistency is easy to understand.
Fig. 14 shows the variation of tensile strength of concrete with the water-cement ratio of mortar. The (5) f st = −8.087w/c + 6.840 relationship between them can be established by linear regression analysis as follows. Fig. 15 shows the variation of the dissipated energy of concrete in the process of tensile fracture with the mortar water-cement ratio. The relationship between them can also be established by linear regression analysis as follows.

Effect of Aggregate Surface Roughness
The aggregate surface roughness significantly affects the interfacial strength (Guinea et al., 2002;Santos & Júlio, 2013). The experimental studies (Ji et al., 2020;Rao & Prasad, 2004) showed that the increase of aggregate surface roughness can significantly improve the tensile strength and fracture energy of ITZ. According to the splitting tensile test on the granite-mortar hybrid specimens (Zou et al., 2021) and the three-point bending test on the granite-mortar hybrid beams (Ji et al., 2020), the tensile strength and fracture energy of ITZ have the following connections with the aggregate surface roughness (Ji et al., 2020;Zou et al., 2021).
The tensile fracture behaviors of concrete were simulated for different cases of the aggregate surface roughness. Fig. 16 shows the load-displacement curves for different aggregate surface roughness. It can be seen that the peak load of concrete increases with the increase of (6) W = −9.402w/c + 6.390 (7) f st = 0.849R a + 1.234 (8) G = 73.696R a + 37.094 aggregate surface roughness, and the load-displacement curves start to separate before the load reaches its peak value. With the increase of displacement, their gap becomes more and more prominent. This is mainly because the interface bonding strength increases with the increase of aggregate surface roughness. The higher the bonding strength of the interface is, the greater the energy required for fracture is. As a result, the area enclosed by the load-displacement curve and the horizontal axis is also larger. Figs. 17 and 18, respectively, show the variations of the tensile strength and the dissipated energy with the aggregate surface roughness. Hong et al. (2014) also studied the influence of coarse aggregate surface roughness on the tensile properties of concrete through experiments. Compared to their results, the curves are consistent in shape and trend. Through linear regression analysis, the relations of the tensile strength and dissipated energy of concrete with the surface roughness of coarse aggregate can be established as follows.
Generally speaking, the effects of aggregate surface roughness on the tensile fracture properties of concrete are very complicated. For the aggregate surface roughness between 0 and 0.4 mm, both the tensile strength and dissipated energy of concrete increase rapidly with increasing aggregate surface roughness. This is because when the aggregate surface roughness is small, the mechanical properties of ITZ are much weaker than those of mortar. With the increase of aggregate surface roughness, the mechanical properties of ITZ gradually approach those of mortar, so that the proportion of exposed aggregate area in the fracture surface decreases gradually. Fig. 19ac, respectively, shows the tensile fracture paths of aggregate surface roughness of 0, 0.2 mm and 0.4 mm.
For the aggregate surface roughness from 0.6 to 1.6 mm, with the increase of aggregate surface  Ji et al. Int J Concr Struct Mater (2023) 17:2 roughness, the tensile strength and dissipated energy of concrete still increase, but very slowly. This is because the mechanical properties of ITZ are almost the same as those of mortar. Fig. 19d-f, respectively, shows the tensile fracture paths for aggregate surface roughness of 0.6 mm, 1 mm and 1.6 mm. It can be seen that the exposed aggregates at the fracture appearance is very few, and the area of the exposed aggregates changes very little when the aggregate surface roughness continues to increase. Therefore, the mechanical properties of concrete are no longer sensitive to the change of aggregate surface roughness. Hong et al. (2014) carried out the tensile fracture experiments on the spherical high boric acid aggregate concrete with different roughness, and found that the exposed aggregate area at the fracture appearance for the smooth aggregate concrete is larger and the tensile strength was lower than that for striated aggregate concrete aggregate with surface incision treatment. This is consistent with the above results.

Effect of Silica Fume Content in Interfacial Agent
Due to the pozzolanic ash effect and micro-aggregate effect (Abu Noaman et al., 2019;Sankar et al., 2020), a layer of interfacial agent containing silica fume on aggregate surface can effectively reduce the porosity of ITZ and improve its mechanical properties. The experimental results from the splitting tensile test on the granite-mortar hybrid specimens (Zou et al., 2021) and the threepoint bending test on the granite-mortar hybrid beams (Ji et al., 2020) showed that the tensile strength and fracture energy of ITZ have the following connections with the silica fume content in interfacial agent. Page 12 of 14 Ji et al. Int J Concr Struct Mater (2023) 17:2 with the case without interfacial agent on aggregate surface, the peak loads of concrete for the cases with interfacial agent containing 0, 3%, 6%, 9% and 12% silica fume increase by 15.2%, 15.8%, 15.3%, 13.5% and 6.0%, respectively. When the content of silica fume is about 3%, the peak load is the highest. Fig. 21 shows the variation of the tensile strength of concrete with the silica fume content in the interface agent. It can be seen that the effect of silica fume content on the tensile strength of concrete is complex. When the silica fume content is between 0 and 3%, the tensile strength of concrete increases slightly with the increase of silica fume content, but when the silica fume content is between 3% and 12%, the tensile strength of concrete decreases significantly with the increase of silica fume content. This is qualitatively consistent with the experimental results of high performance concrete (HPC) (Ji et al., 2020). Smarzewski (Ji et al., 2020) found that the tensile strength of HPC increases with growing silica fume content up to 10%, but reduces with the further increase of silica fume content from 10% to 25%, and fitted their relationship with a quadratic parabolic curve. In addition, it is noticed that, different from the additive of silica fume to interface agent in this study, Smarzewski added silica fume to cement paste in his experiments.
Similarly, the relationship between the tensile strength of concrete and silica fume content in interfacial agent can be fitted by the following quadratic parabolic function. Fig. 22 shows the variation of the dissipated energy of concrete complete fracture with the content of silica fume in interfacial agent. It can be seen that the variation of the dissipated energy varies is very similar to the variation of tensile strength with the content of silica fume in interfacial agent. Their relationship can also be fitted by a quadratic parabolic function.

Conclusions
The parameterization modeling was used to establish the model of heterogeneous concrete with aggregate size and distribution information and the bilinear cohesive zone constitutive model was used to describe the cracking behavior of mortar and ITZ. Through mediating the mechanical properties of ITZ based on the experiments, the influences of the water-cement ratio of mortar, the surface roughness of coarse aggregate and the content of silica fume in interfacial agent on the tensile fracture behavior of concrete were quantitatively studied by the numerical simulation. The following conclusions were drawn.
(1) The tensile fracture properties of concrete are closely related to the water-cement ratio of mortar.
With the decrease of water-cement ratio of mortar, both the tensile strength and dissipated fracture energy of concrete increase significantly. They have a negative linear correlation with the water-cement ratio of mortar.
(2) The surface roughness of coarse aggregate has a great influence on the tensile fracture properties of concrete, and both the tensile strength and dissipated fracture energy have a bilinear relationship with the surface roughness of coarse aggregate. When the aggregate surface roughness increases from 0 to 0.4 mm, both the tensile strength and dissipated energy of concrete increase rapidly. When the aggregate surface roughness increases from 0.6 to 1.6 mm, however, the tensile strength and dissipated energy of concrete increase very slowly. (3) The tensile strength and dissipated fracture energy of concrete have an approximate quadratic parabola relationship with the content of silica fume in the interfacial agent. The best tensile fracture properties of concrete can be obtained by wrapping a layer of interfacial agent with a content of 3.5-3.7% of silica fume on the aggregate surface.
In addition, it must be pointed that the validity of the above results and conclusions should be limited within the discussed scopes of the water-cement ratio of mortar, the surface roughness of coarse aggregate and the content of silica fume in interfacial agent in this paper.