Journal Pre-proof Effects of microstructure on desiccation cracking of a compacted soil
Qing Cheng, Chao-Sheng Tang, Hao Zeng, Cheng Zhu, Ni An, Bin Shi PII:
S0013-7952(19)31742-9
DOI:
https://doi.org/10.1016/j.enggeo.2019.105418
Reference:
ENGEO 105418
To appear in:
Engineering Geology
Received date:
16 September 2019
Revised date:
23 October 2019
Accepted date:
15 November 2019
Please cite this article as: Q. Cheng, C.-S. Tang, H. Zeng, et al., Effects of microstructure on desiccation cracking of a compacted soil, Engineering Geology (2019), https://doi.org/ 10.1016/j.enggeo.2019.105418
This is a PDF file of an article that has undergone enhancements after acceptance, such as the addition of a cover page and metadata, and formatting for readability, but it is not yet the definitive version of record. This version will undergo additional copyediting, typesetting and review before it is published in its final form, but we are providing this version to give early visibility of the article. Please note that, during the production process, errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
© 2019 Published by Elsevier.
Journal Pre-proof
Effects of microstructure on desiccation cracking of a compacted soil Qing Cheng1 , Chao-Sheng Tang*2 , Hao Zeng3 , Cheng Zhu4 , Ni An5,6 , Bin Shi7
Information of the authors 1. Assistant Professor, School of Earth Sciences and Engineering, Nanjing University, Nanjing,
f
Jiangsu Province, China. E-mail:
[email protected]
oo
2. Professor, School of Earth Sciences and Engineering, Nanjing University, Nanjing, Jiangsu
pr
Province, China. E-mail:
[email protected]
e-
3. MPhil student, School of Earth Sciences and Engineering, Nanjing University, Nanjing, Jiangsu
Pr
Province, China. E-mail:
[email protected]
4. Department of Civil and Environmental Engineering, Rowan University, 201 Mullica Hill Road,
al
Glassboro, New Jersey 08028, USA. Email:
[email protected]
rn
5. School of Earth Sciences and Engineering, Nanjing University, Nanjing, Jiangsu Province, China.
Jo u
E-mail:
[email protected]
6. Geoenvironmental Research Centre, Cardiff School of Engineering, Cardiff University, Newport Road, Cardiff CF24 3AA, UK
7. Professor, School of Earth Sciences and Engineering, Nanjing University, Nanjing, Jiangsu Province, China. E-mail:
[email protected] *Corresponding author: Prof. Chao-Sheng Tang School of Earth Sciences and Engineering, Nanjing University, Nanjing, Jiangsu Province, China. E-mail:
[email protected]
1
Journal Pre-proof
ABSTRACT Desiccation cracking has a significant influence on the hydro- mechanical behaviour of soils. Most previous studies focus on the desiccation cracking of slurry soil samples, whereas little attention has been paid to compacted soils. This study aims to investigate the effects of microstructure on the desiccation cracking of a compacted lean clay. Five soil samples are mixed with different water
f
contents including 12.5%, 14.5%, 16.5% (optimum water content), 18.5%, and 20.5%, compacted to
oo
generate different initial soil microstructures. After compaction, the soil samples are subjected to
pr
saturation and then the same drying process. The pore size distribution of each soil sample is
e-
characterized by performing mercury intrusion porosimetry (MIP) test. The change in water content
Pr
and the evolution of surface crack pattern during the drying process are continuously monitored. Experimental results show that the addition of water content during soil compaction significantly
al
influences the microstructure and desiccation cracking behaviour of soils. With increasing
rn
compaction water content from the dry side to the wet side of the optimum water content, soil
Jo u
microstructure transits from an aggregate structure to a dispersed structure, resulting in the change of pore size distribution from bimodal to unimodal. For soils with aggregate structures, the desiccation cracks initiate simultaneously and distribute uniformly throughout the soil body. With decreasing water content, crack geometrical parameters such as surface crack ratio and crack density increase almost linearly. Comparatively, for soils with dispersed structures, more localized growth of primary and secondary cracks are observed and the crack geometrical parameters show a two-stage linear growth during drying. This study provides microstructural interpretations to important aspects of desiccation cracking in compacted clayey soils and may guide the design of clay materials for geotechnical engineering applications.
2
Journal Pre-proof
Jo u
rn
al
Pr
e-
pr
oo
f
Keywords: microstructure; desiccation cracking; compacted soil; crack pattern; water content
3
Journal Pre-proof
1. Introduction During drying process, water evaporation may result in desiccation cracking of soils. The development of cracks has a significant influence on the mechanical and hydraulic properties of soils, leading to a variety of geotechnical and geological engineering problems (Morris et al., 1992; Miller et al., 1998; Yesiller et al., 2000; Albrecht and Benson, 2001; Rayhani et al., 2007; Tang et al., 2008;
f
Lozada et al., 2015; Chaduvula et al., 2017; An et al., 2018). For example, in landfill engineering, the
oo
desiccation cracking in barrier soils may result in more serious rainfall infiltration and landfill gas
pr
emission, further contaminating surrounding ground water and air (Li et al., 2019; Tang et al., 2019).
e-
In pavement engineering, the desiccation cracking in subsoil has a negative influence on the
Pr
long-term serviceability of roadways and endangers the safe and quality rides for the travelling public (DeCarlo and Shokri, 2014; Guo et al., 218). In slope engineering, the desiccation cracking
rn
et al., 2018; Jiang et al., 2019).
al
may increase the risk of rainfall infiltration and further reduce the factor of safety of the slope (Wang
Jo u
In the past few decades, many researchers paid attention to the qualitative and quantitative analysis of crack patterns by using different methods, such as image analysis and fractal analysis (Miller et al., 1998; Vogel et al., 2005; Tang et al., 2008; Lakshmikantha et al., 2009; Li et al., 2011; Liu et al., 2013; Lu et al., 2016; Zeng et al., 2019). Many efforts have also been devoted to investigate various influencing factors of the desiccation cracking of soils. The influencing factors can be separated into two main groups, namely soil internal factors and external environmental factors. The internal factors influencing desiccation cracking mainly include mineral composition, fine content, water retention characteristics, tensile strength, size and thickness of soil samples (Morris et al., 1992; Yesiller et al., 2000; Mitchell and Soga, 2005; Tang et al., 2008; Lakshmikantha
4
Journal Pre-proof
et al., 2012; Shorlin et al., 2000; Tollenaar et al., 2017). Temperature, boundary condition, relative humidity, wetting-drying cycles and other external environmental factors also affect desiccation cracking (Yesiller et al., 2000; Rodríguez et al., 2007; Groisman and Kaplan, 2007; Tang et al., 2010; DeCarlo and Shokri, 2014; Lakshmikantha et al., 2018; Levatti et al., 2019). In most of the previous studies, the experiments were conducted on slurry samples. However, in terms of engineering
f
practice, compacted soils, instead of slurry soils, are more often used as construction materials,
oo
which highlight the importance of this study and its potential benefits in guiding future engineering
pr
applications.
e-
Geotechnical properties of compacted soils can be significantly influenced by the initial
Pr
compaction state characterized by compaction water content and dry density (Nowamooz et al., 2013; Wang et al., 2014; Otalvaro et al., 2016; Jiang et al., 2017 ). This is because different compaction
al
states may result in different soil microstructures and further affect the hydro- mechanical behaviour
rn
of compacted soils (Delage et al., 1996; Mitchell and Soga, 2005; Romero and Simms, 2008 ). In
Jo u
terms of effects of compaction water content on soil microstructure, a major conclusion from previous studies is that an aggregate structure is well developed in the soil sample compacted on the dry side of the optimum water content. For the wet sample, clay particles forming a matrix envelop the silt grain and fill the inter-granular voids (Delage et al., 1996; Oualmakran et al., 2016; He et al., 2017; Zhang et al., 2018; Yuan et al., 2019). In terms of effects of dry density, a larger dry density may refer to the compression of inter-aggregate pores whereas the intra-aggregate pores almost remain unchanged (Delage and Lefebvre, 1984; Griffiths and Joshi, 1989; Wang and Xu, 2007; Yu et al., 2016). Obviously, desiccation cracking of compacted soils can also be affected by soil microstructure. However, so far, few of previous studies has investigated the effects of
5
Journal Pre-proof
microstructure induced by different compaction water contents on desiccation cracking of compacted soil samples. The principal objective of this study is to investigate the effects of different microstructure induced by various compaction water contents on desiccation cracking of a compacted soil. A series of desiccation cracking tests were carried out on compacted soil samples with various compaction
f
water contents. Moreover, mercury intrusion porosimetry (MIP) technique was used to investigate
oo
the pore size distribution (PSD) of each soil sample. Effects of soil microstructure on crack evolution
e-
pr
and final pattern of desiccation cracking were analysed and quantified.
Pr
2. Materials and methods 2.1 Testing materials
al
Xiashu soil taken from Nanjing, China is tested. It is an aeolian sediment and widely distributed
rn
in the Middle-Lower Yangtze Plain. The fractions of sand, silt and clay are 2%, 76% and 22%,
Jo u
respectively. The plastic and liquid limits are 19.5% and 36.5%, respectively. Index properties of the soil are summarized in Table 1. It is classified as a lean clay (CL) according to the Unified Soil Classification System (ASTM D-2487, 2017).
2.2 Sample preparation and testing procedure The tested soil was air-dried, crushed using a rubber pestle and passed through a 2-mm sieve. De-aired water was added to the sieved, air-dried soil to reach the target water contents of 12.5%, 14.5%, 16.5%, 18.5% and 20.5%, respectively. Then, the soil was sealed inside a plastic bag and kept for 48 hours to ensure the moisture equalization. Five soil samples (160 mm × 160 mm × 5 mm)
6
Journal Pre-proof
were prepared at different water contents. Static compaction was conducted using a loading machine, in which the axial displacement rate is controlled and fixed at 1.0 mm/min. The dry density is prepared to be 1.7 Mg/m3 by controlling the compaction height. After compaction, a desired amount of de-aired water was added to the compacted soil sample for saturation. The sample was then sealed using a plastic wrap and kept for another 48 hours for moisture equalization. After that, the plastic wrap was removed and the soil sample was setting up in the testing apparatus and air dried under
oo
f
room temperature (20±2o C). In addition, a slurry sample with an initial water content of 80% was
pr
also prepared for comparison in terms of the crack propagation manner.
e-
Moreover, small cubes with side length of 5 mm were prepared for the MIP tests to investigate
Pr
effects of compaction water content on pore size distribution. For the MIP tests, freeze-drying method was adopted (Delage et al., 1996; Tang et al., 2011c; Yu et al., 2016). The small cubic
al
samples were submerged in liquid nitrogen for a few minutes. Then, the frozen samples were put into
rn
a freeze-dryer with a temperature of -50℃ for 48 hours. During the MIP test, an absolute pressure p
Jo u
is applied on mercury to enter the soil pores (Romero and Simms, 2008). According to Washburn equation (Washburn, 1921), the pore diameter d can be calculated as follows: 𝑑=−
4𝜎𝐻𝑔 𝜃𝑛𝑤 𝑝
(1)
where 𝜎𝐻𝑔 is the surface tension of mercury (0.485 N/m at room temperature); 𝜃𝑛𝑤 is contact angle between mercury and the pore wall (162°, as suggested by Penumadu and Dean (2000)).
2.3 Testing apparatus The experimental setup used in Tang et al. (2010) was adopted in this study. The soil sample was placed on a balance to measure the weight with an accuracy of 0.01 g for monitoring change in water
7
Journal Pre-proof
content during the drying process. A digital camera was fixed above the sample to capture the images of surface crack pattern evolution with elapsed drying time.
2.4 Cracking image analysis The images taken by the digital camera were analysed using the software Crack Image Analysis
f
System (CIAS, available at www.climate-engeo.com). First of all, the image pre-processing
oo
including grayscale (Fig. 1(b)), binalization (Fig. 1(c)), denoising (Fig. 1(d)), skeletonising (Fig. 1(e))
pr
and crack identification (Fig. 1(f)) is carried out (Fig. 1). Then, the next step is quantitative analysis.
e-
Geometric indexes (e.g., length and width of each crack segment, angle of each intersection), general
Pr
indexes (e.g., number of intersections and crack segments, average crack length, average crack width, surface crack ratio) and statistical indexes (e.g., probability density function of crack length, crack
al
width and intersection angle) can be obtained through this step. In this study, some quantitative
rn
indexes are reported. Note that only the central part (150 mm × 150 mm) was used for the image
Jo u
analysis in order to minimize the boundary effects on cracking observations. The edge shrinkage area of each soil sample is counted when determining surface crack ratio (Zeng et al., 2019). More details of the procedure of the quantitative cracking image analysis are available in Tang et al. (2008) and Liu et al. (2013).
3
Experimental results and discussion
3.1 Microstructure of soil samples with different compaction water contents Fig. 2 shows the pore size distributions obtained from MIP tests for soil samples compacted at various water contents. As can be seen from the figure, for the soil samples compacted at water
8
Journal Pre-proof
contents of 12.5% and 14.5% (on the dry side of the standard Proctor optimum water content) and at water content of 16.5% (standard Proctor optimum water content), bimodal pore size distributions with two distinct pore domains can be observed. The two distinct pore domains are named as micro pore (or intra-aggregate pore) and macro pore (inter-aggregate pore) (Delage et al., 1996; Romero et al., 2011; Ng et al., 2016). On the dry side and at optimum water content, clay particles assemble
f
together and an aggregate structure is well developed in the soil. Both micro pores and macro pores
oo
can be detected by the MIP technique. The corresponding diameter at the peak point of the micro
pr
pores is about 0.6 μm and that of the macro pores is in the range from 10 to 30 μm. With increasing
e-
compaction water content, the amount of micro pores increases whereas the amount of macro pores
Pr
decreases because of the reduced aggregation of clay particles (Delage et al., 1996). However, for the soil samples compacted on the wet side of the standard Proctor optimum water content (at water
al
contents of 18.5% and 20.5%), the pore size distributions show a unimodal form. This is attributed to
rn
the assembly of clay-particle matrix around silt grains and the gradual formation of a dispersed
Jo u
structure (Delage et al., 1996). The corresponding diameter at the peak point is about 0.8 μm and there is no obvious separation between micro pores and macro pores. At various compaction water content levels, the soil samples with different microstructures are obtained. According to Alonso et al. (2013), a simple microstructural state variable ξm is adopted to quantitatively describe the microstructure of soil samples: 𝜉𝑚 = 𝑒𝑚 ⁄𝑒
(2)
where e is the total void ratio and em is the microstructural void ratio. The total void ratio e consists of the microstructural void ratio em and the macrostructural void ratio eM (Romero et al., 2011). In this study, the delimiting pore diameter, at which the derivative of PSD curve becomes zero, is used
9
Journal Pre-proof
to separate micro pores and macro pores (Ng et al., 2016). The delimiting pore diameters are found to be 1.9, 1.9, 2.4, 5.8 and 5.8 μm for the soil samples compacted at water contents of 12.5%, 14.5%, 16.5%, 18.5% and 20.5%, respectively. Fig. 3 shows the microstructural state variable ξm of each soil sample compacted at various compaction water contents, which is calculated based on the pore size distributions shown in Fig. 2. It can be found that with increasing compaction water content, the
f
microstructural state variable increases, corresponding to a less aggregated structure of clay particles
oo
(Delage et al., 1996). In this study, the MIP tests were conducted on the compacted soil samples. For
pr
a lean clay, the saturation process performed under zero stress does not significantly alter the PSD
Pr
e-
(Oualmakran et al., 2016).
3.2 Crack initiation and propagation in compacted samples with different microstructure
al
Under a drying condition, pore water evaporates from the surface of saturated samples. As the
rn
water content decreases, matric suction is developed and the suction-dependent tensile stress
Jo u
increases. Cracking occurs when the tensile stress exceeds the tensile strength (Tang et al., 2011a). With the elapsed drying time, the desiccation cracking further propagates. Fig. 4 shows the crack initiation and propagation of each soil sample with different microstructure. For the soil samples compacted on the dry side and at the optimum water content, crack initiates evenly on the soil surface and then widens gradually (Fig. 4 (a)-(c)). Despite of small crack width, the initiation of cracks occurs randomly and almost simultaneously on the surface of soil sample. In contrast, cracks evolve progressively in slurry samples, with the onset of primary cracks followed by the branching of sub-cracks (Fig. 5). This is because cracks usually initiate at the surface defects where shrinkage distortion and stress concentration occurs (Zabat et al., 1997; Weinberger, 1999; Tang et al., 2011b).
10
Journal Pre-proof
The soil samples compacted on the dry side and at the optimum water content have an aggregate structure, as discussed in previous section. During evaporation, soil volume shrinks as the inter-aggregate pores shrink and there is a rearrangement of aggregates (Nowamooz and Masrouri, 2010; Burton et al., 2015; Li et al., 2018). The strain energy build-up prior to cracking within a soil element subjected to tensile stress can be given as (1.299σ2 S2 t)/E, where σ is the tensile stress, S is
f
the pore diameter, t is the layer thickness and E is the Young’s modulus (Costa et al., 2013). The
oo
strain energy is proportional to the pore size. The macro pores between aggregates has a larger pore
pr
size, leading to higher strain energy prior to cracking. Thus, the places between aggregates can be
e-
regarded as the surface weak zones, resulting in crack initiation. Several crack initiation points
Pr
appearing almost at the same time, as marked by the red circles in Fig. 4 (a)-(c), can be observed. However, the slurry samples in the literature (Tang et al., 20 11a) have a more homogeneous
al
microstructure and shrink uniformly during the drying process. Desiccation cracking and strain
rn
energy concentration/release process can be maintained in an orderly manner (Tang et al., 2008) (Fig.
Jo u
5). Moreover, from one crack initiation point, three cracks propagate in three directions with an angle of 120° simultaneously, results in a “Y” shape crack intersection. This can be explained by adopting the concept of Griffith’s fracture energy balance. The “Y” shape crack intersection is the most energy efficient, where the crack length required is the minimum for expending the same amount of energy (Kodikara et al., 2000; Costa et al., 2013). With increasing compaction water content, soil microstructure transits into a less aggregated structure, resulting in a more localized cracking pattern. Furthermore, in soil samples compacted to the dry side of the optimum water content, a minor extent of sub-crack growth is observed (Fig. 4(a)-(c)), which is different from the slurry samples in the literature (Tang et al., 2011a) as well. This is likely because the intra-aggregate pores are almost
11
Journal Pre-proof
unchanged during drying (Nowamooz and Masrouri, 2010; Burton et al., 2015; Li et al., 2018). Moreover, no new cracks form and the additional strain is accommodated by further opening of existing cracks when a certain crack intensity is reached (Bai et al., 2000, Tang et al., 2006; Yin, 2010). For the soil samples compacted on the wet side with a dispersed structure, primary cracks initiate first (Fig. 4 (d)-(e)). The primary cracks of the soil sample compacted at water content of
f
18.5% initiate both at the inner surface and the edge, as shown in Fig. 2(d). With a higher
oo
compaction water content of 20.5%, the cracks initiates from the edge of soil sample, as shown in
pr
Fig. 4(e). This is because with a higher compaction water content, the microstructure is closer to the
e-
slurry samples. Therefore, the crack propagation manner of the soil samples with a dispersed
Pr
structure is more similar to that of the slurry samples (as shown in Fig, 5). After the primary cracks initiating, sub-cracks initiate at existing primary cracks and almost perpendicular to the primary
al
cracks (Fig. 4 (d)-(e)). This can be explained by adopting the maximum stress release criterion
rn
(Lachenbruch, 1962; Morris et al., 1992). When the primary cracks initiate, the internal tensile stress
Jo u
perpendicular to the existing primary cracks are released and the newly maximum tensile stress is parallel to the existing primary crack. Hence, the sub-cracks grow at the local maximum tensile stress and are perpendicular to the primary cracks. Quantitative analysis is carried out for the images showing development of desiccation cracking during the drying process. Fig. 6 shows the change of surface crack ratio Rsc with water content during drying. The surface crack ratio is defined as the ratio of the surface area of cracks to the total surface area of soil sample. It should be noted that the water content at the onset of cracking (indicated by the intersection point between surface crack ratio curve and water content axis) is defined as the cracking water content ωc (Tang et al., 2011b). It is found that with increasing
12
Journal Pre-proof
compaction water content from 12.5% to 20.5%, the cracking water content decreases from about 25% to 19%. The cracking water contents of the samples compacted on the dry side and at the optimum water content is higher than those of the samples compacted on the wet side. This is likely because the cracks in the samples compacted on the dry side and at the optimum water content appears in the void between aggregates while the cracks in the samples compacted on the wet side appears in the
f
void between soil particles, as discussed previously. According to the Young-Laplace equation, there
4𝑇𝑠 cos 𝜃 𝑑
(3)
pr
𝑠=
oo
is an inverse relationship between pore diameter d and suction s:
e-
where Ts is the surface tension of water (mN/m); θ is the contact angle (°). Larger pore diameter
Pr
develops smaller suction, resulting in higher amount of water retained in the soil. As a result, a higher cracking water content is observed in the soil samples compacted on the dry side. Moreover, it
al
can be seen from the figure that with decreasing water content during drying, the surface crack ratios
rn
of all the soil samples with various microstructure show an increasing trend. Specifically, for the soil
Jo u
samples compacted on the dry side and at the optimum water content with bimodal pore size distributions, Rsc increases almost linearly with decreasing water content. This is likely because for the soil samples compacted on the dry side and at optimum water content, all the cracks initiate almost simultaneously and the generated cracks gradually widen with further drying. For these three samples, the changing rate of Rsc and the final Rsc after reaching an equilibrium state decreases with decreasing compaction water content. The final Rsc of soil sample with compaction water content of 12.5% (0.056) is about twice than the soil sample with compaction water content of 16.5% (0.029). However, for the soil samples compacted on the wet side with unimodal pore size distributions, the relationship between surface crack ratio and water content before reaching the equilibrium state can
13
Journal Pre-proof
be regarded as consisting of two approximately linear segments. The changing rate of the first segment is much smaller than that of the second one. Actually, the first segment corresponds to the initiation of primary cracks and sub-cracks. The second segment corresponds to widening the existing cracks after the geometric crack pattern tends to stabilize (Fig. 4). Fig. 7 shows the change of crack density with water content during drying. Note that the crack
f
density (/cm2 ) represents number of crack segments per unit area. In general, the crack density
oo
increases with decreasing water content before reaching the equilibrium state. Similarly, for the soil
pr
samples compacted on the dry side and at the optimum water content with bimodal pore size
e-
distributions, the crack density increases almost linearly from the cracking water content to water
Pr
content of about 15%. For the soil samples compacted on the wet side with unimodal pore size distributions, the crack density increases slightly in the water content range from cracking water
al
content to 15%. With further drying to water content of 10%, a more dramatic increase in the crack
rn
density is observed. The corresponding water content at the intersection point is larger than that in
Jo u
Fig. 6 because the initiation of sub-cracks result in a significant increase in the crack density while do not count much in the surface crack ratio. The corresponding water content when the crack density becomes stabilised is larger than that for surface crack ratio (Fig. 6). This is because when the geometric structure of the crack network tended to stabilize and no new cracks were initiated, cracks become widen with further drying (Tang et al., 2011b). Moreover, the final crack density at the equilibrium state decreases with increasing compaction water content.
3.3 Final crack pattern of compacted samples with different microstructure Fig. 8 shows the final crack pattern of each soil sample. It can be found that there exist more fine
14
Journal Pre-proof
and thin cracks in samples with lower compaction water contents and vice versa. With increasing compaction water content, the geometric crack pattern tends to be simpler. This represents that there are less weak zones exist in a soil sample with a higher compaction water content. Compacted samples in this study exhibit similar geometrical shapes of soil clods (e.g., quadrangle, pentagon, and hexagon) as those observed in desiccated slurry samples (Yesiller et al., 2000; Rodríguez et al., 2007;
f
Groisman and Kaplan, 2007; Tang et al., 2010; DeCarlo and Shokri, 2014; Lakshmikantha et al.,
oo
2018; Levatti et al., 2019). In previous studies investigating desiccation cracking of slurry samples
pr
(Vogel et al., 2005; Péron et al., 2009; Tang et al., 2011a), the researchers get a common view that
e-
the crack segments are smooth and the soil surface is split into separate clods by the crack network.
Pr
However, as shown in Fig. 8, the crack segments are rough and the clods of compacted samples are not fully separated. This can be attributed to the homogenous structure of slurry samples and
al
heterogeneous structure of compacted samples.
rn
The final crack pattern of each sample is quantitatively analysed. The relationships between
Jo u
change of crack density as well as average width and microstructural state variable are plotted in Fig. 9. With increasing microstructural state variable, the crack density decreases almost linearly whereas the average width show an opposite trend. Moreover, it is found that on the wet side, with a given change in microstructural state variable, the change in average width is about 4.5 times more than that on the dry side. Crack width is one of the important indices to characterize the crack pattern. It is closely related to the crack propagation. For example, the widest cracks at the end of the drying process are the primary cracks (Tang et al., 2011b). Fig. 10 shows the probability-density function of crack width of each soil sample at the end of the desiccation test. Generally, with increasing compaction water
15
Journal Pre-proof
content, the probability density function of crack width changes from unimodal to bimodal, indicating that the grading characteristics of primary cracks and sub-cracks becomes more obvious. It can be seen from the figure that for the two samples compacted on the dry side (with water contents of 12.5% and 14.5%), the range of the crack width is relatively narrow, primarily distributed in the range from 0 to 0.10 cm. This confirms that for the soil sample with an aggregate structure, there is
f
no sub-cracks. For the soil sample compacted at optimum water content of 16.5%, the majority of
oo
crack width is still distributed in the range from 0 to 0.10 cm. However, a fraction of crack width is
pr
wider (0.10~0.25 cm), which is also observed in Fig. 8. This is likely because even though the soil
e-
sample compacted on the optimum water content has less obvious aggregates and a similar coating
Pr
clayey fraction to that found on the dry side (Delage et al., 1996). For the soil samples compacted on the wet side (with water contents of 18.5% and 20.5%), the crack with is mainly distributed in the
al
range from 0 to 0.20 cm and the distribution is relatively average. Extra-wide cracks (with width
rn
over 0.30 cm) appear on the surface of the sample compacted at water content of 20.5%. The
Jo u
observed phenomenon also confirms that both wide primary cracks and relatively thin sub-cracks exist in soil samples with a dispersed structure. Crack intersection angle is another important index for describing the crack pattern and also influenced by the crack propagation. The probability-density function of crack intersection angle is presented in Fig. 11. The probability-density function of crack intersection angle of each soil sample is basically a normal distribution. For the soil samples compacted on the dry side and at optimum water content, the intersection angles are mainly distributed in the range from 100 to 140°. This corresponds to the observed phenomenon that for soil samples with an aggregate structure, only primary cracks exist and they are formed simultaneously at the beginning stage of desiccation.
16
Journal Pre-proof Primary cracks always have “Y” shape intersections with an intersection angle of about 120° (Vogel et al., 2005; Péron et al., 2009; Tang et al., 2011a). For the soil samples compacted on the wet side with a dispersed structure, most of the intersection angles are in the range of 80~100° and the peak value of the probability-density function is at around 90°. This corresponds to the observation that for soil samples with a dispersed structure, sub-cracks initiates after the formation of primary cracks.
f
Sub-cracks normally show a “T” or “+” intersections with an intersection angle of about 90° (Vogel
pr
oo
et al., 2005; Péron et al., 2009; Tang et al., 2011a).
e-
3.4 General effects of microstructure on desiccation cracking of compacted soil
Pr
Based on the abovementioned phenomenon and interpretation, the crack propagations of soil samples with an aggregate structure and a dispersed structure show different manners. Fig. 12 and
al
Fig. 13 provide the schematic view of the desiccation cracking process in soils with aggregate and
rn
dispersed microstructures, respectively. For soil samples with an aggregate structure, the cracking
i.
Jo u
propagation can be summarized into two stages: Independent “Y” shaped cracks between aggregates initiate evenly on the soil surface with soil shrinkage during evaporation (Fig. 12 (a)-(b)). These cracks extend gradually until reaching the edge or intersecting with each other and the geometric structure of the crack pattern is formed (Fig. 12 (b)-(c)). ii.
The geometric structure of the crack pattern tends to stabilise when the water content continues to decrease. The existing cracks keep widening until the desiccation ceased (Fig. 12 c)-(d)).
For soil samples with a dispersed structure, desiccation cracking takes place in three stages:
17
Journal Pre-proof
i.
During the drying process, water evaporates and soil shrinks. Independent cracks between soil particles initiate at the inner surface or the edge (Fig. 13 (a)-(b)) and propagates until reaching the edge or intersecting with each other (Fig. 13 (b)-(c)). These initial independent cracks are considered as primary cracks.
ii.
Perpendicular to the existing primary cracks, some thin cracks initiate and propagate when
f
the water content keeps decreasing (Fig. 13 (c)-(d)). These post-developed cracks are
oo
regarded as sub-cracks. The primary cracks and the sub-cracks form the geometric structure
With further decreasing water content, the geometric structure of the crack pattern stays
e-
iii.
pr
of the crack pattern.
Pr
nearly unchanged while the primary cracks and the sub-cracks widen until the desiccation
Conclusions
rn
4
al
ceased (Fig. 13 (d)).
Jo u
Effects of soil microstructure on desiccation cracking of a compacted lean clay are investigated experimentally in this study. Soil samples with different microstructures induced by various compaction water contents are tested. The water content evolution, surface crack initiation and crack propagation are monitored during the drying process. Quantitative analysis of crack evolutions with the elapsed drying time leads to the following major conclusions: 1.
The pore size distributions of soil samples compacted on the dry side and at optimum water content are bimodal because of the aggregate structure, whereas the pore size of soil samples compacted on the wet side follows a unimodal distribution due to the dispersed structure. With increasing compaction water content, the microstructural state variable increases.
18
Journal Pre-proof
2.
For soil samples with an aggregate structure, independent cracks initiate evenly on the soil surface at the same time and then widen gradually. These cracks mostly intersect at an angle of 120°. For soil samples with a dispersed structure, primary cracks initiate first followed by the growth of sub-cracks along the orthogonal directions. The primary cracks formed “Y” shapes with intersection angles of 120° and the sub-cracks show “T” shapes with intersection angles of
The cracking water contents of the samples compacted on the dry side and at optimum water
oo
3.
f
90°.
pr
content is higher than those compacted on the wet side. Because the crack in the samples with an
e-
aggregate structure is likely to appear in the void between aggregates while the crack in the
Pr
samples with a dispersed structure appears in the void between soil particles. The macro pore has a larger diameter and develops smaller suction, resulting in higher amount of water retained
During the process of desiccation cracking, the surface crack ratio and the crack density of soil
rn
4.
al
in the soil.
Jo u
samples with an aggregate structure increase almost linearly with decreasing water content before stabilising. Soil samples with dispersed structures show a two-stage linear growth with decreasing water content. 5.
In terms of final crack pattern, with increasing compaction water content, the crack pattern of the soil sample becomes less complicated. There are more fine and thin cracks for the soil sample with a smaller microstructural state variable and vice versa.
Acknowledgements This work was supported by National Natural Science Foundation of China (Grant No.
19
Journal Pre-proof
41572246, 41772280, 41322019, 41902271), Natural Science Foundation of Jiangsu Province (BK20171228, BK20170394), and the Fundamental Research Funds for the Central Universities.
References Albrecht, B.A., Benson, C.H., 2001. Effect of desiccation on compacted natural clay. Journal of
f
Geotechnical and Geoenvironmental Engineering 127 (1), 67-75.
oo
Alonso, E.E., Pinyol, N.M., Gens, A., 2013. Compacted soil behaviour: initial state, structure and
pr
constitutive modelling. Géotechnique 63 (6), 463-78.
e-
An, N., Tang, C.S., Xu, S.K., Gong, X.P., Shi, B., Inyang, H.I., 2018. Effects of soil characteristics
Pr
on moisture evaporation. Engineering Geology 239, 126-135. ASTM D-2487, 2017. Standard practice for classification of soils for engineering purposes (Unified
al
Soil Classification System). American Society for Testing and Materials, West Conshohocken.
Jo u
403 (6771), 753-756.
rn
Bai, T., Pollard, D.D., Gao, H., 2000. Explanation far fracture spacing in layered materials. Nature
Burton, G.J., Pineda, J.A., Sheng, D., Airey, D., 2015. Microstructural changes of an undisturbed, reconstituted and compacted high plasticity clay subjected to wetting and drying. Engineering Geology 193, 363-373. Chaduvula, U., Viswanadham, B.V.S., Kodikara, J., 2017. A study on desiccation cracking behavior of polyester fiber-reinforced expansive clay. Applied Clay Science 142, 163-172. Costa, S., Kodikara, J., Shannon, B., 2013. Salient factors controlling desiccation cracking of clay in laboratory experiments. Géotechnique 63 (1), 18-29. Decarlo, K.F., Shokri, N., 2014. Effects of substrate on cracking patterns and dynamics in
20
Journal Pre-proof
desiccating clay layers. Water Resources Research 50 (4), 3039-3051. Delage, P., Audiguier, M., Cui, Y., Howat, M.D., 1996. Microstructure of a compacted silt. Canadian Geotechnical Journal 33 (1), 150-158. Delage, P., Lefebvre, G., 1984. Study of the structure of a sensitive Champlain clay and of its evolution during consolidation. Canadian Geotechnical Journal 21 (1), 21-35.
f
Griffiths, F.J., Joshi, R.C., 1989. Change in pore size distribution due to consolidation of clays.
oo
Géotechnique 39 (1), 159-167.
pr
Groisman, A., Kaplan, E., 2007. An experimental study of cracking induced by desiccation.
e-
Europhysics Letters 25 (6), 415-420.
Pr
Guo, Y., Han, C.J., Yu, X., 2018. Laboratory characterization and discrete element modeling of shrinkage and cracking in clay layer. Canadian Geotechnical Journal 55 (5), 680-688.
al
He, Y., Cui, Y.J., Ye, W.M., Conil, N., 2017. Effects of wetting-drying cycles on the air permeability
rn
of compacted Teguline clay. Engineering Geology 228, 173-179.
Jo u
Jiang, N.J., Tang, C.S., Yin, L.Y., Xie, Y.H., Shi, B., 2019. Applicability of Microbial Calcification Method for Sandy-Slope Surface Erosion Control. Journal of Materials in Civil Engineering, 31 (11), 04019250.
Jiang, Y., Chen, W.W., Wang, G.H., Sun, G.P., Zhang, F.Y., 2017. Influence of initial dry density and water content on the soil-water characteristic curve and suction stress of a reconstituted loess soil. Bulletin of Engineering Geology and the Environment 76 (3), 1085-1095. Kodikara, J.K., Barbour, S.L., Fredlund, D.G., 2000. Desiccation cracking of soil layers. Proceedings of the Asian conference on unsaturated soils: From theory to practice, Singapore, pp. 693-698. Lachenbruch, A.H., 1962. Mechanics of thermal contraction cracks and ice-wedge polygons in
21
Journal Pre-proof
permafrost. Geological Society of America, New York. Lakshmikantha, M.R., Prat, P.C., Ledesma, A., 2009. Image analysis for the quantification of a developing crack network on a drying soil. Geotechnical Testing Journal 32 (6), 505-515. Lakshmikantha, M.R., Prat, P.C., Ledesma, A., 2012. Experimental evidence of size effect in soil cracking. Canadian Geotechnical Journal 49 (3), 264-284.
f
Lakshmikantha, M.R., Prat, P.C., Ledesma, A., 2018. Boundary effects in the desiccation of soil
oo
layers with controlled environmental conditions. Geotechnical Testing Journal 41 (4), 675-697.
pr
Levatti, H.U., Prat, P.C., Ledesma, A., 2019. Numerical and experimental study of initiation and
e-
propagation of desiccation cracks in clayey soils. Computers and Geotechnics 105, 155-167.
Engineering Geology 123, 347-358.
Pr
Li, J.H., Zhang, L.M., 2011. Study of desiccation crack initiation and development at ground surface.
al
Li, Z.S., Derfouf, F.M., Benchouk, A., Abou-Bekr, N., 2018. Volume change behavior of two
rn
compacted clayey soils under hydraulic and mechanical loadings. Journal of Geotechnical and
Jo u
Geoenvironmental Engineering 144 (4), 04018013. Li, H.D., Tang, C.S., Cheng, Q., Li, S.J., Gong, X.P., Shi, B., 2019. Tensile strength of clayey soil and the strain analysis based on image processing techniques. Engineering Geology, 253, 137-148. Liu, C., Tang, C.S., Shi, B., Suo, W.B., 2013. Automatic quantification of crack patterns by image processing. Computers & Geosciences 57, 77-80. Lozada, C., Caicedo, B., Thorel, L., 2015. Effects of cracks and desiccation on the bearing capacity of soil deposits. Géotechnique Letters 5 (3), 112-117. Lu, Y., Liu, S.H., Weng, L.P., Wang, L.J., Li, Z., Xu, L., 2016. Fractal analysis of cracking in a
22
Journal Pre-proof
clayey soil under freeze-thaw cycles. Engineering Geology 208, 93-99. Miller, C.J., Mi, H., Yesiller, N., 1998. Experimental analysis of desiccation cracking propa gation in clay liners. Journal of the American Water Resources Association 34 (3), 677-686. Mitchell, J.K., Soga, K., 2005. Fundamentals of soil behaviour, 3rd edn. John Wiley, Sons, Inc, New Jersey.
f
Morris, P.H., Graham, J., Wiliams, D.J., 1992. Cracking in drying soils. Canadian Geotechnical
oo
Journal 29 (2), 263-267.
pr
Ng, C.W.W., Sadeghi, H., Hossen, S.K.B., Chiu, C.F., Alonso, E.E., Baghbanrezvan, S., 2016. Water
Pr
Geotechnical Journal 53 (8), 1258-1269.
e-
retention and volumetric characteristics of intact and re-compacted loess. Canadian
Nowamooz, H., Masrouri, F., 2010. Relationships between soil fabric and suction cycles in
al
compacted swelling soils. Engineering Geology 114, 444-455.
rn
Nowamooz, H., Jahangir, E., Masrouri, F., 2013. Volume change behaviour of a swelling soil
Jo u
compacted at different initial states. Engineering Geology 153, 25-34. Otalvaro, I.F., Neto, M.P.C, Delage, P., Caicedo, B., 2016. Relationship between soil structure and water retention properties in a residual compacted soil. Engineering Geology 205, 73-80. Oualmakran, M., Mercatoris, B.C.N., Francois, B., 2016. Pore-size distribution of a compacted silty soil after compaction, saturation, and loading. Canadian Geotechnical Journal 53 (12), 1902-1909. Penumadu, D., Dean, J., 2000. Compressibility effect in evaluating the pore-size distribution of kaolin clay using mercury intrusion porosimetry. Canadian Geotechnical Journal 37 (2), 393-405.
23
Journal Pre-proof
Péron, H., Herchel, T., Laloui, L., Hu, L.B., 2009. Fundamentals of desiccation cracking of fine-grained soils: experimental characterization and mechanisms identification. Canadian Geotechnical Journal 46 (10), 1177-1201. Rayhani, M.H.T., Yanful, E.K., Fakher, A., 2007. Desiccation- induced cracking and its effect on the hydraulic conductivity of clayey soils from Iran. Canadian Geotechnical Journal 44 (3),
f
276-283.
oo
Rodríguez, R., Sánchez, M., Ledesma, A., Lloret, A., 2007. Experimental and numerical analysis of
pr
desiccation of a mining waste. Canadian Geotechnical Journal 44 (6), 644-658.
e-
Romero, E., Della Vecchia, G., Jommi, C., 2011. An insight into the water retention properties of
Pr
compacted clayey soils. Géotechnique 61 (4), 313-328. Romero, E., Simms, P.H., 2008. Microstructure investigation in unsaturated soils: a review with
al
special attention to contribution of mercury intrusion porosimetry and environmental scanning
rn
electron microscopy. Geotechnical and Geological Engineering 26 (6), 705-727.
Jo u
Shorlin, K.A., de Bruyn, J.R., Graham, M., Morris, S.W., 2000. Development and geometry of isotropic and directional shrinkage-crack patterns. Physical Review E 61 (6), 6950-6957. Tang, C.A., Zhang, Y.B., Liang, Z.Z., Xu, T., Tham, L.G., Lindqvist, P.A., Kou, S.Q., Liu, H.Y., 2006. Fracture spacing in layered materials and pattern transition fro m parallel to polygonal fractures. Physical review E 73 (5) 49-72. Tang, C.S., Cui, Y.J., Tang, A.M., Shi, B., 2010. Experiment evidence on the temperature dependence of desiccation cracking behavior of clayey soils. Engineering Geology 114 (3-4), 261-266. Tang, C.S., Shi, B., Liu, C., Gao, L., Inyang, H., 2011a. Experimental investigation of the desiccation cracking behavior of soil layers during drying. Journal of Materials in Civil Engineering 23 (6),
24
Journal Pre-proof
873-878. Tang, C.S., Shi, B., Liu, C., Suo, W.B., Gao, L., 2011b. Experimental characterization of shrinkage and desiccation cracking in thin clay layer. Applied Clay Science 52 (1-2), 69-77. Tang, C.S., Shi, B., Liu, C., Zhao, L.Z., Wang, B.J., 2008. Influencing factors of geometrical structure of surface shrinkage cracks in clayey soils. Engineering Geology 101 (3-4), 204-217.
f
Tang, C.S., Tang, A.M., Cui, Y.J., Delage, P., Schroeder, C. and Shi, B., 2011c. A study of the
oo
hydro- mechanical behaviour of compacted crushed argillite. Engineering geology, 118 (3-4),
pr
93-103.
e-
Tang, C.S., Zhu, C., Leng, T., Shi, B., Cheng, Q., Zeng, H., 2019. Three-dimensional characterization
Engineering Geology, 255, 1-10.
Pr
of desiccation cracking behavior of compacted clayey soil using X-ray computed tomography.
al
Tollenaar, R.N., van Paassen, L.A., Jommi, C., 2017. Observations on the desiccation and cracking
rn
of clay layers. Engineering Geology 13 (1), 23-31.
Jo u
Volgel, H.J., Hofmann, H., Roth, K., 2005. Studies of crack dynamics in clay soil I: Experimental methods, results, and morphological quantification. Geoderma 125 (3), 203-211. Wang, L.L., Tang, C.S., Shi, B., Cui, Y.J., Zhang, G.Q., Hilary, I., 2018. Nucleation and propagation mechanisms of soil desiccation cracks. Engineering Geology, 238, 27-35. Wang, Q., Cui, Y.J., Tang, A.M., Li, X.L., Ye, W.M., 2014. Time- and density-dependent microstructure features of compacted bentonite. Soils and Foundations 54 (4), 657-666. Wang, Y.H., Xu, D., 2007. Dual porosity and secondary consolidation. Journal Geotechnical and Geoenvironmental Engineering ASCE 133 (7), 793-801. Washburn, E.W., 1921. A note on a method of determining the distribution of pore sizes in a porous
25
Journal Pre-proof
material. Proceedings of the National Academy of Sciences of the United States of America 7, 115-116. Weinberger, R., 1999. Initiation and growth of cracks during desiccation of stratified muddy sediments. Journal of Structural Geology 21, 379-386. Yesiller, N., Miller, C.J., Inci, G., Yaldo, K., 2000. Desiccation and cracking behavior of three
f
compacted landfill liner soils. Engineering Geology 57 (1-2), 105-121.
oo
Yin, H.M., 2010. Fracture saturation and critical thickness in layered materials. International Journal
pr
of Solids and Structures 47 (7-8), 1007-1015.
e-
Yu, C.Y., Chow, J.K., Wang, Y.H., 2016. Pore-size changes and responses of kaolinite with different
Pr
structures subject to consolidation and shearing. Engineering Geology 202, 122-131. Yuan, S.Y., Liu, X.F., Buzzi, O., 2019. Effects of soil structure on the permeability of saturated
al
Maryland clay. Géotechnique 69 (1), 72-78.
rn
Zabat, M., Vayer-Besançon, M., Harba, R., Bonnamy, S., Van Damme, H., 1997. Surface topography
Jo u
and mechanical properties of smectite films. Progress in Colloid and Polymer Science. Springer, Berlin, pp. 96-102.
Zeng, H., Tang, C.S., Cheng, Q., Inyang, H.I., Rong, D.Z., Lin, L., Shi, B., 2019. Coupling effects of interfacial friction and layer thickness on soil desiccation cracking behavior. Engineering Geology, 260, 105220. Zhang, T.W., Cui, Y.J., Lamas-Lopez, F., Calon, N., D'Aguiar, S.C., 2018. Compacted soil behaviour through changes of density, suction, and stiffness of soils with remoulding water content. Canadian Geotechnical Journal 55 (2), 182-190.
26
Journal Pre-proof
Figure captions Fig. 1 Crack image processing procedure Fig. 2 Pore size distribution of soil samples with various compaction water contents Fig. 3 Microstructural state variable determined for soil samples with various compaction water contents Fig. 4 Crack patterns captured at varying desiccation stages in soil samples compacted at (a)
oo
f
ω0 =12.5%, (b) ω0 =14.5%, (c) ω0 =16.5%, (d) ω0 =18.5% and (e) ω0 =20.5%
pr
Fig. 5 Progressive growth of primary cracks and sub-cracks in a slurry sample
e-
Fig. 6 Change of surface crack ratio with water content during drying
Pr
Fig. 7 Change of crack density with water content during drying Fig. 8 Final crack pattern of soil samples with various compaction water contents
al
Fig. 9 The evolution of crack density and average crack width with microstructural state variable in
rn
five soil samples
water contents
Jo u
Fig. 10 Probability-density function of crack width for five soil samples with different compaction
Fig. 11 Probability-density function of crack intersection angle for five soil samples with different compaction water contents Fig. 12 Schematic drawing of the development of desiccation cracks in the soil sample with an aggregate structure Fig. 13 Schematic drawing of the development of desiccation cracks in the soil sample with a dispersed structure
27
Journal Pre-proof
Table 1. Index properties of test soil
Property
Value
Specific gravity, Gs
2.73 2
Silt content (%)
76
Clay content (%)
22
oo
f
Sand content (%)
Liquid limit, ωl (%)
36.5
pr
Plastic limit, ωp (%)
Pr
Shrinkage limit, ωsl (%)
e-
Plasticity index, Ip (%)
19.5 17.0 10.5 16.5
Maximum dry density (Mg/m3 )
1.7
Jo u
rn
al
Optimum water content, ωopt (%)
28
Journal Pre-proof
Compaction water content affects microstructure and desiccation cracking of soils.
With increasing compaction water content, the microstructural state variable increases.
For soils with aggregate structures, the desiccation cracks initiate simultaneously and distribute uniformly.
For soils with dispersed structures, localized growth of primary and secondary cracks are
oo
With increasing microstructural state variable, the final crack pattern of the soil sample becomes
rn
al
Pr
e-
pr
less complicated.
Jo u
f
observed.
29