Next Article in Journal
Effect of Nanosilica on Mechanical Properties and Microstructure of PVA Fiber-Reinforced Geopolymer Composite (PVA-FRGC)
Next Article in Special Issue
Verification of the Influence of Loading and Mortar Coating Thickness on Resistance to High Temperatures Due to Fire on Load-Bearing Masonries with Clay Tiles
Previous Article in Journal
Cutting Characteristics for Sugar Maple Using Single Grit with Spherical Cone and Triangular Pyramid Geometries
Previous Article in Special Issue
The Effects of Nano-SiO2 and Nano-TiO2 Addition on the Durability and Deterioration of Concrete Subject to Freezing and Thawing Cycles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Model for Predicting the Tortuosity of Transport Paths in Cement-Based Materials

1
School of Materials Science and Engineering, Shijiazhuang Tiedao University, Shijiazhuang 050043, China
2
Foundational subjects department, Shijiazhuang TieDao University SiFang College, Shijiazhuang 051132, China
3
Key Laboratory of Smart Material and Structures Mechanics, Hebei Province, Shijiazhuang Tiedao University, Shijiazhuang 050043, China
4
School of Foreign Studies, Hebei Normal University, Shijiazhuang 050024, China
*
Author to whom correspondence should be addressed.
Materials 2019, 12(21), 3623; https://doi.org/10.3390/ma12213623
Submission received: 23 September 2019 / Revised: 27 October 2019 / Accepted: 31 October 2019 / Published: 4 November 2019
(This article belongs to the Special Issue Concrete and Construction Materials)

Abstract

:
The tortuosity of the pore structure is an important factor affecting medium (water and harmful ions) transport in cement-based materials. In this study, a new tortuosity model was established to reveal the effect of aggregate size, morphology, and graded media on the transport path in cement-based materials. Based on the stereological principle and the geometric algorithm, the distribution model of the ideal pebble and polygonal aggregate in cement-based materials was given first. Then, based on the image processing technology and MATLAB software, the morphology of the actual aggregate was also characterized to prove the similarity relationship between the ideal aggregate and actual aggregate. The reliability of the tortuosity model was verified by the mercury intrusion porosimetry test and data from other literature. Based on the tortuosity model, the influences of the aggregate particle shape parameters, hydration degree, and water-to-cement ratio on the tortuosity of the transport path were analyzed. Finally, the tortuosity model was further simplified to facilitate engineering application.

1. Introduction

The tortuosity in pore structure is one of the key factors affecting the transport performance of cement-based materials. It is usually defined as the ratio of the average diffusion path of a particle to the shortest distance of its straight line [1]. Many scholars have been devoted to testing pore tortuosity through experimental methods, however, it is relatively difficult to test directly. At present, the prediction of pore tortuosity is based on experimental statistical methods or empirical models and the reliability of the model is indirectly proven by experiments. Zhong et al. [2] established a model of pore tortuosity based on the average pore size measured by experiments. De Larrard et al. [3] established the pore tortuosity model based on the theory of an artificial neural network by using the existing experimental data as input data, combined with the numerical model. Piet et al. [4] established a tortuosity model based on the shape and gradation of spherical aggregate particles. The above models were obtained through the experimental statistical method, and the influence of the change of the aggregate’s own morphology on the transport path in cement-based materials was not fully considered, resulting in a certain under estimation between the theoretical prediction and the actual medium transport.
Alireza et al. [5] digitized the geometry and morphology of cracks by image analysis technology, and then established the tortuosity model of cracks. Akira et al. [6] used the three-dimensional spatial distribution model of calcium silicate hydrate (C–S–H) to establish the tortuosity model. Promentilla et al. [7] characterized the three dimensional pore structure based on x-ray synchrotron radiation micro-tomography technology, and established the pore tortuosity model. Zuo et al. [1] established a two-dimensional geometric model of the tortuosity of the transport path by analyzing an image of cement hardened paste by scanning electron microscopy. Liu et al. [8] and Zhou et al. [9] used fractal theory to analyze the tortuosity of capillary pore in cementitious materials, and established the corresponding tortuosity model. Li et al. [10] used macro and micro models to establish the tortuosity model of nano-channels in concrete. Liu et al. [11] provided a multi-scale transport tortuosity model ranging from the nano-scale through to the micro-scale and macro-scale of concrete, which was based on the porous media mechanics and multi-scale method, combined with numerical calculation and 3D visualization.
However, the above models were obtained by the image analysis method, and the influence of the randomness of aggregate distribution in concrete on the transport path of cement-based materials was not fully considered, resulting in a certain gap between the theoretical prediction and actual medium transport. In this study, based on the stereological principle and geometric method, tortuosity models of cement-based materials containing ideal pebbles and gravel polygonal aggregates were proposed, considering the morphology and gradation of aggregate. Second, the morphology of the actual aggregate was quantitatively characterized by image processing technology considering the effect of actual aggregate morphology on tortuosity, and finally the tortuosity model was validated by the mercury intrusion porosimetry test and the data reported in previous studies.

2. Theoretical Model of Tortuosity for Ideal Aggregate in Cement-Based

In this paper, mortar can be seen as the result of embedding fine aggregates into a cement slurry, and concrete can be seen as the result of embedding coarse aggregates into mortar. Due to the different sizes of the coarse and fine aggregates, we took the grading of the aggregate size as a weight function to homogenize the whole aggregate to obtain the average size of the aggregate. Then, the average size aggregate with the same size was uniformly embedded in the matrix, and the mortar and concrete model after the aggregate homogenization was obtained. The shape of the aggregate varies, as the fine aggregates are small, it is regarded as sphere; coarse aggregates are mainly divided into pebbles and gravel according to different sources, so it is regarded as ellipsoid and polyhedron, respectively [12], as shown in Figure 1. As the aggregate of the same particle size is evenly distributed in the matrix, the representative volume unit was taken, as shown in Figure 2.
For the pebble aggregate, the red solid line between points AB and CD is expressed in s ¯ , which is the average actual path length of the particles along the aggregate. The red dotted line between points AB and CD is expressed in l ¯ , which is the average displacement length between the start point and the end point. BC is the average distance between the pebble aggregates, expressed in o. For crushed stone aggregates, the red solid line between points DE and FG is expressed in s ¯ , which is the average actual path length of the particles along the aggregate. The red dotted line between point DE and FG is expressed in l ¯ , which is the average displacement length between the starting point and the end point. EF is the average distance between the crushed stone aggregates, expressed in λ ¯ . The average tortuosity is defined as the ratio of the actual path length of particles to the displacement length between the starting point and the end point of the particles. Combined with the literature [1], the specific expression of average tortuosity is shown in Equation (1):
τ c = s ¯ s t + λ s t ¯ l ¯ s t + λ s t ¯ s s a ¯ + λ s a ¯ l ¯ s a + λ s a ¯ τ c p
τ c p = η c ( 1 h a ) τ u h - c p + h a ( 1 + 7 h a ) w / c 0.35 η r τ h - c p
where s s a ¯ and s s t ¯ denote the actual transport path length of the particles moving along the sand aggregate surface and the stone aggregate surface, respectively; l ¯ s a and l ¯ s t denote the linear distance between the starting point and the end point of the sand and the stone, respectively; λ s a ¯ and λ s t ¯ denote the sand and the stone aggregate spacing, respectively; ha is the degree of hydration of the cement; τuh-cp and τh-cp are the tortuosity of the stacked cement particles and the tortuosity of the fully hydrated cement particles, respectively; and φuh-cp and φh-cp are the porosity of the stacked cement particles and the porosity of the fully hydrated cement particles, respectively, with φuh-cp = 0.58; w is the mass of water; c is the mass of cement.

2.1. Derivation of Tortuosity Model

As an indispensable and important part of concrete, the aggregate has a great influence on the corrosion resistance of concrete. The model treated aggregate size as an inclusion embedded in the cement slurry matrix. The shape of the aggregate was mainly pebbles and gravel, as shown in Figure 2.
The pebble aggregate particles are regarded as rotating ellipsoidal particles with different diameters, embedded in the mortar matrix. If the size distribution function fd of an elliptical aggregate is known, the short half-axis of the revolving ellipsoid is b and the long half-axis is a, k is the ratio of length to diameter of a revolving ellipsoid. When k = a/b > 1, the ellipsoid is prolate; when k = a/b < 1, the ellipsoid is oblate.
The volume (V) of the rotating ellipsoid is obtained as follows:
V ( b ) = { 4 π 3 k b 3 k > 1 4 π 3 k 2 b 3 k > 1
The surface area (S) of the rotating ellipsoid is obtained as follows:
S ( b ) = { 4 π k b 2 k > 1 4 π k 2 b 2 k < 1
the equivalent diameter of an aspheric particle is defined as the diameter of a sphere, which equals the volume of the aspheric particle. In other words, the particle size distribution of the non-spherical aggregate particles can be correlated with the particle size distribution of the equivalent spherical particles. The equivalent diameter (Deq), can be expressed as:
D eq = { 2 b k 1 3 k > 1 2 b k 2 3 k < 1
The gravel aggregate particles are regarded as regular octahedron and icosahedron. According to geometry and stereoscopic principles, the relationship between the edge length (a) and its circumscribed sphere radius (rcir) can be expressed as:
a = { 2 r c i r octahedron ( 5 1 ) r c i r 2 icosahedron
The volume (V) of the regular octahedron and icosahedron are given as follows:
V = { 2 3 a 3 octahedron 15 + 5 5 12 a 3 icosahedron
The surface area (S) of the regular octahedron and icosahedron are given as follows:
S = { 2 3 a 2 octahedron 5 3 a 2 icosahedron
For Plato polyhedral aggregate particles, the particle size of the polyhedral aggregate can be correlated with that of the spherical aggregate by means of an equivalent diameter. For a regular octahedron and icosahedron, its equivalent diameter can be expressed as:
D e q = { ( 1 π ) 1 3 r c i r octahedron ( 4 10 + 2 5 π ) 1 3 r c i r icosahedron
We mainly discuss two special grades of aggregate: Fuller grade aggregates and equal volume grade aggregates, which represent the upper and lower limits of aggregate size grading in concrete. Therefore, the size distribution function fd(x) of the polyhedral aggregate can be expressed as [13]:
f d ( x ) = q ( D max   e q q D min   e q q ) D e q q + 1 { q = 2.5 F u l l e r q = 3 E V F
The volume fraction (fst) of aggregate in concrete can be calculated by the concrete mix ratio. When the edge length of concrete blocks is known to be L, the number of aggregates (N) can be obtained.
N = f s t L 3 D min e q D max eq f d ( x ) V ( x ) d x
The average value of the surface spacing between adjacent aggregates can be expressed by the mean free path ( λ ¯ ) [14], which is directly affected by the aggregate volume fraction (fs) and solid specific surface area (SV), as shown in Equations (11) and (12):
S V = D min e q D max eq N f d ( x ) S ( x ) d x L 3
λ ¯ = 4 ( 1 f s ) S V
The average arc length ( l ¯ ) and chord length ( s ¯ ) of the ellipsoidal aggregate are given as follows, respectively.
l ¯ = 2 D e q min D e q max f d ( x ) b d x
s ¯ = π 4 [ k + 1 + 2 ( k 2 + 1 ) ] D e q min D e q max f d ( x ) b d x
The average arc and chord length of the regular octahedron and icosahedron aggregate are given as follows, respectively.
l ¯ = { D e q min D e q max 2 ( 1 π ) 1 3 f d ( x ) x d x octahedron D e q min D e q max 2 ( 4 10 + 2 5 π ) 1 3 f d ( x ) x d x icosahedron
s ¯ = { D e q min D e q max 2 2 ( 1 π ) 1 3 f d ( x ) x d x octahedron D e q min D e q max 5 ( 5 1 ) 2 ( 4 10 + 2 5 π ) 1 3 f d ( x ) x d x icosahedron
By substituting Equations (13)–(16) into Equation (1a), the tortuosity of cement-based materials can be expressed as follows:
τ m = s s a ¯ + λ s a ¯ l ¯ s a + λ s a ¯ τ c p
τ c = s ¯ s t + λ s t ¯ l ¯ s t + λ s t ¯ τ m
where τc, τm, and τcp are the tortuosity of the concrete, mortar, and porosity of cement pastes, respectively, and the expression is given in [1].
Porosity (φ) is also an important parameter affecting ion diffusion performance in concrete. The formula of porosity is as follows [15]:
φ = Max ( f v c w / c 0.39 h a w / c + 0.32 , 0 )
where w is the mass of water; c is the mass of cement; and fvc is the cement volume fraction.

2.2. Numerical Simulation of the Influence of Aggregate Morphology and Grades on Tortuosity

In theoretical simulation, the aggregate morphology is often regarded as regular octahedron and regular icosahedron, spherical, elliptical, and oblate pebble aggregates, as shown in Figure 3. According to GB/T 14865-2011 [16] in China, the length of pebble particles should be less than 2.4 times the average particle size of the corresponding particle size, and the thickness of pebble particles should be more than 0.4 times the average particle size. Therefore, the range of the length–diameter ratio of the pebble aggregate is 0.4~2.4.
According to Equation (17) as established in Section 2.1, for aggregates with Fuller and equal volume distribution, the influence of aggregate morphology including spherical, ellipsoid, octahedron, and icosahedron on tortuosity is shown in Figure 4. It can be seen that the shape, the volume fraction, and the gradation of aggregates are important factors affecting the change in the tortuosity in cement-based materials, and the tortuosity increased with the increase in the aggregate volume fraction.
It can be seen from Figure 4a that when the aggregate gradation and volume fraction were the same, the tortuosity value changes significantly with the morphology of gravel, and its size is sphere > icosahedron > octahedron in turn. Compared with spheres, the tortuosity values of icosahedral and octahedral aggregates were reduced by 11% and 15%, respectively. According to Equations (13)–(16), it is known that the closer the shape the aggregate is to a sphere, the larger the perimeter of the aggregate. With the increase in the perimeter of the aggregate in unit volume concrete, the length of the transport path increases, therefore, the tortuosity value of the concrete increases.
It can also be seen from Figure 4a that when the aggregate gradation and volume fraction are the same, the influence of the pebble morphology on the transport is simulated for the length–diameter ratio k = 0.4, 1, 1.2, and 2.4, representing the oblate, spherical, ellipsoid, and needle-bar aggregates. As can be seen from the figure, the tortuosity value was ellipsoid (k = 0.4) > sphere (k = 1) > ellipsoid (k = 1.2) > ellipsoid (k = 2.4). Compared with spherical aggregates, the tortuosity value of ellipsoid (k = 0.4), ellipsoid (k = 1.2), and ellipsoid (k = 2.4) increased by 18%, 10%, and 15%, respectively. The tortuosity value decreased with the increase in the length–diameter ratio. The main reason is that when the short axis of the ellipsoidal aggregate is the same as the radius of the spherical aggregate, according to Equations (13) and (14), as the length–diameter ratio increases, the quantity of the aggregate in the unit volume concrete increases, thereby increasing the perimeter of the aggregate and extending the transport path, therefore increasing the tortuosity value. Although the ellipsoidal aggregate in concrete with a length–diameter ratio (k) less than 1 has a larger tortuosity value, its compressive strength is lower, so it is used less in engineering.
In addition, from Figure 4a,b, it can be seen that the gradation of the aggregate also has a great influence on the tortuosity values. When the shape and the volume fraction of the aggregates are the same, the tortuosity of cement-based materials with Fuller gradation aggregates is greater than that with an equal volume gradation aggregates. For example, when the rotary ellipsoid and regular icosahedron k = 1.2 changed from 0.4 to 0.75, concrete tortuosity and Fuller gradation at the same volume gradation decreased by 7% and 9%, respectively. The equal volume aggregate gradation and Fuller aggregate gradation represent the upper and lower limits of concrete aggregate gradation, respectively. When the aggregate shape and volume fraction are the same, the amount of aggregate under the Fuller gradation in the unit volume concrete will be larger than that of the equal volume gradation. The more aggregates there are in the concrete, the longer the aggregate circumference sum, resulting in a longer ion transport path and the greater the tortuosity.

3. Morphological Characterization of Actual Aggregate in Cement-Based Materials

The tortuosity model established in Section 2.1 was mainly based on ideal aggregate particles, which has a certain gap with the actual aggregate morphology. The theoretical prediction in Section 2.2 also confirms that the aggregate morphology has a great influence on the pore tortuosity in cement-based materials. Therefore, it is necessary to carry out statistical characterization of the actual aggregate morphology to establish the similarity relationship between the ideal and actual aggregate. There are many parameters to represent the morphological characteristics of aggregates such as flatness, aspect ratio, micro-characteristics, roundness, practical sphericity, roughness, corner parameters, axial coefficient, roundness, etc. [17,18,19,20,21,22]. According to the parameters required for the tortuosity model established in Section 2.1, the parameters of the micro-characteristics were selected to characterize the morphology of the actual pebble and gravel aggregates.

3.1. Expressions of the Microscopic Characteristics of Actual Aggregate Particles

The microscopic characteristics refer to the ratio between the equivalent perimeter (Pe) and the real perimeter (P) of the aggregate particles. In this way, the microscopic characteristics of pebbles can be described as the ratio of the perimeter of the equivalent ellipse circumference of particles to the real perimeter of the particles. The gravels are described as the ratio of the perimeter of a regular decagon tangential to the perimeter of a particle with an equal area.
Thus, the microscopic characteristics (ξ) of pebbles and gravel can be expressed as follows:
ξ = P e P
Thus, the ξ of pebbles can be expressed as follows:
ξ = π 2 [ A π 1 e 2 + A 1 e 2 π + 2 [ A π 1 e 2 + A 1 e 2 π ] ] P
where A represents the sum of the pixels in the boundary area of particles; P represents the sum of the continuous pixels in the boundary area of particles; Pe represents the circumference of the equivalent ellipse of the object; and e represents the eccentricity of the equivalent ellipse.
Gravel is often regarded as octahedrons and icosahedrons, therefore, its projection shape is regular quadrilateral and decagon, respectively. The ξ of gravel is as follows:
ξ = { 2 2 D e q P octahedron   2 D e q P 10 + 2 5 icosahedron
According to Equations (20) and (21), the tortuosity of the actual aggregate concrete can be expressed as:
τ c = s s t ¯ + λ s t ¯ ξ s t ( l s t ¯ + λ s t ¯ ) τ m
τ m = s s a ¯ + λ s a ¯ ξ s a ( l s a ¯ + λ s a ¯ ) τ c p
where ξst is the microscopic characteristics of sand and ξsa is the microscopic characteristics of stone.

3.2. Characterization of Microscopic Properties of Actual Aggregate

The specific processing method is as follows:
(1) Sampling through the quartile method and laying it on a white background plane to avoid overlapping aggregate particles as far as possible.
(2) The images of coarse and fine aggregate particles are taken by a digital camera with more than 5 million pixels in the JPG file format. According to the JPG file format image obtained above, this paper used MATLAB programming software (MathWorks company of Beijing, China) to process the images of coarse and fine aggregate particles.
(3) The image in the JPG file format taken by digital camera is a RGB image that can be transformed into a gray image, which can greatly reduce the amount of image data and facilitate image processing, as shown in Figure 5a. By analyzing the gray histogram of the image to determine the optimal threshold needed for image segmentation, the threshold is used for binarization, and the gray image is transformed into a binary image, reverse binary image, and optimized image, as shown in Figure 5b–d.
(4) Then, the image area, edge circumference, eccentricity of the equivalent ellipse, the long axis of the equivalent ellipse, the short axis of the equivalent ellipse, and the diameter of the equivalent circle are extracted.
(5) Finally, the shape parameters of the aggregates are calculated by the weighted average method:
F = f s ξ 100
In the formula, F represents the shape parameters of the aggregate particles as a whole; fs represents the volume fraction of the aggregates per grain size; and ξ represents the shape parameter value of the corresponding granular aggregate.

3.3. Testing Process and Statistical Results of Micro-Properties of Fine Aggregate

After the obtained sand sample was uniformly stirred, the sample was divided by the quarter method, and after washing with water, six samples of not less than 550 g each were taken. The two samples were respectively placed in an oven at a temperature of 105 ± 5 °C and dried to a constant weight, and then 500 g of the sample was taken. The sieving results of the sample are shown in Figure 6.
According to the definition of the micro-characteristics and considering the actual situation, sand grains below 1.18 mm were regarded as spheres and their shape parameters were 1, where the sand particle size in the range of 1.18~2.36 mm and 2.36~4.75 mm were counted by the MATLAB programming language, according to the image processing steps given in Section 3.2. The statistical results of the micro-characteristic are shown in Table 1 and Table 2, respectively. It can be seen that the micro-characteristics of sand particles with particle sizes of 1.18~2.36 mm and 2.36~4.75 mm were between 0.95~0.98.

3.4. Testing Process and Statistical Result of Micro-Properties of Coarse Aggregate

The coarse aggregate was divided by a quartering method to 22 kg, and then pebble and gravel aggregate samples with particle sizes of 5~10 mm and 10~20 mm were randomly sampled as shown in Figure 7. For the image statistics, each aggregate sample was photographed from seven angles as shown in Figure 8. After image processing with MATLAB, the statistical results of the morphological parameters of pebbles and gravels are shown in Table 3, Table 4, Table 5, Table 6 and Table 7, respectively.
Table 3, Table 4, Table 5, Table 6 and Table 7 present the statistical morphological parameters of each aggregate including the area, circumference, eccentricity of equivalent ellipse, long and short axes, and equivalent diameter of circle by using MATLAB programming software from seven angles. These statistical morphological parameters were brought into Equation (22) and the average value of the shape parameters of the weighted aggregates were calculated and are given in Table 8. As can be seen from Table 8, the shape parameters of the aggregate with particle sizes of 5~10 mm and 10~20 mm ranged from 0.89~0.95.

4. Verification of Tortuosity Model

Due to the complexity of pore tortuosity, it is relatively difficult to measure directly. At present, mercury intrusion porosimetry (MIP) is commonly used to characterize the pore structure characteristics of cement-based materials such as porosity, tortuosity, and pore connectivity in cement-based materials [23]. Therefore, it was introduced in this paper to verify the reliability of the model established.

4.1. Preparation of Samples

Considering the possible errors caused by coarse aggregates in concrete, mortar specimens were selected to validate the model, where the water -o-cement ratio was 0.45, and the volume fraction of sand was 20%, 40%, and 50%, respectively. The gradation of sand is shown in Figure 6. The stirred slurry was first poured into a plastic beaker of 500 mL, shaken for 3 min, and then injected into a PVC (Polyvinyl chloride) tube with a diameter of 16 mm, and shaken for another 2 min. The purpose of vibration was to eliminate the bubbles produce d during mixing to reduce the test error. After three days of curing in the standard room, the sample was taken out from the PVC tube and then both ends were removed. The intermediate sample was taken to a height of about 15 mm and then standardized to the specified age for testing the porosity tortuosity. The sample was immersed in anhydrous ethanol before the test to prevent hydration.

4.2. Testing Process and Results of Pore Tortuosity

AutoPore IV9500 produced by Micrometrics (Ottawa, Canada) was used for MIP, where the maximum pressure can reach 415 MPa, and the pore size of the measurement range was 3 nm to 360 μm. The measurement was conducted in two stages: a manual low pressure from 0.003 MPa to 0.21 MPa, and automatic high pressure from 0.21 MPa to 242 MPa. After completing the low-pressure test, the penetrometer was taken out from the low-pressure chamber and weighed. The high-pressure test was carried out. The contact angle was set at 130 degrees and the balance time was 30 s.
The experimental and theoretical prediction results of tortuosity in cement-based materials are shown in Table 9, where the micro-characteristics of 2.36~4.75 mm particle size are 0.96. By substituting it into Equations (19)–(22), the theoretical results of actual tortuosity can be obtained. As can be seen from Table 9, the theoretical prediction results coincide with that of the MIP test, whose maximum error is 10.41%.
To further verify the reliability of the tortuosity model, the tortuosity value given by Zuo [1] was compared with the predicted results, as shown in Table 10. It can be seen from the table that the relative error between the measured and predicted values was less than 10%, which also prove that the prediction model of tortuosity and the characterization of morphological parameters are reasonable.

5. Sensitivity Analysis of Model Parameters

5.1. Effect of Aggregate Morphology on Tortuosity

Aggregate morphology is mainly expressed by the microscopic characteristics of aggregate particles (ξ), where we assumed that ξ = 0.5–1 and was brought into Equation (22). The relationship of the ξ and tortuosity is demonstrated in Figure 9. It can be seen that the tortuosity values of cement-based with pebble and gravel aggregates decreased with the increase in the ξ. For example, the tortuosity values at ξ = 0.5 decreased by 51% and 50% with respect to that at ξ = 1.0 for the pebble and gravel aggregate, respectively. The main reason is that the ξ values affect the roughness of the aggregate surface. The smaller the ξ values, the rougher the surface and the longer the perimeter of the aggregate, which increase the transport paths of the medium and are conducive to enhancing the durability of cement-based materials.

5.2. Effect of Hydration Degree on Tortuosity

The degree of hydration of cement (ha) has a great influence on the pore tortuosity of cement paste, where we assumed that ha = 0.5~0.9 and was brought into Equation (1b). The relationship of the ha and tortuosity is demonstrated in Figure 10. It can be seen that the tortuosity values of cement-based with pebble and gravel aggregates increased with the increase in the ξ. For example, the tortuosity values at ha = 0.5 decreased by 37.6% and 37.5% with respect to that at ha = 0.9 for the pebble and gravel aggregate, respectively. The main reason is the decrease in porosity and pore connectivity with the increase in hydration degree of cement, which leads to the increase in tortuosity.

5.3. Effect of Water–Cement Ratio on Transport

Water–cement ratio (w/c) has a great influence on the pore structure of cement matrix between the aggregate and matrix, where we assumed that w/c = 0.35~0.55 and was brought into Equation (1b). The relationship of the w/c and tortuosity is demonstrated in Figure 11. It can be seen that the tortuosity values of cement-based with pebble and gravel aggregates decreased with the increase in the ξ. For example, the tortuosity values at w/c = 0.35 decreased by 29.9% and 30% with respect to that at w/c = 0.55 for the pebble and gravel aggregate, respectively. The main reason is the decrease in porosity and pore connectivity with the increase in the hydration degree of cement, which leads to the increase in tortuosity.

6. Simplified Application of Tortuosity Model in Engineering

In practical engineering, the pore tortuosity is an important parameter for predicting the transport of medium in structural concrete. The tortuosity model established above fully considers many factors such as the morphology of the aggregate, volume fraction of the aggregate, aggregate gradation and hydration degree of the slurry, etc. There are many parameters and complex calculation processes, but in practical engineering, there is an urgent need for a simpler tortuosity model in the process. Therefore, the above-mentioned tortuosity model can be simplified, and its expression can be expressed as follows:
In practical engineering, it is necessary to simplify the tortuosity formula to facilitate the practical application of engineering. Table 11 is calculated according to Equation (1b). In Table 12, t1 represents the ratio of the circumference of the maximum section of the sphere, the ellipsoid of rotation, the regular octahedron, and the regular icosahedron to the circumference of the circumferential circle. t2 represents the ratio of the area of the maximum section of sphere, ellipsoid of rotation, regular octahedron, and regular icosahedron to the area of the circumscribed circle. t3 represents the ratio of the volume of sphere, ellipsoid of rotation, regular octahedron, and regular icosahedron to the volume of the circumscribed sphere.
The assumptions were as follows:
(1) The shape of the fine aggregate is a spherical shape; the shape of the coarse aggregate is a spheroidal shape and a regular polyhedron;
(2) Through the screening experiment, the grading of coarse and fine aggregates was statistically calculated, and the aggregate was weighted according to the grading and screening particle size to obtain the average particle size of the aggregate D s ¯ .
The average particle size of the fine aggregate is D s a ¯ . The average particle size of the coarse aggregate is D s t ¯ . R is the radius of the average particle size.
(1) The volume fraction of the aggregate fs is calculated by the concrete mix ratio.
(2) The volume fraction of fine aggregate is fsa. The volume fraction of the coarse aggregate is fst;
(3) In order to simplify the calculation, the regular polyhedron is unified into a rotating ellipsoid in combination with Table 12.
By introducing Equations (26) and (27) into (23) and (24), the simplified tortuosity formula can be expressed as follows in Equation (28):
λ = 4 ( 1 f s ) S V
S V = N t 2 4 k π R 2 ¯ L 3
N = 4 L 3 f s 3 π k t 3 R 3 ¯
s = t 1 2 k π D
τ c = t 1 D s t ¯ + 2 ( 1 f s t ) t 3 D s t ¯ 3 f s t t 2 ( D s t ¯ + 2 ( 1 f s t ) t 3 D s t ¯ 3 f s t t 2 ) π d ¯ + 2 ( 1 f s a ) D s a ¯ 3 f s a ( d ¯ + 2 ( 1 f s a ) D s a ¯ 3 f s a ) τ c p

7. Conclusions

(1) Based on the stereological principle, the tortuosity model of ideal pebble and gravel aggregates in cement-based materials was established by the geometric method. Then, the microscopic characteristics of the actual sand and coarse aggregate were counted by using the image processing software MATLAB. The results show that the microscopic characteristics of particle size were 0.95, 0.96, 0.93, and 0.92 in the range of 2.36–4.75 mm, 5–10 mm, and 10–20 mm, respectively. On this basis, the prediction model of the tortuosity degree of actual aggregate in the actual cement-based materials was established. The reliability of the tortuosity model was proven by the mercury intrusion test data and experimental data reported in the literature. The maximum error of the model was 8.59%.
(2) The effects of the aggregate particle shape parameters, aggregate volume fraction, and water cement ratio on the tortuosity of the transport path in concrete were analyzed. The tortuosity decreased with the increase in the aggregate particle shape parameters, increased with the increase in the cement hydration degree, and increased with the increase in the water cement ratio. Combined with engineering practice and considering the morphological characteristics of the aggregate, a simplified tortuosity model was presented, which facilitates engineering applications.

Author Contributions

Writing—original draft preparation, T.C.; writing—review and editing, G.S.; visualization, L.Z.; supervision, C.W.; methodology, Y.Z.; software, N.Y.; translation A.X.

Funding

This research was funded by the National Natural Science Foundations of China (No. 51478278), and the Higher Education Science and Technology Research Project of Hebei Province (No. ZD2016065).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zuo, X.B.; Sun, W.; Li, H.; Zhou, W.J. Geometrical model for tortuosity of transport paths in hardened cement pastes. Adv. Cem. Res. 2012, 24, 145–154. [Google Scholar] [CrossRef]
  2. Zhong, R.; Xu, M.; Vieira Netto, R.; Wille, K. Influence of pore tortuosity on hydraulic conductivity of pervious concrete: Characterization and modeling. Constr. Build. Mater. 2016, 125, 1158–1168. [Google Scholar] [CrossRef]
  3. De Larrard, T.; Benboudjema, F.; Colliat, J.B.; Torrenti, J.M.; Deleruyelle, F. Concrete calcium leaching at variable temperature: Experimental data and numerical model inverse identification. Comput. Mater. Sci. 2010, 49, 35–45. [Google Scholar] [CrossRef] [Green Version]
  4. Stroeven, P. A stereological approach to roughness of fracture surfaces and tortuosity of transport paths in concrete. Cem. Concr. Compos. 2000, 22, 331–341. [Google Scholar] [CrossRef]
  5. Akhavan, A.; Shafaatian, S.-M.-H.; Rajabipour, F. Quantifying the effects of crack width, tortuosity, and roughness on water permeability of cracked mortars. Cem. Concr. Res. 2012, 42, 313–320. [Google Scholar] [CrossRef]
  6. Hatanaka, A.; Elakneswaran, Y.; Kurumisawa, K.; Nawa, T. The Impact of Tortuosity on Chloride Ion Diffusion in Slag-Blended Cementitious Materials. J. Adv. Concr. Tech. 2017, 15, 426–439. [Google Scholar] [CrossRef] [Green Version]
  7. Promentilla, M.A.B.; Sugiyama, T.; Hitomi, T.; Takeda, N. Quantification of tortuosity in hardened cement pastes using synchrotron-based X-ray computed microtomography. Cem. Concr. Res. 2009, 39, 548–557. [Google Scholar] [CrossRef] [Green Version]
  8. Liu, X.W.; Yang, J.B.; Xia, K.Q.; Zhang, P.; Li, Z.G. Capillary Absorption Dynamics for Cementitious Material Considering Water Evaporation and Tortuosity of Capillary Pores. Adv. Mater. Res. 2013, 821, 1213–1218. [Google Scholar] [CrossRef]
  9. Zhou, Z.; Du, Z.; Chen, Z. Grout dispersion considering effect of pore tortuosity. Chin. J. Nonferrous Metals 2016, 26, 1721–1727. [Google Scholar]
  10. Li, W.; Guo, L.; Gao, L. Molecular dynamics simulation of chloride diffusion in concrete considering constrictivity and tortuosity of nanochannels. J. Southeast Univ. 2015, 45, 109–114. [Google Scholar] [CrossRef]
  11. Liu, Z.; Chen, W.; Zhang, Y.; Lv, H. A three-dimensional multi-scale method to simulate the ion transport behavior of cement-based materials. Constr. Build. Mater. 2016, 120, 494–503. [Google Scholar] [CrossRef]
  12. Roul, P.; Schinner, A.; Kassner, K. Mechanical properties of non-cohesive polygonal particle aggregates. Granul. Matter 2010, 13, 303–317. [Google Scholar] [CrossRef]
  13. Xu, W.; Chen, H.; Liu, L. Evaluation of Mesostructure of Particulate Composites by Quantitative Stereology and Random Sequential Packing Model of Mono-/Polydisperse Convex Polyhedral Particles. Ind. Eng. Chem. Res. 2013, 52, 6678–6693. [Google Scholar] [CrossRef]
  14. Xu, W.; Chen, W.; Chen, H.; Tian, X.; Zhao, H. Determination of overlapping degree of interfacial layers around polydisperse ellipsoidal particles in particulate composites. Phys. A Stat. Mech. Appl. 2014, 399, 126–136. [Google Scholar] [CrossRef]
  15. Clifton, J.R.; Bentz, D.P.; Ponnersheim, J.M. Sulfate Diffusion in Concrete; US Department of Commerce Technology Administration, National Institute of Standards and Technology: Gaithersburg, MD, USA, 1994.
  16. China National Institute of Standardization. Pebble and Crushed Stone for Construction; GB/T 14865-2011; CNIS: Beijing, China, 2011. [Google Scholar]
  17. Bowman, E.T.; Soga, K.; Drummond, W. Particle shape characterisation using Fourier descriptor analysis. Géotechnique 2001, 51, 545–554. [Google Scholar] [CrossRef]
  18. Ehrlich, R.; Weinberg, B. An exact method for characterization of grain shape. J. Sediment. Petrol. 1970, 40, 205–212. [Google Scholar]
  19. Gori, U.; Mari, M. The correlation between the fractal dimension and internal friction angle of different granular materials. Soils Found. 2001, 41, 17–23. [Google Scholar] [CrossRef]
  20. Hyslip, J.P.; Vallejo, L.E. Fractal analysis of the roughness and size distribution of granular materials. Eng. Geol. 1997, 48, 231–244. [Google Scholar] [CrossRef]
  21. Hentschel, M.L.; Page, N.W. Selection of Descriptors for Particle Shape Characterization. Part. Part. Syst. Charact. 2003, 20, 25–38. [Google Scholar] [CrossRef]
  22. Mora, C.F.; Kwan, A.K.H. Sphericity, shape factor, and convexity measurement of coarse aggregate for concrete using digital image processing. Cem. Concr. Res. 2000, 30, 351–358. [Google Scholar] [CrossRef]
  23. Laskar, M.A.I.; Kumar, R.; Bhattacharjee, B. Some aspects of evaluation of concrete through mercury intrusion porosimetry. Cem. Concr. Res. 1997, 27, 93–105. [Google Scholar] [CrossRef]
Figure 1. Aggregate shape. (a) Pebble aggregate; (b) gravel aggregate.
Figure 1. Aggregate shape. (a) Pebble aggregate; (b) gravel aggregate.
Materials 12 03623 g001
Figure 2. Definition of tortuosity. (a) Pebble aggregate; (b) gravel aggregate.
Figure 2. Definition of tortuosity. (a) Pebble aggregate; (b) gravel aggregate.
Materials 12 03623 g002aMaterials 12 03623 g002b
Figure 3. Aggregate shape in the tortuosity model. (a) Octahedron; (b) Icosahedron; (c) Sphere; (d) Prolate ellipsoid; (e) Oblate ellipsoid.
Figure 3. Aggregate shape in the tortuosity model. (a) Octahedron; (b) Icosahedron; (c) Sphere; (d) Prolate ellipsoid; (e) Oblate ellipsoid.
Materials 12 03623 g003
Figure 4. Numerical simulation results of the pore tortuosity. (a) Tortuosity under equal volume gradation. (b) Tortuosity under Fuller gradation.
Figure 4. Numerical simulation results of the pore tortuosity. (a) Tortuosity under equal volume gradation. (b) Tortuosity under Fuller gradation.
Materials 12 03623 g004
Figure 5. Analysis and calculation process of the aggregate parameters by the MATLAB program. (a) Gray image; (b) binary image; (c) reverse binary image; (d) optimized image.
Figure 5. Analysis and calculation process of the aggregate parameters by the MATLAB program. (a) Gray image; (b) binary image; (c) reverse binary image; (d) optimized image.
Materials 12 03623 g005aMaterials 12 03623 g005b
Figure 6. Sieving results of the sample.
Figure 6. Sieving results of the sample.
Materials 12 03623 g006
Figure 7. Actual aggregate morphology. (a) Actual pebble aggregate; (b) actual gravel aggregate.
Figure 7. Actual aggregate morphology. (a) Actual pebble aggregate; (b) actual gravel aggregate.
Materials 12 03623 g007
Figure 8. Aggregate projection shapes at different angles. (a) Pebble; (b) gravel.
Figure 8. Aggregate projection shapes at different angles. (a) Pebble; (b) gravel.
Materials 12 03623 g008
Figure 9. Effect of aggregate particle shape parameters on concrete tortuosity.
Figure 9. Effect of aggregate particle shape parameters on concrete tortuosity.
Materials 12 03623 g009
Figure 10. Effect of hydration degree on concrete tortuosity.
Figure 10. Effect of hydration degree on concrete tortuosity.
Materials 12 03623 g010
Figure 11. Effect of water-to-binder ratio on the tortuosity of concrete.
Figure 11. Effect of water-to-binder ratio on the tortuosity of concrete.
Materials 12 03623 g011
Table 1. Statistical results of sand particles for 1.18~2.36 mm sizes.
Table 1. Statistical results of sand particles for 1.18~2.36 mm sizes.
No.S1 #S2 #S3 #S4 #S5 #S6 #
Average value0.980.960.950.950.970.98
Standard deviation0.030.040.030.040.030.02
Table 2. Statistical results of sand particles for 2.36~4.75 mm sizes.
Table 2. Statistical results of sand particles for 2.36~4.75 mm sizes.
No.S1 #S2 #S3 #S4 #S5 #S6 #
Average value0.970.980.950.960.970.96
Standard deviation0.040.030.030.020.040.03
Table 3. Statistical results of the aggregate from different angles.
Table 3. Statistical results of the aggregate from different angles.
Aggregate ShapeNumber Area/pxPerimeter/pxEccentricityMajorAxis Length/pxMinor Axis Length/pxEquiv. Diameter/px
pebble11289133.510.7751.1932.5240.51
21814186.070.5053.4746.1748.05
32130169.040.6159.0746.6852.07
41660152.430.6453.2740.8745.97
51271244.920.7453.8936.2840.22
61668163.380.8766.2032.8446.08
71003140.430.8954.2524.3535.73
gravel12124177.610.8471.2538.5052.01
22077169.940.8167.0639.7251.42
31614181.470.7959.2936.2345.33
41478188.160.6852.0637.9943.38
51229140.760.7649.4532.1739.55
61481152.260.7453.3735.7143.42
71338133.730.5745.7337.4441.27
Table 4. Statistical results of the pebble aggregate for 5~10 mm particle sizes.
Table 4. Statistical results of the pebble aggregate for 5~10 mm particle sizes.
NumberArea/pxPerimeter/pxEccentricityMajorAxis Length/pxMajorAxis Length/pxEquiv. Diameter
13262226.270.4368.4161.8564.44
21725265.290.8169.1640.7946.86
32937201.790.6671.1953.2361.15
41888156.860.6255.6543.8449.02
53375228.080.8387.9649.6565.55
62108178.840.7464.2243.1451.80
71126126.660.8150.3929.3637.86
81518148.940.5448.8341.0243.96
91395143.560.6147.6337.7942.14
101350200.770.6548.9437.3241.45
Table 5. Statistical results of the pebble aggregate for 10~20 mm particle sizes.
Table 5. Statistical results of the pebble aggregate for 10~20 mm particle sizes.
NumberArea/pxPerimeter/pxEccentricityMajorAxis Length/pxMajorAxis Length/pxEquiv. Diameter
160,920961.320.67323.96241.61278.51
264,3361015.770.67334.66248.11286.21
361,8451530.840.71338.64238.64280.61
462,2771776.390.79386.19233.72281.59
554,9451221.350.69268.62226.81264.49
654,781975.180.76329.07214.08264.10
776,6261403.140.74389.22262.78262.35
864,8331719.120.71353.57250.07287.26
968,5521485.870.51323.22277.23295.44
1053,980939.8440.74322.72214.81262.16
1166,2001007.140.80378.02226.63290.32
1256,2751034.770.74328.29221.63267.68
1378,8871386.090.69375.78271.08266.93
1468,7711826.980.68353.69256.89295.91
Table 6. Statistical results of the gravel aggregate for 5~10 mm particle sizes.
Table 6. Statistical results of the gravel aggregate for 5~10 mm particle sizes.
NumberArea/pxPerimeter/pxEccentricityMajorAxis Length/pxMajorAxis Length/pxEquiv. Diameter
12335193.770.5059.4351.6454.52
21242128.940.5944.8436.2739.76
31267165.470.9369.6124.9140.16
41380150.430.5746.5738.3641.91
51590159.130.8461.8333.6544.99
62380251.170.8478.3242.5755.04
72504214.600.7167.6747.8056.46
81342184.980.5446.0638.6341.33
92127225.690.8371.5139.6252.04
101547158.880.7354.4637.1044.38
111688175.300.7960.8735.3346.08
Table 7. Analysis results of the gravel aggregate for 10~20 mm particle sizes.
Table 7. Analysis results of the gravel aggregate for 10~20 mm particle sizes.
NumberArea/pxPerimeter/pxEccentricityMajorAxis Length/pxMajorAxis Length/pxEquiv. Diameter
174,7681419.450.71368.95261.22308.54
297,9953081.720.64418.67321.58353.23
388,3531863.010.85464.77247.19335.40
470,3261515.840.83400.65224.80299.25
578,8791556.750.81426.32249.97266.91
684,3651292.030.79426.08256.52327.75
755,1382146.550.76346.99223.76264.96
870,0023779.760.76414.74266.80298.54
954,1351629.620.67325.22240.87262.54
1056,826908.230.70320.44228.80268.99
1152,974913.380.71266.23222.47259.71
1292,1251280.870.81454.03269.28342.48
1351,4722047.560.78340.04209.88256.00
1459,4541157.890.81367.35213.16275.13
Table 8. Average statistical results of the aggregate particle shape parameters.
Table 8. Average statistical results of the aggregate particle shape parameters.
Aggregate TypeAggregate SizeShape Parameter
pebble 5–10 mm0.92
pebble 10–20 mm0.95
gravel 5–10 mm0.91
gravel 10–20 mm0.89
Table 9. Comparison of the test and model prediction results.
Table 9. Comparison of the test and model prediction results.
NumberMeasured Average Porosity (%)Measured Average Tortuosity Simulated Average TortuosityError (%)
119.7835.7139.4310.41
215.4711.8312.919.12
312.3810.519.856.27
Table 10. Comparison of the test results and model results.
Table 10. Comparison of the test results and model results.
NumberWater Cement RatioMeasured Average Porosity (%)Measured Average TortuositySimulated Average TortuosityError (%)
SLZ0.350.3510.1916.1115.493.84
SLZ0.450.4512.8414.0712.868.59
SLZ0.550.5512.2614.2113.554.64
SLZ0.650.6510.5311.7210.867.33
Table 11. τcp ranges.
Table 11. τcp ranges.
w/c0.3–0.40.4–0.50.5–0.6
τ c p 17.47–14.3714.36–11.8211.82–9.73
Table 12. Range of values of the aggregate parameters ti for different shapes.
Table 12. Range of values of the aggregate parameters ti for different shapes.
Shapet1t2t3
Sphere111
Ellipsoid1~1.771~6.521~6.52
Octahedron0.90.550.32
Icosahedron0.980.760.61

Share and Cite

MDPI and ACS Style

Cao, T.; Zhang, L.; Sun, G.; Wang, C.; Zhang, Y.; Yan, N.; Xu, A. Model for Predicting the Tortuosity of Transport Paths in Cement-Based Materials. Materials 2019, 12, 3623. https://doi.org/10.3390/ma12213623

AMA Style

Cao T, Zhang L, Sun G, Wang C, Zhang Y, Yan N, Xu A. Model for Predicting the Tortuosity of Transport Paths in Cement-Based Materials. Materials. 2019; 12(21):3623. https://doi.org/10.3390/ma12213623

Chicago/Turabian Style

Cao, Tongning, Lijuan Zhang, Guowen Sun, Caihui Wang, Ying Zhang, Na Yan, and Aoxue Xu. 2019. "Model for Predicting the Tortuosity of Transport Paths in Cement-Based Materials" Materials 12, no. 21: 3623. https://doi.org/10.3390/ma12213623

APA Style

Cao, T., Zhang, L., Sun, G., Wang, C., Zhang, Y., Yan, N., & Xu, A. (2019). Model for Predicting the Tortuosity of Transport Paths in Cement-Based Materials. Materials, 12(21), 3623. https://doi.org/10.3390/ma12213623

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