5.1. Determining the Mesoscale Parameters
The discrete element simulation of recycled lump concrete involves modeling the following eight kinds of contacts:
New aggregate contacting new aggregate (denoted as new aggregate–new aggregate contact);
New mortar–new mortar contact;
New aggregate–new mortar contact;
Old aggregate–old aggregate contact;
Old mortar–old mortar contact;
Old aggregate–old mortar contact;
Old aggregate–new mortar contact;
Old mortar–new mortar contact.
The former six contacts describe the new and old concrete. The approach for determining their mesoscale parameters is similar to that illustrated in
Table 4. The seventh and eighth types of contacts mainly describe the cohesion between new concrete and the old concrete. The mesoscale parameters employed in describing those contacts are based on the following rules:
(i) The seventh contact mainly reflects the bonding between old aggregate and the new mortar. There should be no notable difference between this type of contact and the aggregate–mortar contact in the new concrete. Therefore, the parameters employed in the seventh contact are the same as those used to represent the third (new aggregate–new mortar) contact.
(ii) The cohesive behavior of the ITZs between the old and new mortar can refer to published test results describing old mortar–new mortar interfaces.
A group led by Xiao published [
13] the results of nanoscale indentation tests on recycled aggregate concrete. They showed that the properties of the old mortar–new mortar interface are tightly linked with those of the new mortar. As the microscale modulus of the new mortar varies, the micro-modulus of the old mortar–new mortar ITZs changes in tandem, maintaining the ratio of the ITZ’s modulus to that of the new mortar’s close to 85%. Liu and his colleagues also published [
57] the results of similar nanoscale indentation tests. They found that the microscale hardness of the old mortar–new mortar ITZs is the same as the hardness of the new mortar. They attributed this result to the contribution of the recycled coarse aggregate’s characteristics. Recycled aggregate is much rougher than typical natural aggregate, and its surface has many more microscale pores. Both of those features enhance the bonding between the new mortar and old mortar.
In this study’s simulations, the macroscale properties of the new mortar–old mortar contacts were conservatively assumed to be 85% of those of new mortar–new mortar contacts.
Table 7 summarizes the detailed parameters of the eight types of contacts used in the discrete element model of the recycled lump concrete.
5.2. Crack Development
Traditional measuring techniques have difficulty expressing the initialization and propagation of cracks within concrete, but the DEM can provide a primary representation. In this study, three recycled lump concrete samples under axial compressive loading were simulated. The DCL replacement ratio simulated was 33%, and the characteristic lump size was 75 mm. The compressive strengths of the new and old concrete represented varied, but the strength combinations were as follows:
(1) The strength of the new concrete fpr, new was taken as 45.29 MPa, distinctly higher than that of the old concrete fpr, old—28.43 MPa;
(2) The compressive strengths of both the new and old concrete were taken as 45.29 MPa;
(3) The new concrete’s strength was taken as 28.43 MPa, distinctly lower than that of the old concrete—45.29 MPa.
The details of the mesoscale parameters of the new and old concrete are illustrated in
Section 4.2.1.
Figure 10a,c present the numerical crack development process around a waste concrete lump when the applied load lies within 30% to 60% of the peak load. The crack distribution for the whole concrete sample at 40% of the bearing capacity given by the DEM is also shown in
Figure 10d. The combined strength of the new and old concrete clearly influenced the recycled lump concrete’s behavior in uniaxial compression. When the new concrete was stronger than the old, most of the damage during the initial loading was located at the waste concrete lumps, especially localized at the old aggregate–old mortar interfaces. Failure of the interfaces between new and old concrete also happened, but its amount was relatively small in terms of the overall damage. Thus, the interfaces between the new and old concrete were not an obvious weak point in this situation.
In the case where the new concrete’s strength was equal to that of the old concrete, the crack distribution was relatively uniform over the entire specimen. The amount of damage in the old concrete lumps was much less than in the first case, but the damage in the new concrete was more serious.
When the new concrete was distinctly weaker than the old, the microscale damage was localized in the new concrete, especially at the new aggregate–new mortar interfaces. As compared with the former two strength combinations, the amount of contact failure at the ITZs between the new and old concrete was significantly greater.
5.3. Influence of the Distribution of DCLs and Coarse Aggregates
The influence of the spatial distribution of DCLs and aggregates is discussed in this section using the mesoscale multiphase models together with the discrete element method. However, before studying such an effect, it is vital to determine the throw procedure times. This is because the statistical characteristics of the concrete’s mechanical properties are unstable if too few random throws are simulated. However, if too many random throws are simulated, the computation cost is exorbitant.
Figure 11 lists the average peak strengths, the stresses at axial strains of 2500 and 3000
με, and their standard deviations against the number of random samples. All of the samples in the figure include ordinary concrete, as well as recycled lump concrete with replacement ratios of 11%, 22%, 33%, and 44%. For the recycled lump concrete, the compressive strengths of new and old concrete were taken as 45.29 and 28.43 MPa, respectively (here, it indicates that the new and old concretes adopt the mesoscale parameters of the concretes with the strengths of 45.29 MPa and 28.43 MPa in
Table 5). As
Figure 11 shows, beyond 20 samples, the concrete average strength and its standard deviation changed only slightly. Thus, the number of random throws was set at 30 in these simulations.
A total of 420 concrete samples were simulated to study the influence of the random distribution of DCLs and coarse aggregates on recycled lump concrete’s performance in compression. The main variables were the relative strength of the new and old concrete, and the lump replacement ratio.
Table 8 presents the details of these variables.
Figure 12 displays the calculated compressive stress–strain curves of several recycled lump concrete samples. In these samples, the compressive strength of the new concrete was obviously higher than that of old concrete. To better express the influence of the random distribution of DCLs and coarse aggregates on the stress–strain curve, the figure also presents the average stress at full load and its coefficient of variation (COV).
As can be seen from
Figure 12, the mean compressive strength and elastic modulus of the recycled lump concrete decreased significantly with increasing lump replacement. This agrees with the observations reported from previous studies [
3]. The main reason was due to the incorporation of the weak old concrete. Additionally, it is also easy to see that the variation trend in the stress–strain curve of the recycled lump concrete caused by the random distribution of DCLs and aggregates was similar to that of ordinary concrete induced by the coarse aggregates alone.
(1) Up to 90% of the peak load, the stress–strain curves of the various samples coincided completely, indicating that the random distribution of mesoscale phases affected the compressive performance of the concrete only slightly. It also may be deduced that the elastic modulus was nearly uninfluenced by the spatial locations of the DCLs and coarse aggregate pieces.
(2) As the applied load increased further, an impact of the random distribution of DCLs and coarse aggregates gradually became apparent. The variability of the stress at a specific axial strain also increased.
(3) When the axial load reached 50% to 60% beyond the peak load, the variability of the concrete’s stress–strain curve increased notably and, with continued loading, it continued to increase quickly.
Table 9 shows the statistical results of the peak strength and the stresses at axial strains of 2500 and 3000
με for all of the concrete samples listed in
Figure 12. Two findings are described below.
(1) The random distribution of coarse aggregates generated variability of about 1.27% in the ordinary concrete’s strength. Du’s group conducted a similar numerical study, and they found that the variability of the concrete’s strength was about 0.91% influences by the spatial location of the coarse aggregates [
58]. That agrees well with the present result. However, it should be pointed out that both of those predictions are significantly lower than the tested strength variability from actual experiments (about 5% to 6%) [
59,
60]. This may because many factors such as initial pores, aggregate shape, and aggregate content influence the compressive strength of concrete, but only the distribution of the coarse aggregates was considered in these simulations.
(2) When the replacement ratios of DCLs were 11%, 22%, 33% and 44%, the corresponding variability in recycled lump concrete’s peak strength was 2.27%, 1.81%, 2.26% and 2.80%, respectively. It can be seen evidently that the strength variability for the recycled lump concrete was generally 0.54–1.53 percentage points greater than that of the ordinary concrete, but this variability had no monotonic trend with an increase in the lump replacement ratio. That higher strength variability of recycled lump concrete may mainly have been due to the coexistence of random distribution of coarse aggregates and random distribution of DCLs. For the same reason, the variability of stress at an axial strain of 2500 or 3000 με was also greater for recycled lump concrete than for ordinary concrete.
In the case where the strength of the new concrete was equal to or less than that of the old concrete, the statistical characteristics of concrete stress are presented in
Table 9. It is easy to observe the following:
(1) When the two strengths were about equal, the variation in the recycled lump concrete’s compressive strength was relatively small as the replacement ratio of DCLs increased. In that situation, the strength variability of the recycled lump concrete was 0.13–0.27 percentage points stronger than that of ordinary concrete. This was mainly because, even though the attributes of the constituents were equal, the distribution of interfaces between new and old concrete grew with increasing amounts of the DCLs.
(2) When the new concrete was distinctly weaker than the old concrete, the compressive strength of recycled lump concrete increased with the amount of DCLs included. The variabilities of the compressive strength and the stress at an axial strain of 2500 με were both slightly higher for recycled lump concrete than for ordinary concrete. However, the variability in the stress at an axial strain of 3000 με for recycled lump concrete intercrossed with that for ordinary concrete.
(3) Overall, the strength variability of recycled lump concrete was smallest when the strengths of the new concrete and old concrete were about equal.
The following formulas were proposed in previous studies [
3,
61] for predicting the compressive strength of recycled lump concrete prisms:
where the subscripts “
RLC”, “
new”, and “
old” denote the recycled lump concrete, the new concrete, and the old concrete, respectively.
Here, Equation (6) was applied to predict the compressive strengths of the 420 recycled lump concrete samples listed in
Table 8. The calculated results agree well with the strength values given by the discrete element method. The relative errors were mostly within the range of −7.56% to 8.51%, and the correlation coefficient
R2 was 0.982. The utility of this empirical formulas applied for predicting the compressive strength of recycled lump concrete prisms was re-confirmed in the current discrete element analyses.
5.4. The Influence of Lump Shape
A real crushing process produces DCLs with various shapes, but scholarly research mainly focused on roughly spherical lumps. The relevance of its conclusions for other shapes remains to be verified. This study took some steps in that direction through applying the discrete element method.
Figure 13 shows the circular and elliptical DCLs in the simulations. These DCLs had the same area, but they had length ratios of their long axis to short axis of 1.0, 1.5, or 2.0 (see
Figure 13a–c). In practice, of course, DCLs with various shapes are used together. The simulations treated four such mixed-use cases (see
Figure 13d–g). During these calculations, the compressive strengths of the new and old concrete were taken as 45.29 MPa and 28.43 MPa, respectively (i.e., the new and old concrete adopted the mesoscale parameters of the concretes with the strengths of 45.29 MPa and 28.43 MPa in
Table 5). The replacement ratio of DCLs was 33% in all cases. The DCLs and coarse aggregates were thrown 30 times in each case of
Figure 13. In this process, the inclination angle of a DCL was not fixed.
Figure 14 displays the simulated compressive stress–stain curves for the seven cases listed in
Figure 13. The mean stress–strain curves and the coefficients of variation for the stresses at axial strains of 2500 and 3000
με are also presented. The plots show the following:
(1) The predicted compressive strength and elastic modulus of recycled lump concrete were only slightly influenced by the shape of the DCLs included. The maximum differences for the strength and modulus were only 5.1% and 5.5%, respectively. This is quite different from the results reported in literature [
62], which observed that the morphology of the coarse aggregates affects the concrete strength and stiffness notably.
(2) The bonding area between new and old concrete was relatively higher for specimens containing DCLs with a larger length ratio of their long axis to short axis. Comparing the seven cases listed in
Figure 13 and
Figure 14, the compressive strength and elastic modulus were higher for Cases A and D in which there was less new concrete–old concrete interfacial area. Cases C and F illustrate the situation with a higher number of interfaces, and the mechanical properties were slightly smaller. As illustrated previously, these relationships arose because the old mortar–new mortar interfaces were relatively weak. When casting recycled lump concrete specimens containing small-perimeter DCLs, the interfacial area bonding the old mortar to the new mortar is minimized, which increases the concrete’s strength and elastic modulus.
A Kolmogorov–Smirnov test [
58] showed that the compressive strengths for the seven cases conform to the normal distribution.
Figure 14h presents the average value
μ, the variance
s of the compressive strength, and the strength
v with their 95% confidence limits. The figure shows that, as the DCL shape changed, the coefficient of variation for the compressive strength fluctuated by 1.89% to 2.46%, but there was no monotonic trend. Therefore, in actual applications of recycled lump concrete, it is feasible to mix in DCLs with a ratio of long axis to short axis smaller than 2.0.