Comparisons of Tensile Fracturing Behaviors of Hydraulic Fully Graded and Wet-Screened Concretes: A Mesoscale Study

'ewidely usedwet-screeningmethod 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.'e 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 wetscreening 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 wetscreening, 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.


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].us, 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.erefore, 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 di erences between the mechanical parameters of HFGC and those of the corresponding wet-screened concrete, related experimental comparisons are still sparse due to the practical di culties of performing extensive large-scale tests.Pioneered by Blanks 2 Advances in Materials Science and Engineering and McNamara [12], limited comparative experimental results can be found in [13][14][15][16][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 e ects of curing age, by Serra et al. [19] who studied the in uence of the wet-screening of the dam concrete on the development of the modulus of elasticity based on a speci c 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 di erent strain rates.On the whole, although it has been well accepted that the di erences 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 e ect of these three factors, which can provide fundamental knowledge of the underlying di erentiation 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 ner 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][26][27][28] and the lattice model [29][30][31].In the continuum model, concrete is usually characterized by a continuum composite material with each component discretized by nite 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) nite 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).

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

4
Advances in Materials Science and Engineering a composite consisting of the homogeneous mortar matrix representing cement paste and ne 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 con guration of coarse aggregates, given the size and spatial distributions, volume fraction, and geometric characteristics of the coarse aggregate.
Unlike ordinary concrete usually adopting the wellknown Fuller's curve, the size distribution of coarse aggregates of HFGC is commonly described by three or four grading segments with xed 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  Advances in Materials Science and Engineering 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 threeand four-graded concretes, respectively.While for the wetscreened 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  6 Advances in Materials Science and Engineering directly calculated through dividing the total weight of coarse aggregate used in a unit volume of concrete by its unit weight.However, for the simpli ed 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].e existing studies indicate that the shape of coarse aggregate a ects 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 simpli ed 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 Advances in Materials Science and Engineering 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 A F 50% 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.

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 8 Advances in Materials Science and Engineering 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.

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 nite element software ABAQUS.en, 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 nal FE mesh discretization with ITZ elements highlighted and the zooming image of a part of the FE mesh are depicted in Figure 3.

Constitutive Modeling.
Concerning concrete fracturing, it is well recognized that mesoscale cracking under loading rst 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.us, 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 brie y summarized below.
In the CDP model, two independent hardening variables, that is, equivalent compressive and tensile plastic strains ( ε p c and ε p t ), are employed in order to consider compressive crushing and tensile cracking, respectively.en, two independent damage variables d c ( ε p c ) and d t ( ε p t ) are introduced to characterize the compressive and tensile damage states.In addition, to represent the overall damage in an isotropic manner, a scale variable d is de ned as where s t and s c are the functions of the stress state which are used to represent sti ness recovery e ects associated with stress reversals [37].us, the damaged elastic modulus E related to di erent failure mechanisms under tension and compression can be calculated by where E 0 represents the initial elastic modulus.
Based on the concept of damage mechanics, the e ective stress σ can be obtained as where σ is the Cauchy stress.e yield function of the CDP model is given in the e ective stress space as where ε p [ ε p t • ε p c ] T ; p and q are the e ective hydrostatic pressure and the e ective Mises equivalent deviatoric stress, respectively; σ max is the algebraically maximum eigenvalue of σ; the brackets AEae are used in Macaulay sense; σ c ( ε p c ) is the uniaxial compressive e ective strength; α and c 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 β( ε p ) can be calculated by where σ t ( ε p t ) is the uniaxial tensile e ective strength.Figure 4 illustrates the yield surface in the case of plane stress.

Advances in Materials Science and Engineering
In order to describe the dilatancy reasonably, the nonassociated ow rule is adopted in the CDP model, and the ow potential is formulated as where ϵ is the parameter de ning the rate at which the function approaches the asymptote; σ t0 is the uniaxial tensile stress at failure; ψ is the dilation angle measured in the p − q plane at high con ning pressure.As stated earlier, the material softening under tension is de ned by the relationship between the uniaxial tensile e ective strength and equivalent tensile plastic strain (5), which means mesh sensitivity will be encountered when applying the CDP model in FE simulations.erefore, a stress-displacement relation is used in this study to de ne the tensile softening behavior for alleviating the in uence of mesh sensitivity on the simulation results.

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 e ect inevitably exists, and therefore, its in uence on the solution of a quasistatic problem should be small enough to be neglected.In order to minimise the dynamic e ect, the loading duration should be large enough, while on the other hand, the computational e ort increases proportionally with the increase of loading duration.Hence, a balance has to be made between the computational efciency and simulation accuracy, which can be achieved through comparing the results under di erent loading durations (or loading rates).

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 coe cient of variation (CoV), which can be used as the basis of the subsequent comparative study.

Results and Discussion
e uniaxial tensile fracturing behaviors of hydraulic threegraded concrete with A F 50% 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 wetscreened concrete are compared in detail in this section.

Numerical Specimens and Mechanical Properties.
As stated earlier, the di erences 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 e ects 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  12 Advances in Materials Science and Engineering specimens.us, by comparing the MCS results of Types I and II, the e ect 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 e ect of variation of gradation and volume fraction of coarse aggregate can be individually discussed.e numerical specimens in Types IV, V, and VI employ the same mesostructures as those of the numerical specimens in Type III, and the only di erence 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 90%, while for Types V and VI, the ratios are set to 80% and 70%, respectively.us, the e ect 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 e ect of three factors on the di erences of the investigated macroscopic mechanical parameters can be analyzed, and the governing factor(s) corresponding to each of these parameters can also be identi ed.
For each numerical specimen, uniaxial tensile fracturing simulation is performed.In all FE simulations, the left end of the specimen is xed 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 Advances in Materials Science and Engineering ITZ in the wet-screened concrete are independent and should be de ned 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 75% 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.

E ect of Specimen Size
Variation.e 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 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 t 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 e ect on fracture energy.With respect to each investigated macroscopic parameter, similar probability distributions can be found in Types I and II, Advances in Materials Science and Engineering which indicates the in uence of the specimen size on the probability distributions of these parameters is not evident.Moreover, with the increase of the specimen size, di erent 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 (3.95% and 3.21%, respectively), whereas small MV reductions are shown for peak stress and peak strain (1.76% and 5.48%, 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 e ect 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.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.16 Advances in Materials Science and Engineering

E ect of Variation of Gradation and Volume Fraction of Coarse Aggregate.
e 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 t 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 18.32%), which can be attributed to the existence of large coarse aggregates (much sti er 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 3.65%, 18.73%, and 8.85%, respectively).Overall, the variation of gradation and volume fraction of coarse aggregate is found Advances in Materials Science and Engineering to give rise to big di erences 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 e ects on the variation of mean fracture energy.e 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.
To better understand the underlying di erentiation 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

E ect of the Weaker ITZs Surrounding Large Coarse
Aggregates. e 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 signi cant 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.
e 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 Advances in Materials Science and Engineering 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 su er 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 6.76%, 23.53%, 17.03%, and  Furthermore, to reveal the di erentiation 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 di erent mechanical properties of surrounding large coarse aggregates, which belong to Types IV, V, and VI, respectively, are sketched in Figure 26.Although the long weaker ITZs de nitely 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.

Discussion on the Combined E ect and the Governing Factor(s).
e 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.
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 wetscreened 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 di erence 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 e ect of the weaker ITZs surrounding large coarse aggregates is much more signi cant 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 e ect 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 di erence 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 signi cant change in response to wet-screening, di erent governing factor(s) can be identi ed for di erent mechanical parameters.
Regarding the dispersion of the investigated parameters, it is shown in Figure 28 that the e ects of these factors di er from each other.Speci cally, the e ect of specimen size variation is quite small, and the e ects of the other two factors are relatively larger.
e 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 e ect 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 e ect 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.variation of gradation and volume fraction of coarse aggregate, and the weaker ITZs surrounding large coarse aggregates.e 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) e 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) e 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.

Figure 1 :Figure 2 :
Figure 1: Typical coarse aggregate grading curves for hydraulic fully graded concrete.

Figure 9 :
Figure 9: Probability density of macroscopic mechanical parameters related to tensile fracturing behaviors of Type I. (a) Tensile elastic modulus.(b) Peak stress.(c) Peak strain.(d) Fracture energy.

Figure 10 :
Figure 10: Probability density of macroscopic mechanical parameters related to tensile fracturing behaviors of Type II.(a) Tensile elastic modulus.(b) Peak stress.(c) Peak strain.(d) Fracture energy.

Figure 11 :
Figure 11: In uence of the sample number on CoVs of macroscopic mechanical parameters of Type I. (a) Tensile elastic modulus.(b) Peak stress.(c) Peak strain.(d) Fracture energy.
3, and 0.4, respectively, is illustrated in Figure 2(b).e wet-screened concrete specimens corresponding to the above two examples are shown in Figures 2(c) and 2(d), respectively.

Figure 12 :
Figure 12: In uence of sample number on CoVs of macroscopic mechanical parameters of Type II.(a) Tensile elastic modulus.(b) Peak stress.(c) Peak strain.(d) Fracture energy.

Figure 13 :
Figure 13: Macroscopic stress-strain curves and their mean curve of Type III.

Figure 14 :
Figure 14: Probability density of macroscopic mechanical parameters related to tensile fracturing behaviors of Type III.(a) Tensile elastic modulus.(b) Peak stress.(c) Peak strain.(d) Fracture energy.

Figure 15 :
Figure 15: In uence of thw sample number on CoV of macroscopic mechanical parameters of Type III.(a) Tensile elastic modulus.(b) Peak stress.(c) Peak strain.(d) Fracture energy.

Figure 16 :
Figure 16: Distribution of tensile plastic strain on the mesoscale.At peak stress: (a) Type II; (c) Type III.At the maximum displacement: (b) Type II; (d) Type III.

Figure 17 :Figure 18 :Figure 19 :
Figure 17: Macroscopic stress-strain curves and their mean curve of Type IV.

Figure 20 :
Figure 20: Probability density of macroscopic mechanical parameters related to tensile fracturing behaviors of Type IV.(a) Tensile elastic modulus.(b) Peak stress.(c) Peak strain.(d) Fracture energy.

Figure 21 :
Figure 21: Probability density of macroscopic mechanical parameters related to tensile fracturing behaviors of Type V. (a) Tensile elastic modulus.(b) Peak stress.(c) Peak strain.(d) Fracture energy.

Figure 22 :
Figure 22: Probability density of macroscopic mechanical parameters related to tensile fracturing behaviors of Type VI.(a) Tensile elastic modulus.(b) Peak stress.(c) Peak strain.(d) Fracture energy.

Figure 23 :
Figure 23: In uence of the sample number on CoV of macroscopic mechanical parameters of Type IV.(a) Tensile elastic modulus.(b) Peak stress.(c) Peak strain.(d) Fracture energy.

Figure 24 :
Figure 24: In uence of the sample number on CoV of macroscopic mechanical parameters of Type V. (a) Tensile elastic modulus.(b) Peak stress.(c) Peak strain.(d) Fracture Energy.

Figure 25 :
Figure 25: In uence of the sample number on CoV of macroscopic mechanical parameters of Type VI.(a) Tensile elastic modulus.(b) Peak stress.(c) Peak strain.(d) Fracture energy.

Figure 26 :
Figure 26: Distribution of tensile plastic strain on the mesoscale.At peak stress: (a) Type IV; (c) Type V; (e) Type VI.At the maximum displacement: (b) Type IV; (d) Type V; (f) Type VI.

Table 1 :
Mechanical properties of concrete components.