#### Abstract

The widely used wet-screening method in the experimental testing of hydraulic fully graded concrete inevitably results in a gap between the real mechanical parameters of hydraulic fully graded concrete specimens and those of the corresponding wet-screened specimens and therefore necessitates the comparative study on their mechanical behaviors. To this end, a two-dimensional mesoscale modeling methodology is developed for simulating the tensile fracturing behaviors of hydraulic fully graded and wet-screened concretes, and extensive Monte Carlo simulations are performed. The individual effects of specimen size variation, variation of gradation and volume fraction of coarse aggregates, and the weaker interfacial transition zones surrounding the large coarse aggregates to be removed by wet-screening are detailed followed by the discussion on the combined effect of these three main factors. All the mean values of the macroscopic mechanical parameters related to tensile fracturing behaviors are found to show significant change in response to wet-screening, and the underlying differentiation mechanism and governing factor(s) are identified. Furthermore, it is shown that the randomness of the investigated parameters can be roughly described by the Gaussian distribution, and the dispersion of each of the investigated parameters of hydraulic fully graded concrete is higher than that of the corresponding wet-screened concrete.

#### 1. Introduction

Fully graded concrete is widely used in the construction of hydraulic structures such as gravity dam and arch dam [1]. Compared to ordinary concrete with the maximum size of aggregate (MSA) no more than 40 mm, hydraulic fully graded concrete (hereinafter referred to as HFGC) employs larger coarse aggregates. Typically, the MSA reaches to 80 mm in the case of hydraulic three-graded concrete, while for the four-graded case, the MSA is even larger and increases up to 150 mm (or 120 mm). Moreover, HFGC is also featured by high volume fraction of aggregates composed of fine aggregates and coarse aggregates with the cutoff size between them taken to be 5 mm, low cement content, large admixture dosage, and high water-to-cement ratio [2]. As a result, the mechanical behaviors of HFGC are different from those of the extensively studied ordinary concrete [3]. Thus, to obtain the fundamental mechanical parameters especially the parameters related to fracturing behaviors, which are needed by the failure analysis of hydraulic mass concrete structures, a detailed investigation on the mechanical behaviors of HFGC is required [4, 5].

With regard to identifying the mechanical parameters of concrete, the larger aggregate size implies that a larger size of concrete specimens, typically at least three times of MSA, is needed to fulfill the requirement of statistical representations [6]. According to the Chinese code SL 352-2006 (test code for hydraulic concrete), the length of specimens *L* with cube shape for hydraulic three-graded concrete should be set to 300 mm, while for the four-graded case, *L* is further increased to 450 mm. Compared to the standard specimen with *L* equal to 150 mm prepared for ordinary concrete, the larger specimen for HFGC causes particular challenges when conducting experimental testing in the usual concrete laboratory, including large size of the testing machine, inconvenience of testing operation, and high testing expense, which further hampers to a large extent the direct experimental testing on HFGC specimens [7].

As an alternative to carrying out mechanical experiments on large specimens, the well-known wet-screening method is commonly employed in practice [8]. Following this method, experimental testing is performed on standard specimens, which are cast using the wet-screened concrete obtained by removing the coarse aggregates with the size larger than 40 mm from the original fully graded concrete through sieving, enabling the testing under usual laboratory conditions. Although the wet-screening method can avoid the direct testing on large specimens, a gap inevitably exists [9] between the real mechanical parameters of HFGC specimens and those of the corresponding wet-screened specimens; that is, the standard specimens cast using the wet-screened concrete. Therefore, an extrapolation procedure should be executed with the aim of acquiring the accurate characterization of the mechanical properties of HFGC [10], which necessitates a detailed comparative study on the mechanical parameters of hydraulic fully graded and wet-screened concretes [11].

In spite of the awareness of the differences between the mechanical parameters of HFGC and those of the corresponding wet-screened concrete, related experimental comparisons are still sparse due to the practical difficulties of performing extensive large-scale tests. Pioneered by Blanks and McNamara [12], limited comparative experimental results can be found in [13–17], and the most recent representative investigations were carried out by Deng et al. [18] who presented the behaviors of dam and sieved concretes in uniaxial tension and compression together with the effects of curing age, by Serra et al. [19] who studied the influence of the wet-screening of the dam concrete on the development of the modulus of elasticity based on a specific experimental program using in situ creep cells, and by Shi et al. [20] who performed experimental study on uniaxial compression properties of large aggregates and wet-screened concrete at different strain rates. On the whole, although it has been well accepted that the differences mentioned above are mainly attributed to specimen size variation, variation of gradation and volume fraction of coarse aggregate, and the weaker interfacial transition zones (ITZs) surrounding large coarse aggregate particles in HFGC [21], the individual effect of these three factors, which can provide fundamental knowledge of the underlying differentiation mechanism, has been rarely studied in the past experimental research works. Furthermore, being aware of the fact that the mechanical parameters of concrete are intrinsically random, it is necessary to compare the mechanical parameters of HFGC and the corresponding wet-screened concrete in the statistical sense, which, to the best of our knowledge, has not been conducted.

On the contrary, since the mechanical behaviors of concrete on the structural scale (macroscale) are greatly controlled by its components and their interactions taking place on a finer scale (mesoscale) [22, 23], several mesoscale models have been developed to provide tools for a better understanding of the extremely complicated mechanical behaviors of concrete, especially fracturing [24]. Roughly, there are two types of mesoscale models, namely, the continuum model [25–28] and the lattice model [29–31]. In the continuum model, concrete is usually characterized by a continuum composite material with each component discretized by finite elements, while for the lattice model, a discrete system composed of lattice elements is used to represent concrete. On the whole, both of them are capable of realizing reasonable simulations of microcracking, coalescence of multiple distributed cracks into localized cracks, and fracture propagation, provided each component in concrete is modeled with use of a well-suited mechanical model. In contrast to the continuum model, the lattice model is considered computationally less demanding as concrete mesostructure is roughly represented by a discrete system with relatively less degrees of freedom and meanwhile can still possess the ability to capture the most important aspects of concrete fracturing. However, it is hard to investigate the interactions of concrete components in a real sense as the actual concrete mesostructure is not fully taken into account. Consequently, mesoscale numerical simulation using the continuum model can be considered as a promising complement to experimental testing in the comparative study on the mechanical parameters of HFGC and wet-screened concrete.

Considering this, a two-dimensional (2D) finite element (FE) mesoscale modeling framework for hydraulic fully graded and wet-screened concretes is proposed in this study, in which concrete is considered as a three-phase composite composed of coarse aggregate, mortar, and ITZ, and comprehensive comparisons of tensile fracturing behaviors of HFGC and the wet-screened concrete are performed based on extensive Monte Carlo simulations (MCS).

#### 2. Mesostructure Generation of HFGC and Wet-Screened Concrete

To explicitly model the components and their interactions of both HFGC and the wet-screened concrete on the mesoscale, the internal material structure of concrete should be generated. In this regard, concrete is considered to be a composite consisting of the homogeneous mortar matrix representing cement paste and fine aggregates, randomly distributed coarse aggregates, and ITZs surrounding coarse aggregates. For a specimen, the key issue in generating its internal material structure is to determine the configuration of coarse aggregates, given the size and spatial distributions, volume fraction, and geometric characteristics of the coarse aggregate.

Unlike ordinary concrete usually adopting the well-known Fuller’s curve, the size distribution of coarse aggregates of HFGC is commonly described by three or four grading segments with fixed maximum and minimum aggregate sizes for each segment. With respect to hydraulic three-graded concrete, coarse aggregates are divided into three groups: the small coarse aggregate with size ranging from 5 to 20 mm, the medium coarse aggregate with size ranging from 20 to 40 mm, and the large coarse aggregate with size ranging from 40 to 80 mm, while regarding the four-graded case, another group of coarse aggregate with even larger size ranging from 80 to 150 mm (or 120 mm), namely, the extra-large coarse aggregate, is used besides the above three groups. In addition, it is usually assumed that the size of the coarse aggregate in a certain segment follows the uniform distribution as the size distribution of the coarse aggregate for each segment is commonly not precisely controlled during the production of HFGC. Hence, on the basis of the given coarse aggregate weight of each grading segment, the size distribution of coarse aggregate can be determined with ease. Figure 1 depicts two typical coarse aggregate grading curves corresponding to hydraulic three- and four-graded concretes, respectively. While for the wet-screened concrete, the coarse aggregate gradation can be calculated in a straightforward way, by removing grading segment(s) with size larger than 40 mm from the coarse aggregate gradation of HFGC.

Concerning the spatial distribution of coarse aggregates, it is well accepted that coarse aggregates can be considered to be randomly distributed in a certain specimen. Consequently, the location of an individual coarse aggregate needed to be placed in the specimen, which can be represented by the coordinates of its geometric center, is assumed to be uniformly distributed throughout the specimen in this study.

For a given concrete mix design, the volume fraction of coarse aggregate in the three-dimensional (3D) sense can be directly calculated through dividing the total weight of coarse aggregate used in a unit volume of concrete by its unit weight. However, for the simplified two-dimensional case, a conversion of the volume fraction of coarse aggregate from 3D to 2D should be conducted. To this end, the well-known Walraven’s conversion equation used by several researchers is employed in the present work [32].

The existing studies indicate that the shape of coarse aggregate affects the mechanical behaviors of concrete, especially the mesoscale fracture mechanism [24]. Regarding HFGC, there are mainly two types of coarse aggregates: the gravel stone and the crushed stone. In the 2D case, an ellipse can be utilized to represent the gravel stone, while the crushed stone can be simplified to a polygon. Due to the fact that the crushed stone is more commonly used in practice compared to the gravel one, the shape of the coarse aggregate in this study is modeled by the polygon.

After determining the gradation and volume fraction of the coarse aggregate along with the shape and size of the concrete specimen, the required concrete mesostructure can be randomly generated by employing the widely used parameterized modeling approach based on the take-and-place method [33]. In this study, the versatile mesostructure generator for concrete (MGC), developed using MATLAB, is used to generate the following required HFGC specimens and the wet-screened concrete specimens, and the detailed procedure and implementation of MGC can be referred to our previous work [34]. Figure 2(a) sketches an example of a hydraulic four-graded concrete specimen with the coarse aggregate volume fraction and the proportions of the small, medium, large, and extra-large coarse aggregates set to 0.25, 0.25, 0.2, and 0.3, respectively, while an example of a hydraulic three-graded concrete specimen with the same coarse aggregate volume fraction and the proportions of the small, medium, and large coarse aggregates set to 0.3, 0.3, and 0.4, respectively, is illustrated in Figure 2(b). The wet-screened concrete specimens corresponding to the above two examples are shown in Figures 2(c) and 2(d), respectively.

#### 3. Finite Element Modeling Methodology and Monte Carlo Simulations

Provided the mesostructure of a concrete specimen, the corresponding computational model is needed for performing the following mesoscale study. A 2D mesoscale FE modeling methodology for tensile fracturing simulations is developed in this section. Moreover, the Monte Carlo method is employed to take into account the randomness of concrete mesostructure and to obtain statistical characteristics of the macroscopic mechanical parameters related to tensile fracturing behaviors, serving as the basis of the subsequent comparative study.

##### 3.1. Mesh Discretization

In order to realize the precise mesh discretization of the generated HFGC and wet-screened concrete specimens, a two-step meshing approach is proposed, by which the accurate representation of concrete mesostructure in the FE model can be achieved. Firstly, an original FE mesh including coarse aggregates and mortar, rigidly obeying the given arrangement of coarse aggregates, is generated by exploiting the powerful preprocessing modules provided by the commercial finite element software ABAQUS. Then, the thin-layer four-node elements with a uniform thickness approximately set to 100 *μ*m according to the experimental observation [35], which are used to model ITZs, are automatically inserted between coarse aggregate elements and their surrounding mortar elements. Following this strategy, a mesh generator is developed using MATLAB and Python, with the detailed procedure and implementation described in the previous study [34]. Moreover, in order to ensure the accurate representation of the geometric characteristics of coarse aggregates and meanwhile to generate a computationally less demanding FE model, the distance between two neighbouring seeds along the edges of coarse aggregates and the boundary of the specimen, which controls the average element size (or mesh density) to a great extent, is set to 0.4 times of the minimum size of the coarse aggregate. An example of the final FE mesh discretization with ITZ elements highlighted and the zooming image of a part of the FE mesh are depicted in Figure 3.

##### 3.2. Constitutive Modeling

Concerning concrete fracturing, it is well recognized that mesoscale cracking under loading first appears in ITZs owing to their weaker properties. Afterwards, the existing cracks propagate into mortar and additional cracks may initiate within mortar during the process of further loading, while coarse aggregates commonly behave elastically. Thus, the mechanical behavior of coarse aggregates is simulated herein by the isotropic linear elastic model, whereas a continuum damaged plasticity (CDP) model implemented in ABAQUS [36] is used to describe the mechanical behaviors of both mortar and ITZ, which is briefly summarized below.

In the CDP model, two independent hardening variables, that is, equivalent compressive and tensile plastic strains ( and ), are employed in order to consider compressive crushing and tensile cracking, respectively. Then, two independent damage variables and are introduced to characterize the compressive and tensile damage states. In addition, to represent the overall damage in an isotropic manner, a scale variable is defined aswhere and are the functions of the stress state which are used to represent stiffness recovery effects associated with stress reversals [37].

Thus, the damaged elastic modulus related to different failure mechanisms under tension and compression can be calculated bywhere represents the initial elastic modulus.

Based on the concept of damage mechanics, the effective stress can be obtained aswhere is the Cauchy stress.

The yield function of the CDP model is given in the effective stress space aswhere ; and are the effective hydrostatic pressure and the effective Mises equivalent deviatoric stress, respectively; is the algebraically maximum eigenvalue of ; the brackets are used in Macaulay sense; is the uniaxial compressive effective strength; and are the dimensionless material constants, which can be determined by comparing the initial equibiaxial and uniaxial compressive yield stress and by comparing the yield conditions along the tensile and compressive meridians, respectively; and can be calculated bywhere is the uniaxial tensile effective strength. Figure 4 illustrates the yield surface in the case of plane stress.

In order to describe the dilatancy reasonably, the nonassociated flow rule is adopted in the CDP model, and the flow potential is formulated aswhere is the parameter defining the rate at which the function approaches the asymptote; is the uniaxial tensile stress at failure; is the dilation angle measured in the plane at high confining pressure.

As stated earlier, the material softening under tension is defined by the relationship between the uniaxial tensile effective strength and equivalent tensile plastic strain (5), which means mesh sensitivity will be encountered when applying the CDP model in FE simulations. Therefore, a stress-displacement relation is used in this study to define the tensile softening behavior for alleviating the influence of mesh sensitivity on the simulation results.

##### 3.3. Numerical Solution Algorithm

Owing to the highly nonlinear and softening behavior of concrete in the process of tensile fracturing, the ABAQUS/Explicit solver is employed in the present work in order to capture the entire fracturing process.

Regarding the explicit FE modeling, the dynamic effect inevitably exists, and therefore, its influence on the solution of a quasistatic problem should be small enough to be neglected. In order to minimise the dynamic effect, the loading duration should be large enough, while on the other hand, the computational effort increases proportionally with the increase of loading duration. Hence, a balance has to be made between the computational efficiency and simulation accuracy, which can be achieved through comparing the results under different loading durations (or loading rates).

##### 3.4. Monte Carlo Simulations and Statistical Analysis

Due to the random arrangement of coarse aggregates, the mechanical behaviors of concrete obtained from both experimental testing and numerical modeling vary with specimens having the same size, especially in the nonlinear regime, which consequently necessitates the analysis of concrete mechanical behaviors in a statistical sense in order to get more conclusive and accurate results. Hence, the Monte Carlo simulations are conducted in this study aiming at obtaining the statistical characteristics of the mechanical parameters related to tensile fracturing behaviors of both HFGC and wet-screened concrete.

Within the context of MCS, the sample number (the number of concrete specimens), which controls the accuracy of statistical characteristics and therefore should be large enough in order to meet the requirement of statistical convergence, can be taken as the smallest one needed to stabilize the values of statistical parameters with respect to the number of concrete specimens. On the contrary, for each specimen in the MCS, the corresponding macroscopic mechanical parameters related to concrete tensile fracturing behaviors, including tensile elastic modulus, peak stress (tensile strength), peak strain, and fracture energy, can be determined with ease based on the uniaxial tensile stress-strain (or displacement) curve obtained from the mesoscale FE modeling, and then, statistical analysis can be performed for each parameter mentioned above, providing the statistical characteristics of these parameters, including mean value (MV), standard deviation (SD), and coefficient of variation (CoV), which can be used as the basis of the subsequent comparative study.

#### 4. Results and Discussion

The uniaxial tensile fracturing behaviors of hydraulic three-graded concrete with in the sense of 2D and the proportions of small, medium, and large coarse aggregate set to 0.3, 0.3, and 0.4, respectively, and the corresponding wet-screened concrete are compared in detail in this section.

##### 4.1. Numerical Specimens and Mechanical Properties

As stated earlier, the differences in macroscopic mechanical parameters related to tensile fracturing behaviors of HFGC and the wet-screened concrete specimens are mainly attributed to three factors: specimen size variation, variation of gradation and volume fraction of coarse aggregate, and the weaker ITZs surrounding large coarse aggregates. In order to study the effects of these three factors, six types of numerical concrete specimens are generated in this study. For each type of numerical specimens, MCS is performed, and the specimen number is taken as 50 to satisfy the requirement of statistical convergence.

Numerical specimens in Type I with standard dimensions of 150 mm 150 mm are used to represent the wet-screened concrete specimens, while numerical specimens in Type II with dimensions of 300 mm 300 mm are generated to represent the wet-screened concrete specimens with the same size as the hydraulic three-graded concrete specimens. Thus, by comparing the MCS results of Types I and II, the effect of specimen size variation can be extracted. Furthermore, by conducting the MCS of Type III consisting of the hydraulic three-graded concrete specimens with dimensions of 300 mm 300 mm and with all ITZs assumed to having the same mechanical parameters and comparing the MCS results of Types I and II, the effect of variation of gradation and volume fraction of coarse aggregate can be individually discussed. The numerical specimens in Types IV, V, and VI employ the same mesostructures as those of the numerical specimens in Type III, and the only difference lies in the mechanical properties of the ITZs surrounding large coarse aggregates (to be removed in the wet-screening process). For Type IV, the ratio of the mechanical properties of ITZs surrounding coarse aggregates with size larger than 40 mm to those of ITZs surrounding coarse aggregates with size no more than 40 mm is set to , while for Types V and VI, the ratios are set to and , respectively. Thus, the effect of the long weaker ITZs can be roughly evaluated by comparing the MCS results of Types III, IV, V, and VI. Finally, by comparing the MCS results of Types I, II, III, IV, V, and VI, the combined effect of three factors on the differences of the investigated macroscopic mechanical parameters can be analyzed, and the governing factor(s) corresponding to each of these parameters can also be identified.

For each numerical specimen, uniaxial tensile fracturing simulation is performed. In all FE simulations, the left end of the specimen is fixed in the horizontal direction, while the opposite end is subjected to a uniformly distributed horizontal displacement up to 0.18 mm corresponding to 1200 microstrains for Type I or 600 microstrains for Types II, III, IV, V, and VI, namely, a displacement-controlled loading scheme is used. Following the strategy discussed in Section 3.3, the loading time is set to 0.036 s, which corresponds to a loading rate 5 mm/s.

According to the setup of numerical specimen type, only the mechanical properties of coarse aggregate, mortar, and ITZ in the wet-screened concrete are independent and should be defined in the mesoscale FE simulations. However, due to the lack of experimental results for ITZ, the compressive and tensile strengths and elastic modulus of ITZ are assumed to be of those of mortar since ITZ is considered to be weaker than mortar, while the other mechanical properties of ITZ are taken as the same as those of mortar. Table 1 lists the mechanical properties of coarse aggregate, mortar, and ITZ adopted in this study. It is noted that the mechanical properties of mortar is directly obtained from the Chinese code GB 50010-2002 (the code for designing concrete structures), and the compression hardening curve and the tension softening curve are shown in Figures 5 and 6, respectively.

##### 4.2. Effect of Specimen Size Variation

The macroscopic stress-strain curves of all 50 specimens of Type I with the mean curve are plotted in Figure 7, and Figure 8 depicts the corresponding results of Type II. Overall, it can be observed that all the concrete specimens approximately exhibit linear elastic responses on the macroscale in the prepeak stage, whereas in the postpeak phase, the macroscopic stress decreases in a nonlinear way (softening) with the increase of applied displacement, which can be attributed to the mesoscopic cracking taking place in ITZs and mortar. Furthermore, it is worth noting that although the elastic responses of all the specimens of Type I or Type II are quite close to each other, dispersion clearly appears during tensile softening, indicating the necessity of MCS.

With the peak stress, peak strain and fracture energy extracted from each stress-strain (displacement) curve, the statistical analysis of the investigated macroscopic parameters can be executed. Figure 9 presents the probability densities and the best fit Gaussian probability density functions (PDFs) of the tensile elastic modulus, peak stress, peak strain, and fracture energy of Type I, together with MVs, SDs, and CoVs, whereas the statistical results of Type II are depicted in Figure 10. It is shown that the randomness of the investigated parameters can be roughly characterized by the Gaussian distribution. However, fracture energy is distinguished from elastic modulus, peak stress, and peak strain by its large CoV, which suggests that the random mesostructure has a bigger effect on fracture energy. With respect to each investigated macroscopic parameter, similar probability distributions can be found in Types I and II, which indicates the influence of the specimen size on the probability distributions of these parameters is not evident. Moreover, with the increase of the specimen size, different variations of the MVs of the investigated parameters can be observed. For the tensile elastic modulus and fracture energy, their MVs slightly increase with increasing specimen size ( and , respectively), whereas small MV reductions are shown for peak stress and peak strain ( and , respectively). On the whole, it can be concluded that the increase of the specimen size does result in the variations of the macroscopic mechanical parameters related to tensile fracturing behaviors, which is consistent with the existing research results concerning the size effect of concrete [38]. However, in view of the large variation amplitudes of the investigated parameters obtained by limited experimental studies [11] (e.g., the typical ratio of the tensile strength of HFGC specimen to that of the corresponding wet-screened concrete specimen is ranging from 0.65 to 0.85), the variation of the specimen size is not a governing factor.

**(a)**

**(b)**

**(c)**

**(d)**

**(a)**

**(b)**

**(c)**

**(d)**

Figures 11 and 12 illustrate the variations of CoVs of tensile elastic modulus, peak stress, peak strain, and fracture energy with the sample number for Types I and II, respectively, from which it is noted that 50 random specimens are enough to achieve statistical convergence.

**(a)**

**(b)**

**(c)**

**(d)**

**(a)**

**(b)**

**(c)**

**(d)**

##### 4.3. Effect of Variation of Gradation and Volume Fraction of Coarse Aggregate

The macroscopic stress-strain curves of all 50 specimens of Type III are drawn in Figure 13, along with the mean curve, and it is shown that the characteristic of the curves is similar to those of Type II.

Figure 14 demonstrates the probability densities and the best fit Gaussian probability density functions (PDFs) of the tensile elastic modulus, peak stress, peak strain, and fracture energy of Type III, with the calculated MVs, SDs, and CoVs. Also, it can be found that the randomness of the macroscopic parameters can be generally described by the Gaussian distribution, and the CoV of fracture energy is much larger than that of tensile elastic modulus, peak stress, and peak strain. Compared to Type II, all calculated CoVs of Type III tend to decrease, which may be due to the increase of coarse aggregate volume fraction. For tensile elastic modulus, it is shown that its MV of Type III is clearly bigger than that of Type II (increased by ), which can be attributed to the existence of large coarse aggregates (much stiffer than mortar and ITZ) in HFGC. On the contrary, the MVs of peak stress, peak strain, and fracture energy of Type III are smaller than those of Type II (reduced by , , and , respectively). Overall, the variation of gradation and volume fraction of coarse aggregate is found to give rise to big differences in the investigated parameters except peak stress and consequently can be considered to play an important role in the variations of tensile elastic modulus, peak strain, and fracture energy. Additionally, it is worth noting that specimen size variation and the variation of gradation and volume fraction of coarse aggregate have opposite effects on the variation of mean fracture energy.

**(a)**

**(b)**

**(c)**

**(d)**

The variations of CoVs of elastic modulus, peak stress, peak strain, and fracture energy with the number of samples for Type III are shown in Figure 15, from which it can also be concluded that the adopted sample number meets the requirement of statistical convergence.

**(a)**

**(b)**

**(c)**

**(d)**

To better understand the underlying differentiation mechanism on the mesoscale with respect to gradation and volume fraction, the elements with nonzero tensile plastic strain at peak stress and the maximum displacement of two typical specimens (one for Type II and the other for Type III) are visualized in Figure 16, in which the element with equivalent tensile plastic strain bigger than 100 microstrains is highlighted in red. It can be found that the majority of tensile cracks on the mesoscale is initiated in ITZs, and the subsequent propagation of mesoscale cracking is prone to take place and develop around the end of long ITZs (surrounding large coarse aggregates) due to their lower capacity of resisting fracture compared to the short ITZs (around other coarse aggregates), leading to the lower peak stress of Type III. Moreover, it is shown that ITZs account for a bigger proportion of the macroscopic crack in the case of Type III compared to the Type II case, resulting in the lower fracture energy of Type III.

**(a)**

**(b)**

**(c)**

**(d)**

##### 4.4. Effect of the Weaker ITZs Surrounding Large Coarse Aggregates

The macroscopic stress-strain curves with respect to Types IV, V, and VI, along with the corresponding mean curves, are illustrated in Figures 17–19, respectively. It can be observed that the weakened mechanical properties of ITZs surrounding large coarse aggregates do not lead to significant change of characteristic of the macroscopic stress-strain relation, which can still be roughly divided into two stages: the prepeak linear stage and the postpeak softening stage.

The statistical results of the tensile elastic modulus, peak stress, peak strain, and fracture energy of Types IV, V, and VI are given in Figures 20–22, respectively. It is shown that the Gaussian distribution generally remains suitable for describing the randomness of these macroscopic parameters, and the dispersion of fracture energy is still much higher than that of the other three parameters. Moreover, the dispersion of the investigated parameters tends to increase as the mechanical properties of the long ITZs gradually weakened, whereas tensile elastic modulus, peak stress, peak strain, and fracture energy all suffer decreased MV. Compared to Type III, the MVs of tensile elastic modulus, peak stress, peak strain, and fracture energy of Type VI are decreased by , , , and , respectively. On the whole, all the investigated parameters show significant change in response to the weaker ITZs surrounding large coarse aggregates and consequently should be considered to play a key role in the differences in the investigated parameters.

**(a)**

**(b)**

**(c)**

**(d)**

**(a)**

**(b)**

**(c)**

**(d)**

**(a)**

**(b)**

**(c)**

**(d)**

Again, it can be observed from Figures 23–25 that the number of specimens (50) is sufficient to reach converged results for Types IV, V, and VI.

**(a)**

**(b)**

**(c)**

**(d)**

**(a)**

**(b)**

**(c)**

**(d)**

**(a)**

**(b)**

**(c)**

**(d)**

Furthermore, to reveal the differentiation mechanism induced by the weaker ITZs surrounding large coarse aggregates, the elements with nonzero tensile plastic strain at peak stress and the maximum displacement of three specimens with the same mesostructure as the one illustrated in Figures 16(c) and 16(d) but different mechanical properties of ITZs surrounding large coarse aggregates, which belong to Types IV, V, and VI, respectively, are sketched in Figure 26. Although the long weaker ITZs definitely give rise to the decrease of tensile elastic modulus, the resulting lower load bearing capacity still reaches to its peak value (represented by peak stress) at a smaller macroscopic strain (peak strain), which is responsible for fewer aggregate-mortar interfacial cracks on the mesoscale (Figures 16(c), 26(a), 26(c), and 26(e)). On the contrary, as the long ITZs become weaker, the number of mesoscale cracks at the maximum displacement gradually decreases (Figures 16(d), 26(b), 26(d), and 26(f)), which means the existence of long weaker ITZs will suppress to a certain degree the development of mesoscale cracking in other regions and consequently leads to less energy dissipation.

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

##### 4.5. Discussion on the Combined Effect and the Governing Factor(s)

The mean stress-strain curves of Types I, II, III, IV,V, and VI are plotted in Figure 27, and the PDFs of these types with respect to tensile elastic modulus, peak stress, peak strain, and fracture energy, with MVs, SDs, and CoVs, are shown in Figure 28.

**(a)**

**(b)**

**(c)**

**(d)**

Even though the consideration of the long weaker ITZs in HFGC results in the decrease of its tensile elastic modulus, it is shown in Figure 27 that the mean tensile elastic modulus of HFGC is still bigger than that of the corresponding wet-screened concrete. Furthermore, since the increase of tensile elastic modulus caused by increased specimen size is relatively small, the variation of gradation and volume fraction of coarse aggregate should be viewed as the governing factor of the difference in the mean tensile elastic modulus. Unlike tensile elastic modulus, the mean peak stress of HFGC is much smaller than that of the corresponding wet-screened concrete, and all the three factors are responsible for this observation in the same direction. However, the effect of the weaker ITZs surrounding large coarse aggregates is much more significant than that of the other two factors and therefore plays a governing role. Similarly, the mean peak strain of HFGC is much smaller compared to that of the corresponding wet-screened concrete, and each of the three factors plays a positive role in this regard. As the effect of increased specimen size on the decrease of mean peak strain is small, the other two factors, especially the variation of aggregate gradation and coarse aggregate content, dominate the difference in the mean peak strain. In addition, although the increased specimen size leads to a small increase in the mean fracture energy (obtained from the stress-displacement curve rather than the stress-strain curve shown in Figure 27), it is found that the mean fracture energy of HFGC is lower than that of the corresponding wet-screened concrete. Since either the variation of aggregate gradation and coarse aggregate content or the long weaker ITZs decrease the mean fracture energy to a relatively large extent, these two factors are viewed to be dominant in this respect. On the whole, although the MVs of tensile elastic modulus, peak stress, peak strain, and facture energy caused by the wet-screened method indeed show significant change in response to wet-screening, different governing factor(s) can be identified for different mechanical parameters.

Regarding the dispersion of the investigated parameters, it is shown in Figure 28 that the effects of these factors differ from each other. Specifically, the effect of specimen size variation is quite small, and the effects of the other two factors are relatively larger. The dispersions of the investigated parameters of HFGC are found to be lower than those of the corresponding wet-screened concrete when taking into consideration only the variation of gradation and volume fraction of coarse aggregate, while an opposite effect of the weaker ITZs surrounding large coarse aggregates on the dispersions can be observed. In view of the fact that the dispersion of each of the investigated parameters of HFGC is clearly higher than that of the corresponding wet-screened concrete due to the combined effect of three main factors, the weaker ITZs surrounding large coarse aggregates are considered to be the key factor controlling the dispersion variations of the investigated parameters.

#### 5. Conclusions

A finite element mesoscale modeling methodology for hydraulic fully graded and wet-screened concretes is first proposed in the present work. Then, based on the mesoscale FE modeling and extensive Monte Carlo simulations, a detailed comparative study on the tensile fracturing behaviors of fully graded and wet-screened concretes is carried out with prime attention placed on the individual and combined effects of three main factors, namely, specimen size variation, variation of gradation and volume fraction of coarse aggregate, and the weaker ITZs surrounding large coarse aggregates. The work results in the following conclusions:(i)It is shown that the tensile fracturing behaviors of hydraulic fully graded and wet-screened concretes are clearly sensitive to the internal material structure, which necessitates the extensive Monte Carlo simulations performed in this study.(ii)All the investigated macroscopic mechanical parameters related to tensile fracturing behaviors, that is, tensile elastic modulus, peak stress, peak strain, and fracture energy, show significant change in response to wet-screening, even though the characteristics of the stress-strain (or displacement) relation of hydraulic fully graded concrete are similar to those of the wet-screened concrete.(iii)The mean tensile elastic modulus of hydraulic fully graded concrete is found to be larger than that of the corresponding wet-screened concrete, which is mainly due to the variation of gradation and volume fraction of coarse aggregate.(iv)Owing to the existence of the long weaker ITZs in hydraulic fully graded concrete, the mean peak stress of hydraulic fully graded concrete is smaller than that of the corresponding wet-screened concrete.(v)Both the variation of gradation and volume fraction of coarse aggregate and the long weaker ITZs are considered to lead to a smaller mean peak strain of hydraulic fully graded concrete compared to the corresponding wet-screened concrete.(vi)It is found that the mean fracture energy of hydraulic fully graded concrete is lower than that of the corresponding wet-screened concrete and the variation of gradation and volume fraction of coarse aggregate and the long weaker ITZs are identified to be the governing factors.(vii)The degree of dispersion of the investigated parameters of HFGC is found to be higher than that of the corresponding wet-screened concrete, and the long weaker ITZs can be considered to play a dominant role in this respect.

#### Data Availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

This work was supported by projects of the National Natural Science Foundation of China (Grant nos. 51739006 and U1765204) and by “the Fundamental Research Funds for the Central Universities” (Grant no. 2018B11414).