Next Article in Journal
Numerical Prediction of Strength of Socket Welded Pipes Taking into Account Computer Simulated Welding Stresses and Deformations
Next Article in Special Issue
Salt Scaling Resistance of Variable w/c Ratio Air-Entrained Concretes Modified with Polycarboxylates as a Proper Consequence of Air Void System
Previous Article in Journal
High-Entropy Borides under Extreme Environment of Pressures and Temperatures
Previous Article in Special Issue
Effect of Ornamental Stone Waste Incorporation on the Rheology, Hydration, Microstructure, and CO2 Emissions of Ordinary Portland Cement
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Prediction of Ultrasonic Pulse Velocity for Cement, Mortar, and Concrete through a Multiscale Homogenization Approach

College of Civil Engineering and Architecture, Zhejiang University, Hangzhou 310058, China
*
Author to whom correspondence should be addressed.
Materials 2022, 15(9), 3241; https://doi.org/10.3390/ma15093241
Submission received: 6 March 2022 / Revised: 26 April 2022 / Accepted: 28 April 2022 / Published: 30 April 2022
(This article belongs to the Special Issue Multi-Scale Structural Characterization of Cement-Based Composites)

Abstract

:
Ultrasonic testing (UT) is an important method for concrete, and ultrasonic pulse velocity is commonly used to evaluate the quality of concrete materials in existing studies. The ultrasonic pulse velocity of concrete materials is affected by many factors; therefore, it is necessary to establish a quantitative prediction model for the ultrasonic pulse velocity of concrete materials. Based on the multiscale homogenization method, concrete material is divided into different scales of homogenized materials, namely cement paste, mortar, and concrete. Then, a multiscale ultrasonic pulse velocity model is established through a combination of elasticity formulation and the hydration model. At the three scales of cement paste, mortar, and concrete, the elastic parameters and ultrasonic pulse velocity were predicted with the water-to-cement ratio of 0.35, 0.5, and 0.65, respectively. The ultrasonic pulse velocity of concrete with different water-to-cement ratios and different ages were measured in the test and predicted by the model. The results show that the predicted value of ultrasonic pulse velocity is within the error range of ±1.5% of the measured ultrasonic pulse velocity, suggesting that the established prediction model of ultrasonic pulse velocity can reliably predict the velocity change in concrete materials.

1. Introduction

Ultrasonic testing (UT) is a commonly used concrete testing method, with advantages including fast testing speed, convenient operation procedure, and a non-destructive nature. The strength and other parameters of concrete can be evaluated based on its ultrasonic pulse velocity (UPV) [1,2,3]. It is critical to quantify the correlation of the required parameters with the ultrasonic pulse velocity of concrete. Some scholars have obtained empirical formulas by fitting a large amount of data, which are used to determine the strength of concrete relying on ultrasonic pulse velocity [4]. However, empirical correlation formulas cannot be applied to concrete of different conditions and properties. For special types of concrete, such as fiber reinforced concrete [5] and recycled aggregates concrete [6], the correlation between parameters needs to be re-given through experimental data analysis. In addition, many scholars have used ultrasonic pulse velocity to evaluate the damage of concrete [7], and have further completed the imaging of non-uniform damage [8,9]. Consequently, accurate ultrasonic pulse velocity models are needed to predict the UPV of the concrete for quantifying concrete performance. In some studies, the prediction model has been established to predict the strength [10] and porosity [11] of the concrete. Moreover, the forecasting progress can be used to investigate the development of the microstructure in the early hydration process [12].
In order to optimize the evaluation of concrete performance based on ultrasonic pulse velocity detection, it is necessary to study the ultrasonic pulse velocity prediction model of concrete materials. However, the ultrasonic pulse velocity of concrete is affected by various factors, such as the properties of materials [13], mix proportion [14], age, and moisture content [15]. Therefore, the accurate prediction of the ultrasonic pulse velocity of concrete materials remains an important research topic.
The mathematical model for predicting ultrasonic pulse velocity of concrete has been established by analysis of multiple influencing parameters [16]. Furthermore, with the model in the area of data analysis, a satisfactory prediction model can be obtained for ultrasonic pulse velocity [17]. In addition, in the existing studies, based on the elastic theory of homogeneous and isotropic materials, theoretical formulas for the correlation of elastic parameters and ultrasonic pulse velocity can be obtained through theoretical deduction [18]. In order to apply the elasticity formulation for homogeneous isotropic materials on non-homogeneous concrete, the homogenization approach of continuum micromechanics can provide a solution. In previous works, many scholars have studied the micromechanics model to obtain the uniform properties of the material from low-scale to high-scale, such as thermodynamic parameters [19,20], diffusion coefficient [21,22], and freezing behavior [23,24,25]. Many results in continuum micromechanics can be used for elastic parameter prediction of cement-based materials. Constantinides et al. [26] studied the elastic properties of two types of calcium silicate hydrates in cement-based materials through a micro-mechanical model and could predict the macro elasticity of cement paste with high accuracy. Ulm et al. [27] analyzed the applicability of the micromechanical model in concrete materials and studied the poroelasticity of multiscale and multiphase materials. Bernard et al. [28] successfully predicted the elastic modulus of concrete from the elastic parameters of microscopic hydrates by combining micromechanics and hydration models. Related studies on the prediction of elastic properties of cement-based materials through micromechanical models include Hellmich and Mang [29] and Pichler et al. [30]. Moreover, many scholars have studied the early elastic properties of cement-based materials, monitoring the hydration process by nonlinear elastic waves [31,32] and ultrasonic testing [33]. In this paper, the ultrasonic pulse velocity model of concrete materials was established by studying changes of elastic parameter in multiphase concrete materials.
In this study, based on the multiscale homogenization method and the elasticity formulation of homogenized multiphase material, a multiscale ultrasonic pulse velocity model for multiphase concrete material was established. Combined with the prediction of hydration process, the elastic modulus and ultrasonic pulse velocity could be predicted by the model in three-scale homogenized materials of cement paste, mortar, and concrete. The accuracy of the prediction model was verified by comparing the measured ultrasonic pulse velocity with the predicted value of the model. Based on this prediction model, the change of ultrasonic pulse velocity during the hydration process of concrete can be calibrated for future studies. Moreover, the concrete damage parameters can be introduced based on the established model for calculating the ultrasonic pulse velocity of concrete, which provides a theoretical basis for the study of ultrasonic pulse velocity of concrete damage.

2. Micromechanical Homogenization

Based on the elastic theory of homogeneous and isotropic materials, the formula for the correlation of elastic parameters such as ultrasonic pulse velocity and elastic modulus can be obtained through the derivation of the theoretical formula. The ultrasonic wavelength used in concrete ultrasonic testing is larger than the size of the discontinuities existing inside concrete. Therefore, the ultrasonic waves do not interact with these inhomogeneities, and concrete can reasonably be regarded as a homogeneous material [18]. The correlation theoretical formula can be applied in the ultrasonic testing of concrete:
V = E ρ 1 ν ( 1 + ν ) ( 1 2 ν ) = K + 4 / 3 G ρ
where V is the ultrasonic pulse velocity of the homogeneous material; E is the elastic modulus of the material; ρ is the density of the material; ν is the Poisson’s ratio of the material; and K and G are the bulk and shear modulus of the material, respectively.
However, concrete is a composite of different components, the overall properties of which cannot be directly determined when viewed as a homogeneous material. Therefore, we need to determine the K, G, and ρ of the concrete first.
The focus is then on how to predict K, G and ρ, or the effective properties of multiscale multiphase concrete composite (Khom, Ghom, ρhom). In this paper, the iterative homogenization method is used in concrete such a multiphase material to predict the parameters of K and G at different scales.
The homogenization method describes the overall elastic properties of representative volume elements (RVE) at different scales. The size of the representative volume unit defined within it needs to be much smaller than the overall analysis, and the characteristic length of the inhomogeneity in the unit needs to be much smaller than the size of the unit. In this scheme, the material in the RVE can be considered as uniform at the statistical level [34]. The elastic properties of higher-scale RVEs are described by the uniform structure at the micro scale. By repeatedly performing the homogenization scheme, the overall elastic properties of multiphase inhomogeneous materials can be predicted. For concrete, the multiscale division and homogenization scheme is shown in Figure 1.
According to the homogenization method, RVE selected at different scales can be regarded as uniform, so that the elastic parameters of cement-based materials at each scale are calculated according to the theory of micromechanics. Since each component is considered isotropic, the homogenized material is isotropic as well.
Combined with research results of micromechanics on the elasticity formulation of homogeneous multiphase materials, the authors in [35] summarize the following formulations for bulk modulus and shear modulus of homogeneous multiphase materials:
{ K hom = r = 1 n f r K r P r / r = 1 n f r P r G hom = r = 1 n f r G r Q r / r = 1 n f r Q r ρ hom = r = 1 n f r ρ r
where K hom , G hom , ρ hom , K r , G r , and ρ r are the bulk modules, shear modulus, density of homogenized material, and phase r, respectively; P r , Q r , and f r are the compressibility, shear compliance, and volume fraction of each phase r, respectively; and n is the number of phases in the composite.
The elastic modulus and Poisson’s ratio can be calculated according to the following Equation (3):
{ E hom = 9 K hom G hom 3 K hom + G hom v hom = 3 K hom 2 G hom 6 K hom + 2 G hom
where E hom and ν hom are the elastic modulus and Poisson’s ratio of homogenized material.
The compressibility P and shear compliance Q are determined by the shape and elastic properties of each phase in the material. In this study, it is considered that each component exists in a nearly spherical state in undamaged concrete. According to Berryman [36], the compressibility P and shear compliance Q are calculated considering spherical inclusions, as seen in Equation (4):
{ P r s p h = K 0 + 4 G 0 / 3 K r + 4 G 0 / 3 Q r s p h = G 0 + F 0 G r + F 0
where K 0 and G 0 are the compressibility and shear compliance of the matrix phase, which represents the lower scale medium. The F 0 can be calculated as:
F 0 = ( G 0 / 6 ) [ ( 9 K 0 + 8 G 0 ) / ( K 0 + 2 G 0 ) ]
According to Equations (2)–(5), the compressibility and shear compliance of homogeneous media at different scales can be obtained, depending on a suitable iteration scheme.
Three scales are considered in concrete materials: (1) cement paste composed of CSH skeleton and pore substances; (2) mortar composed of cement paste and fine aggregates; and (3) concrete composed of mortar and coarse aggregates. The Mori-Tanaka scheme [37] and self-consistent scheme [38] were selected for iterative calculation of the homogenization process of media at different scales. These two approaches are commonly used to predict the overall elastic properties of homogeneous and isotropic composite materials from the micromechanics point of view. There are some connections between them; however, different homogenization processes need corresponding schemes to evaluate the effective properties of composites.
Usually, the Mori-Tanaka scheme defines a medium as reference medium, which has a large volume fraction in the mixture. The self-consistent scheme takes each medium as the reference medium for step-by-step iterative calculation, which is suitable for a closer volume fraction of each phase. Significantly, it leads to different values being substituted in Equation (4) when different schemes are used. K 0 and G 0 are equal to K hom and G hom of the RVE at the lower scale when using the Mori-Tanaka scheme. In the self-consistent scheme, K 0 and G 0 refer to K hom and G hom of the overall medium at the previous iteration, which has multiple iterations at the same scale.
Thus, the self-consistent scheme is used for cement paste due to its complex components, including more overlaps among the components and more pores. Obviously, using the Mori-Tanaka scheme for mortar is more suitable because the volume fraction of cement paste is so large that it can be considered as the main body, and the fine aggregate as spherical inclusions. In the same way, the mortar can be considered as the main body and the coarse aggregate as spherical inclusions in the concrete level. We chose the Mori-Tanaka scheme for concrete. The prediction of elastic parameter and ultrasonic pulse velocity in this paper refers to the above selection of the iterative scheme.

3. Ultrasonic Pulse Velocity Model for Cement Paste

3.1. Hydration Model and Volume Fraction

The RVE of cement paste consists of unhydrated cement clinker, water, hydrates, and air voids [28,29,30,39]. From the hydration degree of 0 (initial) to the maximum (completed), the volume fraction of each material constantly evolves. Following the Powers-Acker model [40], the volume fraction of clinker (clin), water (H2O), hydrates (hyd), and air voids as functions of the hydration degree ξ and the water-to-cement ratio w/c are:
f clin ( ξ ) = 1 ξ 1 + ρ clin ρ H 2 O ( w / c ) = 20 ( 1 ξ ) 20 + 63 ( w / c ) 0
f H 2 O ( ξ ) = ρ clin [ ( w / c ) 0.42 ξ ] ρ H 2 O [ 1 + ρ clin ρ H 2 O ( w / c ) ] = 63 [ ( w / c ) 0.42 ξ ] 20 + 63 ( w / c ) 0
f hyd ( ξ ) = 1.42 ρ clin ξ ρ hyd [ 1 + ρ clin ρ H 2 O ( w / c ) ] = 43.15 ξ 20 + 63 ( w / c )
f air ( ξ ) = ( 1 + 0.42 ρ clin ρ H 2 O 1.42 ρ clin ρ hyd ) ξ 1 + ρ clin ρ H 2 O ( w / c ) = 3.31 ξ 20 + 63 ( w / c )
where ξ is the hydration degree; w/c is the water-to-cement ratio; f clin , f H 2 O , f hyd , and f air refer to the volume fraction of unhydrated cement clinker, water, hydrates, and air voids, respectively; and ρ clin , ρ H 2 O , and ρ hyd refer to the mass densities of clinker, water, and hydrates, respectively, referencing the parameters in Table 1.
The maximum hydration degree is affected by the water-to-cement ratio. If the water-to-cement ratio is below 0.42, the cement clinker cannot be fully hydrated. The hydration degree reaches the maximum when the water is exhausted. If the water-to-cement ratio is larger than 0.42, the maximum hydration degree is equal to 1 when the cement clinker is fully consumed [30,39]. The maximum hydration degree follows as:
max ξ = { ( w / c ) 0.42 for ( w / c ) 0.42 f clin ( max ξ ) 0 , f H 2 O ( max ξ ) = 0 , 1 for ( w / c ) 0.42 f clin ( max ξ ) = 0 , f H 2 O ( max ξ ) 0 .
where ξ is the hydration degree and w/c is the water-to-cement ratio. The hydration degree evolves with age, influenced by the water-to-cement ratio, aggregate properties, and so on. In order to evaluate the proportion of each component in the composites, we need a scheme to predict the hydration degree.
With the calculation of pore distribution and hydration heat by the DuCOM model [42], the increase of the hydration degree ξ with a specific water-to-cement ratio can be quantified. Different hydration processes with water-to-cement ratios of 0.35, 0.5, and 0.65 are shown in Figure 2.
Obviously, the hydration degree shows a notable acceleration since initiation, and a slight increase in the hydration degree was found in the stable period. Based on the analysis of the predicted results, the mixture with a larger water-to-cement ratio has a higher hydration degree at the same age before the process is ended.

3.2. Velocity Prediction

At the scale of cement paste (Level II), the prediction of ultrasonic pulse velocity requires correlation parameters of clinker, water, hydrates, and air voids to predict the elastic parameter, as shown in Table 1.
According to Equations (2)–(5), bulk modulus K, shear modulus G, elastic modulus E, and Poisson’s ratio υ of homogenized material and phase were calculated by the model with the self-consistent scheme. The correlation curves between different parameters and the hydration degree are shown in Figure 3.
At the same degree of hydration, a higher water-to-cement ratio shows a lower bulk modulus, shear modulus, and elastic modulus. With the growing hydration degree, the bulk modulus, shear modulus, and elastic modulus of the cement paste all show an increasing trend, and the differences between these parameters also increase gradually. For cement paste with a high water-to-cement ratio in the initial stage of hydration, the shear modulus and elastic modulus increase from around 0, which is also close to the corresponding parameters of water. Poisson’s ratio decreases with the increase of the hydration degree, and Poisson’s ratios of cement paste with different water-to-cement ratios approach 0.25 when the hydration is completed.
The ultrasonic pulse velocity (UPV) of the cement paste can be calculated using Equation (1) with the elastic parameters. The UPV of cement paste as a function of the hydration degree is shown in Figure 4.
The predicted ultrasonic pulse velocity of the cement paste also increases with the degree of hydration, displaying the same trend as the bulk modulus, shear modulus, and elastic modulus. In the case of a high water-to-cement ratio, the initial value of the wave speed is about 1.5 km/s, which is close to the ultrasonic speed in water. When the hydration is completed, the influence of the water-to-cement ratio on the final UPV is generally linear with a negative correlation. The influence on the initial UPV is more sensitive when the water-to-cement ratio is less than 0.4. Once the water-to-cement ratio is greater than 0.4, its increase has little effect on the initial UPV. The main influencing factor on the final UPV is the proportion of hydrates in the cement paste. When the volume fraction of water, which mainly affects the initial UPV, is larger in the mixture, the ultrasonic pulse velocity will continue to approach the ultrasonic pulse velocity in water.
Considering different schemes, the predicted results are shown in Figure 5. The ultrasonic pulse velocity obtained by the Mori-Tanaka scheme increases linearly with an increasing degree of hydration, starting from about 3.2 km/s. The velocity calculated by the self-consistent scheme gradually increases from about 1.5 km/s, and this value is close to the value of the ultrasonic pulse velocity in water. The cement paste includes hydration products, water, air voids, and unhydrated clinker. Before the hydration starts, the theoretical value of the wave speed is supposed to be close to the wave speed in water when the water-to-cement ratio is relatively large. Meanwhile, the internal components of the cement paste are complex, with impurities and multitudes of and pores, meaning that there is no component with a large volume ratio in the early stage of hydration. Therefore, the self-consistent scheme is adopted for calculation.

4. Ultrasonic Pulse Velocity Model Upgraded to Mortar and Concrete

4.1. Velocity Model for Mortar

When the ultrasonic pulse velocity model upgrades to mortar and concrete, the volume fraction of the aggregate in each scale material should be determined. According to the mix proportion and density of the aggregate, the volume fraction can be calculated using Equations (11) and (12).
{ f s a n d   = M s a n d   / ρ s a n d   1 M g r a v   / ρ g r a v   f c e m   = 1 f s a n d  
{ f g r a v = M g r a v ρ g r a v f m o t   = 1 f g r a v  
where f s a n d   and f c e m   are the volume fraction of fine aggregate and cement paste in the mortar, respectively; f g r a v and f m o t   are the volume fraction of coarse aggregate and mortar in the concrete, respectively; M s a n d   and M g r a v   are the amount of the fine aggregate and coarse aggregate in mix proportion, respectively; and ρ s a n d   and ρ g r a v   are the density of the fine and coarse aggregate, respectively.
The quantified prediction of the mortar and concrete also requires the elastic parameters of the cement paste, while the required physical parameters of the aggregate are shown in Table 2.
The parameter can be given for the prediction of ultrasonic velocity, which is calculated according to the commonly used mix proportion, as shown in Table 3.
The calculated results indicate that the whole cement paste is taken as one phase in the two-phase material composed of cement paste and fine aggregate. According to the established prediction model, the curves of the elastic modulus and Poisson’s ratio of mortar are predicted with the Mori-Tanaka scheme, as shown in Figure 6. Then, the predicted UPVs of mortar with different w/c ratios, which are shown in Figure 7, can be easily obtained using Equation (1).
The overall trend of the predicted curve is consistent with the predicted results of the cement paste. Compared with cement paste, the elastic modulus and ultrasonic pulse velocity of mortar are larger, while the Poisson’s ratio is smaller at the same hydration degree.

4.2. Velocity Model for Concrete

According to the preceding homogenization process of mortar scales, the mortar is considered one phase to the composite, considering the concrete as a two-phase material consisting of mortar and coarse aggregate. With the calculating parameters of the concrete scale shown in Table 2 and Table 3, the predicted results can be obtained by the preceding calculated process with the Mori-Tanaka scheme.
The predicted results for the elastic modulus, Poisson’s ratio, and UPV of concrete as functions of the hydration degree are shown in Figure 8 and Figure 9.
The overall trend of the predicted curves is consistent on the three scales of cement paste, mortar, and concrete. However, the aggregate, occupying the bulk of the volume, has stable properties and basically unchanged relevant parameters with the hydration process. This means that with the progress of hydration, the differences between the prediction values of the elastic modulus, Poisson’s ratio, and ultrasonic pulse velocity corresponding to different water-to-cement ratios are decreasing. Due to the addition of the aggregate phase, the ultrasonic pulse velocity of the homogenized materials at the scale of mortar and concrete increases at the same degree of hydration, compared with the scale of cement paste. In the case of the same water-to-cement ratio from the initial stage of hydration to completed hydration, the change of ultrasonic pulse velocity is also smaller than that of the cement paste scale.

5. Experimental Verification

5.1. Experimental Setup

To verify the accuracy of the multiscale ultrasonic pulse velocity prediction model, a validation experiment was designed.
Three series of concrete specimens were prepared with two different target ages (7 and 28 days) and three different w/c ratios (0.37, 0.42, and 0.55). Ordinary Portland cement with a grade of P.O 32.5 was used. The fine aggregate used for the concrete was local river sand with a fineness modulus of 2.3 and a relative density of 2.65. The coarse aggregate used for concrete was local limestone, with sizes of 5–20 mm and a relative density of 2.62. The hydration degree of the three concrete mixtures at the ages of 7 days and 28 days were predicted by the DuCOM model, which is shown in Table 4, as well as the volume fraction of components at each age.
The concrete specimens with dimensions of 100 mm × 100 mm × 100 mm were demolded after 48 h. Each sample was placed in water to cure at curing room temperature (i.e., 20 ± 2 °C) for 7 days or 28 days.
The specimens and UPV testing tools are shown in Figure 10; the UPV of the specimen from three concrete types were measured at the age of 7 days and 28 days.
The “ZBL-U5200” device was adopted to measure the UPV of the concrete specimen. The measurement of UPV was carried out immediately after the specimen left the curing condition to ensure the stability of its moisture content and other states. The transducers were placed on two parallel surfaces of the specimen, with the Vaseline coupling between the transducers and the specimen. Five measuring points were designed on two parallel surfaces, with three instances of measurements on each measuring point. The ultrasonic pulse velocity detection test adopts the direct transmission, as shown in Figure 11.

5.2. Experimental Results and Comparison

The testing results of the UPV are shown in Figure 12.
With the increase of age, the ultrasonic pulse velocities of groups A, B, and C showed an overall increasing trend. At the same age, the ultrasonic pulse velocity of the specimens in the low water-to-cement ratio group was significantly higher than that in the high water-to-cement ratio group. The error bar shows that the standard error of the measured value of each specimen in the test was also controlled within a reasonable range.
According to the experimental material parameters, the elastic modulus and UPV of concrete can be predicted based on the multiscale ultrasonic pulse velocity prediction model, as shown in Figure 13.
It can be observed that the trend of ultrasonic pulse velocity is consistent with that of the hydration degree. In the early stage of hydration, due to severe hydration reaction, the ultrasonic pulse velocity and elastic modulus increase rapidly, reaching about 80% of the maximum value within 5 days of age. The development trend then slows down. At 28 days of age, the speed is extremely slow or even remains unchanged.
For concrete at the same age, the ultrasonic velocity of those with a water-to-cement ratio of 0.55 is the smallest, which develops to about 4300 m/s at 28 d. The ultrasonic velocity of concrete with a water-to-cement ratio of 0.37 is the highest, which finally develops to about 4500 m/s. Although concrete with a high water-to-cement ratio has a higher degree of hydration at the same age, its elastic model and ultrasonic pulse velocity values are still lower than those of low water-to-cement ratio concrete. Therefore, the mix design has a great influence on the elastic parameters and ultrasonic pulse velocity of the concrete material.
Comparing the predicted results with the measured results, as shown in Figure 14, it can be found that the result points are close to the equal line. The numerical simulation is in good agreement with the experiment, indicating the accuracy of the multiscale ultrasonic pulse velocity model.
The standard deviation of the measured values of each group is marked in Figure 14, showing that the dispersion of the experiment results remains within a reasonable range. The error of the predicted UPV relative to the experimental measured UPV is basically within ±1%. In the group of 28 days age, the relative error with a water-to-cement ratio of 0.55 is slightly larger, but is still within the range of ±1.5%.
The predicted ultrasonic pulse velocity values highly fit with the measured values, showing a satisfactory agreement. If the required parameters are available, the elastic parameters and UPV of the multiphase cement-based material can be predicted by this multiscale ultrasonic pulse velocity model.
With more experimental data, further research will be able to verify the applicability of this new calculation model in other cases.
According to Ye’s study [12], the ultrasonic pulse velocity of concrete was calculated by the multiscale theoretical calculation model in this paper with the experimental parameters. Figure 15 shows the results obtained by the theoretical model with different w/c ratios in Ye’s study, compared with the measured results.
The results indicate that the predicted UPV of concrete is close to the measured UPV, especially after reaching a high hydration degree. It is thought that the air entrapped in the fresh cement paste is the reason for random error in the initial stage, which is more pronounced in the case of high w/c ratios. In addition, the effect of the different cement materials’ properties provided in the paper also leads to some random errors. In general, the multiscale ultrasonic pulse velocity model exhibits a good predictive capacity when it is adapted to predict the ultrasonic pulse velocity of concrete.

6. Conclusions

In this paper, a multiscale ultrasonic pulse velocity model for multiphase concrete materials is established based on the homogenization approach combined with the model of the hydration process. During the hydration process, the elastic parameters and ultrasonic pulse velocity of cement paste, mortar, and concrete can be predicted by the model. The accuracy of this model is verified by concrete ultrasonic testing and other research data.
  • Based on the multiscale homogenization method and the elasticity formulation of homogenized multiphase materials, a multiscale ultrasonic pulse velocity model is established which can predict elastic parameters and ultrasonic pulse velocity during the hydration process, according to material parameters, mixture, and age.
  • In this model, the iterative calculation of elastic parameters requires different schemes at different scales. The self-consistent scheme is applied at the scale of cement paste, while the Mori-Tanaka scheme is used at the mortar and concrete scales.
  • The elastic parameters and ultrasonic pulse velocity at three scales of cement paste, mortar, and concrete were predicted by this model with the water-to-cement ratios of 0.35, 0.5, and 0.65, respectively. At the scale of cement paste, the volume fraction of water and the proportion of hydrates are the main influencing factors of the initial ultrasonic pulse velocity and final ultrasonic pulse velocity, respectively. At the scale of mortar and concrete, the aggregates make the influence of the water-to-cement ratio gradually decrease as hydration progresses, with the large volume fraction and stable nature.
  • According to the experiments detailed in this paper, the relative error of the measured ultrasonic pulse velocity and the predicted value is within ±1.5%. According to Ye’s study [12], the ultrasonic pulse velocity of concrete was predicted by the prediction model, and the obtained prediction value agreed well with the measured value. Thus, the model is reliable for predicting the ultrasonic pulse velocity of concrete materials. This has reference value for the prediction of cement-based ultrasonic pulse velocity. However, for special concrete, the elastic parameters of special aggregates are very different from that of ordinary concrete, which cannot be ignored in the prediction of elastic parameters. Research will be carried out on the prediction of ultrasonic pulse velocity for special concrete in the future.

Author Contributions

Validation, formal analysis, data curation, visualization, J.J.; conceptualization, methodology, investigation, and software, F.G. and J.J.; writing—original draft preparation, review and editing, J.J., F.G. and D.Z. (Dian Zhi); project administration, F.G.; resources, supervision, and funding acquisition, F.G. and D.Z. (Dawei Zhang). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Key-area Research and Development Program of Guangdong Province (grant numbers: 2019B111107002), the Natural Science Foundation of Zhejiang Province (grant numbers: LQ20E080016), and the National Natural Science Foundation of China (grant numbers: 52008367, 51820105012).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article. The data presented in this study are available in [12] and this paper.

Acknowledgments

The authors acknowledge the efforts of reviewers to improve the quality of this paper.

Conflicts of Interest

The authors declare no conflict of interest for this study.

References

  1. Lin, Y.; Shih-Fang, K.; Hsiao, C.; Chao-Peng, L. Investigation of Pulse Velocity-Strength Relationship of Hardened Concrete. ACI Mater. J. 2007, 104, 344. [Google Scholar]
  2. Trtnik, G.; Kavčič, F.; Turk, G. Prediction of Concrete Strength Using Ultrasonic Pulse Velocity and Artificial Neural Networks. Ultrasonics 2009, 49, 53–60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Lin, Y.; Lin, Y.; Chan, C. Use of Ultrasonic Pulse Velocity to Estimate Strength of Concrete at Various Ages. Mag. Concr. Res. 2016, 68, 739–749. [Google Scholar] [CrossRef]
  4. Phoon, K.K.; Wee, T.H.; Loi, C.S. Development of Statistical Quality Assurance Criterion for Concrete Using Ultrasonic Pulse Velocity Method. Mater. J. 1999, 96, 568–574. [Google Scholar]
  5. Hedjazi, S.; Castillo, D. Relationships among Compressive Strength and UPV of Concrete Reinforced with Different Types of Fibers. Heliyon 2020, 6, e03646. [Google Scholar] [CrossRef] [PubMed]
  6. Velay-Lizancos, M.; Martinez-Lage, I.; Azenha, M.; Vázquez-Burgo, P. Influence of Temperature in the Evolution of Compressive Strength and in its Correlations with UPV in Eco-Concretes with Recycled Materials. Constr. Build. Mater. 2016, 124, 276–286. [Google Scholar] [CrossRef]
  7. Yan, D.; Tian, Y.; Liu, K.; Chen, S.; Zeng, Q.; Ruan, S. Evaluation of Mechanical Properties of Concrete after Exposure to Elevated Temperatures Using Ultrasonic Pulse Velocity Measurement and a Split-Hopkinson Pressure Bar. J. Mater. Civil Eng. 2021, 33, 04021352. [Google Scholar] [CrossRef]
  8. Wang, Z.; Wang, Y.; Zhang, D.; Ueda, T. Experimental CT Imaging Verification for Non-Uniform Frost Damage Simulation of Concrete in Mesoscale. In Proceedings of the 7th International Conference of Asian Concrete Federation, Hanoi, Vietnam, 30 October–2 November 2016. [Google Scholar]
  9. Wang, Z.; Gong, F.; Zhang, D.; Hayashida, H.; Ueda, T. Mesoscale Simulation of Concrete Behavior with Non-Uniform Frost Damage with Verification by CT Imaging. Constr. Build Mater. 2017, 157, 203–213. [Google Scholar] [CrossRef]
  10. Ashrafian, A.; Taheri Amiri, M.J.; Rezaie-Balf, M.; Ozbakkaloglu, T.; Lotfi-Omran, O. Prediction of Compressive Strength and Ultrasonic Pulse Velocity of Fiber Reinforced Concrete Incorporating Nano Silica Using Heuristic Regression Methods. Constr. Build Mater. 2018, 190, 479–494. [Google Scholar] [CrossRef]
  11. Maalej, S.; Lafhaj, Z.; Bouassida, M. Micromechanical Modelling of Dry and Saturated Cement Paste: Porosity Assessment Using Ultrasonic Waves. Mech. Res. Commun. 2013, 51, 8–14. [Google Scholar] [CrossRef]
  12. Ye, G.; van Breugel, K.; Fraaij, A.L.A. Experimental Study and Numerical Simulation on the Formation of Microstructure in Cementitious Materials at Early Age. Cem. Concr. Res. 2003, 33, 233–239. [Google Scholar] [CrossRef]
  13. Robeyst, N.; Gruyaert, E.; Grosse, C.U.; De Belie, N. Monitoring the Setting of Concrete Containing Blast-Furnace Slag by Measuring the Ultrasonic P-Wave Velocity. Cem. Concr. Res. 2008, 38, 1169–1176. [Google Scholar] [CrossRef]
  14. Abo-Qudais, S.A. Effect of Concrete Mixing Parameters on Propagation of Ultrasonic Waves. Constr. Build Mater. 2005, 19, 257–263. [Google Scholar] [CrossRef]
  15. Lafhaj, Z.; Goueygou, M.; Djerbi, A.; Kaczmarek, M. Correlation between Porosity, Permeability and Ultrasonic Parameters of Mortar with Variable Water/Cement Ratio and Water Content. Cem. Concr. Res. 2006, 36, 625–633. [Google Scholar] [CrossRef]
  16. Lin, Y.; Lai, C.P.; Yen, T. Prediction of Ultrasonic Pulse Velocity (UPV) in Concrete. ACI Mater. J. 2003, 100, 21–28. [Google Scholar]
  17. Ayaz, Y.; Kocamaz, A.F.; Karakoç, M.B. Modeling of Compressive Strength and UPV of High-Volume Mineral-Admixtured Concrete Using Rule-Based M5 Rule and Tree Model M5P Classifiers. Constr. Build Mater. 2015, 94, 235–240. [Google Scholar] [CrossRef]
  18. Chang, T.; Lin, H.; Chang, W.; Hsiao, J. Engineering Properties of Lightweight Aggregate Concrete Assessed by Stress Wave Propagation Methods. Cem. Concr. Compos. 2006, 28, 57–68. [Google Scholar] [CrossRef]
  19. Liu, J.; Xu, S.; Zeng, Q. A Multi-Scale Micromechanical Investigation on Thermal Conductivity of Cement-Based Composites. In Proceedings of the IOP Conference Series: Materials Science and Engineering, Sanya, China, 19–21 November 2016. [Google Scholar]
  20. Niu, Z.; Shen, J.; Wang, L.; Yang, R. Thermo-Poroelastic Modelling of Cement Sheath: Pore Pressure Response, Thermal Effect and Thermo-Osmotic Effect. Eur. J. Environ. Civ. Eng. 2019, 26, 657–682. [Google Scholar] [CrossRef]
  21. Yang, R.; Lemarchand, E.; Fen-Chong, T. A Micromechanics Model for Solute Diffusion Coefficient in Unsaturated Granular Materials. Transp. Porous Med. 2016, 111, 347–368. [Google Scholar] [CrossRef]
  22. Yang, R.; Lemarchand, E.; Fen-Chong, T.; Li, K. Prediction of Permeability of Monodisperse Granular Materials with a Micromechanics Approach. J. Appl. Geophys. 2016, 127, 82–90. [Google Scholar] [CrossRef]
  23. Wang, Y.; Li, J.; Ueda, T.; Zhang, D.; Deng, J. Meso-Scale Mechanical Deterioration of Mortar Subjected to Freeze Thaw Cycles and Sodium Chloride Attack. Cem. Concr. Compos. 2021, 117, 103906. [Google Scholar] [CrossRef]
  24. Yang, R.; Lemarchand, E.; Fen-Chong, T.; Azouni, A. A Micromechanics Model for Partial Freezing in Porous Media. Int. J. Solids Struct. 2015, 75, 109–121. [Google Scholar] [CrossRef]
  25. Zeng, Q.; Fen-Chong, T.; Dangla, P.; Li, K. A Study of Freezing Behavior of Cementitious Materials by Poromechanical Approach. Int. J. Solids Struct. 2011, 48, 3267–3273. [Google Scholar] [CrossRef]
  26. Constantinides, G.; Ulm, F. The Effect of Two Types of C-S-H on the Elasticity of Cement-Based Materials: Results from Nanoindentation and Micromechanical Modeling. Cem. Concr. Res. 2004, 34, 67–80. [Google Scholar] [CrossRef]
  27. Ulm, F.J.; Constantinides, G.; Heukamp, F.H. Is Concrete a Poromechanics Materials?—A Multiscale Investigation of Poroelastic Properties. Mater. Struct. 2004, 37, 43–58. [Google Scholar] [CrossRef]
  28. Bernard, O.; Ulm, F.; Lemarchand, E. A Multiscale Micromechanics-Hydration Model for the Early-Age Elastic Properties of Cement-Based Materials. Cem. Concr. Res. 2003, 33, 1293–1309. [Google Scholar] [CrossRef]
  29. Hellmich, C.; Mang, H. Shotcrete Elasticity Revisited in the Framework of Continuum Micromechanics: From Submicron to Meter Level. J. Mater. Civ. Eng. 2005, 17, 246–256. [Google Scholar] [CrossRef]
  30. Pichler, B.; Hellmich, C.; Eberhardsteiner, J. Spherical and Acicular Representation of Hydrates in a Micromechanical Model for Cement Paste: Prediction of Early-Age Elasticity and Strength. Acta Mech. 2009, 203, 137–162. [Google Scholar] [CrossRef]
  31. Tragazikis, I.K.; Kordatou, T.Z.; Exarchos, D.A.; Dalla, P.T.; Matikas, T.E. Monitoring the Hydration Process in Carbon Nanotube Reinforced Cement-Based Composites Using Nonlinear Elastic Waves. Appl. Sci. 2021, 11, 1720. [Google Scholar] [CrossRef]
  32. Lefever, G.; Snoeck, D.; De Belie, N.; Van Hemelrijck, D.; Aggelis, D.G. Elastic Wave Monitoring of Cementitious Mixtures Including Internal Curing Mechanisms. Sensors 2021, 21, 2463. [Google Scholar] [CrossRef]
  33. Iliopoulos, S.; Vila Neves, R.; Chenu, A.; Angelis, D. Early Age Ultrasonic Testing of Concrete Composed of Glass Spheres as Aggregates. In Proceedings of the Structural Faults and Repair 2016 International Conference Proceedings, Edinburgh, Scotland, 17–19 May 2016. [Google Scholar]
  34. Gong, F.; Ueda, T.; Wang, Y.; Zhang, D.; Wang, Z. Mesoscale Simulation of Fatigue Behavior of Concrete Materials Damaged by Freeze-Thaw Cycles. Constr. Build Mater. 2017, 144, 702–716. [Google Scholar] [CrossRef] [Green Version]
  35. Gong, F.; Wang, Y.; Ueda, T.; Zhang, D. Modeling and Mesoscale Simulation of Ice-Strengthened Mechanical Properties of Concrete at Low Temperatures. J. Eng. Mech. 2017, 143, 04017022. [Google Scholar] [CrossRef]
  36. Berryman, J.G. Long-Wavelength Propagation in Composite Elastic Media II. Ellipsoidal Inclusions. J. Acoust. Soc. Am. 1980, 68, 1820–1831. [Google Scholar] [CrossRef]
  37. Mori, T.; Tanaka, K. Average Stress in Matrix and Average Elastic Energy of Materials with Misfitting Inclusions. Acta Metall. 1973, 21, 571–574. [Google Scholar] [CrossRef]
  38. Kröner, E. Bounds for Effective Elastic Moduli of Disordered Materials. J. Mech. Phys. Solids 1977, 25, 137–155. [Google Scholar] [CrossRef]
  39. Ullah, S.; Pichler, B.; Scheiner, S.; Hellmich, C. Influence of Shotcrete Composition on Load-Level Estimation in NATM-tunnel Shells: Micromechanics-Based Sensitivity Analyses. Int. J. Numer. Anal. Methods Géoméch. 2012, 36, 1151–1180. [Google Scholar] [CrossRef]
  40. Acker, P. Micromechanical Analysis of Creep and Shrinkage Mechanisms. In Creep, Shrinkage and Durability Mechanics of Concrete and Other Quasi-Brittle Materials; Elsevier: Cambridge, MA, USA, 2001; pp. 15–25. [Google Scholar]
  41. Tennis, P.D.; Jennings, H.M. A Model for Two Types of Calcium Silicate Hydrate in the Microstructure of Portland Cement Pastes. Cem. Concr. Res. 2000, 30, 855–863. [Google Scholar] [CrossRef]
  42. Maekawa, K. Multi-Scale Modeling of Structural Concrete; CRC Press: Boca Raton, FL, USA, 2008. [Google Scholar]
Figure 1. Multiscale homogenization scheme for concrete materials.
Figure 1. Multiscale homogenization scheme for concrete materials.
Materials 15 03241 g001
Figure 2. The predicted curve of the hydration degree increasing with age from the DuCOM model.
Figure 2. The predicted curve of the hydration degree increasing with age from the DuCOM model.
Materials 15 03241 g002
Figure 3. The prediction of elastic parameters based on the homogenization process: (a) bulk modulus; (b) sheer modulus; (c) elastic modulus; and (d) Poisson’s ratio.
Figure 3. The prediction of elastic parameters based on the homogenization process: (a) bulk modulus; (b) sheer modulus; (c) elastic modulus; and (d) Poisson’s ratio.
Materials 15 03241 g003
Figure 4. The ultrasonic pulse velocity prediction of cement paste: (a) hydration process and (b) the influence of the water-to-cement ratio.
Figure 4. The ultrasonic pulse velocity prediction of cement paste: (a) hydration process and (b) the influence of the water-to-cement ratio.
Materials 15 03241 g004
Figure 5. The ultrasonic pulse velocity prediction based on the different schemes.
Figure 5. The ultrasonic pulse velocity prediction based on the different schemes.
Materials 15 03241 g005
Figure 6. The predicted elastic parameters of mortar as a function of hydration degree: (a) elastic modulus and (b) Poisson’s ratio.
Figure 6. The predicted elastic parameters of mortar as a function of hydration degree: (a) elastic modulus and (b) Poisson’s ratio.
Materials 15 03241 g006
Figure 7. The ultrasonic pulse velocity prediction of mortar.
Figure 7. The ultrasonic pulse velocity prediction of mortar.
Materials 15 03241 g007
Figure 8. The predicted elastic parameters of concrete as functions of the hydration degree: (a) elastic modulus and (b) Poisson’s ratio.
Figure 8. The predicted elastic parameters of concrete as functions of the hydration degree: (a) elastic modulus and (b) Poisson’s ratio.
Materials 15 03241 g008
Figure 9. The ultrasonic pulse velocity prediction of concrete.
Figure 9. The ultrasonic pulse velocity prediction of concrete.
Materials 15 03241 g009
Figure 10. The measurement of ultrasonic pulse velocity of concrete specimen.
Figure 10. The measurement of ultrasonic pulse velocity of concrete specimen.
Materials 15 03241 g010
Figure 11. The measurement of ultrasonic pulse velocity: (a) layout of measuring points and (b) direct transmission.
Figure 11. The measurement of ultrasonic pulse velocity: (a) layout of measuring points and (b) direct transmission.
Materials 15 03241 g011
Figure 12. The result of ultrasonic testing with standard error.
Figure 12. The result of ultrasonic testing with standard error.
Materials 15 03241 g012
Figure 13. The predicted curves at different ages based on the multiscale ultrasonic pulse velocity prediction model: (a) elastic modulus and (b) ultrasonic pulse velocity.
Figure 13. The predicted curves at different ages based on the multiscale ultrasonic pulse velocity prediction model: (a) elastic modulus and (b) ultrasonic pulse velocity.
Materials 15 03241 g013
Figure 14. The comparison diagram between the predicted value of the multiscale UPV model and the measured value.
Figure 14. The comparison diagram between the predicted value of the multiscale UPV model and the measured value.
Materials 15 03241 g014
Figure 15. The comparison between predicted UPV and measured UPV of the hydration process with different w/c ratios.
Figure 15. The comparison between predicted UPV and measured UPV of the hydration process with different w/c ratios.
Materials 15 03241 g015
Table 1. Physical parameters of different phases for cement paste.
Table 1. Physical parameters of different phases for cement paste.
PhaseBulk Modulus K
(GPa)
Shear Modulus G
(GPa)
Density ρ ( kg/m 3 ) References
Cement (clinker)116.753.83150[29,30,39]
Water2.301000[29,30,39]
Hydrates18.711.82073[26,27,28,41]
Air000[29,30,39]
Table 2. Physical parameters of the aggregate.
Table 2. Physical parameters of the aggregate.
PhaseBulk Modulus K (GPa)Shear Modulus G (GPa) Density ρ ( kg/m 3 ) References
Fine Aggregate36262650[27]
Coarse Aggregate41.619.22620[29,39]
Table 3. The parameters for prediction.
Table 3. The parameters for prediction.
w/c RatioVolume Fraction of Fine Aggregate in MortarVolume Fraction of Coarse Aggregate in Concrete
0.350.4040.5
0.50.3580.5
0.650.3110.5
Table 4. The parameters of experimental materials.
Table 4. The parameters of experimental materials.
GroupAge
(Days)
w/c RatioHydration DegreeVolume Fraction of Fine Aggregate in MortarVolume Fraction of Coarse Aggregate in Concrete
A70.550.7600.4910.482
280.550.8800.4910.482
B70.420.6760.3500.420
280.420.7800.3500.420
C70.370.6200.3190.420
280.370.7100.3190.420
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jiang, J.; Zhang, D.; Gong, F.; Zhi, D. Prediction of Ultrasonic Pulse Velocity for Cement, Mortar, and Concrete through a Multiscale Homogenization Approach. Materials 2022, 15, 3241. https://doi.org/10.3390/ma15093241

AMA Style

Jiang J, Zhang D, Gong F, Zhi D. Prediction of Ultrasonic Pulse Velocity for Cement, Mortar, and Concrete through a Multiscale Homogenization Approach. Materials. 2022; 15(9):3241. https://doi.org/10.3390/ma15093241

Chicago/Turabian Style

Jiang, Jingluo, Dawei Zhang, Fuyuan Gong, and Dian Zhi. 2022. "Prediction of Ultrasonic Pulse Velocity for Cement, Mortar, and Concrete through a Multiscale Homogenization Approach" Materials 15, no. 9: 3241. https://doi.org/10.3390/ma15093241

APA Style

Jiang, J., Zhang, D., Gong, F., & Zhi, D. (2022). Prediction of Ultrasonic Pulse Velocity for Cement, Mortar, and Concrete through a Multiscale Homogenization Approach. Materials, 15(9), 3241. https://doi.org/10.3390/ma15093241

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop