Introduction

The stress–strain response and volume change behavior of sands in shear and interface shear are mainly influenced by the normal stress, method of testing and testing materials [16], morphological characteristics of sand [712] and surface roughness of the geosynthetic material at the interface [1318]. Morphological characteristics of sand grains reflect the formation process, transportation process and material composition. These characteristics include angularity, roundness, sphericity, regularity and roughness. Several studies are available in literature that correlate the morphological characteristics of sand grains to the shear strength. Through systematic analysis of microphotographs of sand grains and oedometer tests, Cho et al. [8] observed that the increase in angularity of grains leads to increase in constant volume critical state friction angle and recommended that particle shape should be characterized and should be included in the classifications, especially for gravels and sands. Lim et al. [19] showed that particle shape and size have significant influence on the engineering properties of sands. Sand of different origins having wide range of particle size and morphological characteristics were studied through direct shear tests and image analysis and it was observed that the peak and ultimate friction angles as well as dilatancy were significantly affected by the particle shape and gradation of the sands. Santamarina and Cascante [20] studied the effect of particle roughness on the wave propagation through resonant column tests and observed that the surface roughness decreases contact stiffness and increases the interparticle friction.

Importance of understanding the interface shear behavior of sand with geosynthetics arises from the need for effective design of geotechnical structures like geosynthetic liners, retaining walls and embankments. Shear and bond stiffness at the interface depend on the properties of the soil and the geosynthetic material as well as the bonding characteristics of these materials. Many studies were reported on the influence of particle size and morphology on the shear behavior of sand–geomembrane interfaces. Fuggle [21] studied the effects of size and gradation of sand particles on the shear behavior of sand–geomembrane interfaces and observed that the large displacement friction angles and dilation angles were affected significantly by the particle size and gradation. Williams and Houlihan [13] showed that surface roughness, particle size, type and composition and the water content are the major parameters that can influence the interface shear behavior of soil-geomembrane. Similar other studies were reported on the effect of particle size [14], surface hardness of the polymer [15], surface roughness [16], and normal stress [17] on the interface shear behavior of sands. There are many studies in which the particle morphology and classifications were done by the image based technologies [2224]. The recent advancements in digital image technology provides an accurate estimation of the size and shape characteristics compared to conventional measurement techniques. Hence the present study adopts image analysis for characterizing the grain size and morphology, shear band thickness and surface roughness as well.

Though the effect of particle size on the shear strength behavior is studied well by earlier researchers, the effects of particle size and morphology are often combined and the individual effect of these parameters cannot be easily separated. Independent effect of particle size devoid of effects due to morphological characteristics was not studied by many researchers. The present study focuses on the effects of particle size on the shear and interfacial shear behavior of sands through direct shear tests, image analyses and roughness studies. Three different sands of different particle size were selected to study the effect of particle size on the shear and interface shear behavior of these sands. These sands were of same origin and hence the difference in the morphological characteristics of grains was minimal.

Materials Used in the Study

Sand

Three different types of sands were selected for this study. These sands were obtained by scalping specific size fractions from river sand of same origin. They are coarse sand (CS) (particle size 4.75–2 mm), medium sand (MS) (particle size 2–0.425 mm) and fine sand (FS) (particle size 0.425–0.075 mm), classified as poorly graded sands (SP) as per Unified Soil Classification System. In this paper, these sands as referred as CS, MS and FS. Grain size distribution curves for these sands are presented in Fig. 1. Figure 2 shows the photographs and scanning electron microscopic (SEM) images of sand for understanding the physical appearance and size variations of these sands. Properties of these three different sands are given in Table 1. The median size (D50) of three sands ranged between 3 and 0.22 mm, thus making them suitable for the study to understand the particle size effects on shear behaviour of sand alone and interface behavior of sand–geomembrane.

Fig. 1
figure 1

Grain size of all three sands used in this study

Fig. 2
figure 2

The photographs and scanning electron microscopic (SEM) images of sand used in study

Table 1 Properties of the sands used in the study

Geomembrane

A geomembrane made of smooth high-density polyethylene (HDPE) is used in present study. This is the most widely used geomembrane due to its favorable engineering properties like high tensile properties at low strain levels. Properties of geomembrane as given by the manufacturer are listed in Table 2.

Table 2 Properties of geomembrane used in the study

Morphological Characteristics of Sand Grains

Morphological characteristics of sand grains represent the overall geometry and shape of the sand grains. Some of these parameters selected for this study include angularity, roundness, sphericity and regularity. Definitions of these parameters are given by Pan [22], Pentland [25], Riley [26], Janoo [27], Wadell [28], Krumbein [29] and Sympatec [30]. Physical representation of these parameters along with the formulae used for their calculation are presented in Table 3.

Table 3 Morphological characteristics’ descriptors and their formulas along with figures and references are listed here

Recent developments in digital vision along with the availability of powerful image processing software has provided the means for the development of an alternate method for soil particle size analysis based on digital image processing. This technology had subsided the ocular based classification and also quantitative judgment made by the individual practicing engineers. In this study the SEM images of 50 particles of each type of sand were taken and used for the determination of their morphological characteristics. This process involves converting the SEM images to binary images and using the pixel information to obtain the geometrical parameters of the grains to calculate the required morphological characteristics as per Table 3. An algorithm was written in MATLAB for the digital image analysis of grains to obtain the morphological characteristics. The morphological characteristics thus determined for all three types of sand are given in Table 4. These parameters were measured from SEM images of sands taken at same magnification. According to the standard particle shape reference charts compiled by Santamarina and Cho [31], all the sands are classified as well rounded. The values of angularity, roundness, sphericity and roughness are nearly same for all the three sands, indicating that the morphological characteristics of these sands are similar.

Table 4 Morphological characteristics of sands

Direct Shear and Interface Direct Shear Tests

Design of Test Setup

Large size symmetric direct shear test (SDST) setup of sample size (300 mm length × 300 mm breadth × 160 mm height) was used in this study. Description of SDST setup is given in literature [1, 2]. Figure 3 presents the schematic diagram of the modified direct shear test set-up used in the present study. Several modifications were done to the conventional direct shear test set-up for conducting the tests presented in this paper. The upper half of the shear box is completely modified to obtain more uniform stresses and strains along the shear plane during the test as described in earlier studies [1, 2]. The visible front side of the upper half and lower half of the shear box were replaced with a transparent rigid acrylic sheet so as to visualize the particle movement and shear band formation during the tests. The vertical supporting rod for the loading yoke in the conventional direct shear box was replaced by a rectangular supporting frame connected to vertical rod to create an obstruction free foreground to the shear box for capturing images and videos (refer part 7 in Fig. 3). This test set-up has a facility to shear the sample up to a larger displacement (75 mm) to understand the post peak behaviour. In interface direct shear tests, the bottom half was replaced with a box consisting of an acrylic sheet attached to a metal sheet, which are supported on rigid concrete blocks (refer to BH in Fig. 3). One load cell and two LVDTs were used to measure the shear force and displacements respectively during the test. The setup is integrated with an automatic data acquisition system and graphical display.

Fig. 3
figure 3

Schematic diagram of the modified large size direct shear test setup

Sample Preparation and Testing

Sand was filled in the direct shear box in four layers of equal height. Relative density of sand was maintained as 70 % in all the tests. Weight of the sand required for each layer was calculated and poured into the shear box and was compacted to the predetermined layer thickness. The compaction was carried out by placing a wooden plank of plan size 250 × 250 mm and thickness of 20 mm on the sand sample on which light hammer blows were applied. Number of blows were gradually increased from bottom layer to top layer of the sample to ensure uniform compaction. To study the movement of particles during shear through image analyses, a few vertical colored sand columns were created along the transparent side of the shear box. To fill the colored sand, thin plastic channels of 2 mm thickness and 20 mm width were temporarily attached to the full-height of the transparent side of the shear box at approximately equal spacing to the shear box prior to the filling of sand in the box. Colored sand was filled in these channels through funnel to create colored sand columns up to a height of about 150 mm so as to surpass the shear band thickness. Relative density of the sand is maintained in the colored bands also. After the compaction is complete and the sample is prepared to the full height, the plastic channels were carefully withdrawn while ensuring that the bands remain intact and distinct from the surrounding sand. In case of interface direct shear tests, bottom half of the shear box was replaced by a rigid base with the provision to fix the geomembrane as explained in the description of test setup. Therefore only upper half of the box was filled with sand in two layers compacted to the required relative density. To avoid the movement the upper half of the shear box during compaction, it was robustly fixed to the platform of the shear box (refer to part 10 in Fig. 3).

After the sample preparation, the load pad is kept on the sample and normal load was applied by a lever arm mechanism. To obtain the symmetric loading conditions in the sample, the load pad was robustly fixed to the shear box with restraining screws available on front and rear sides (refer to part 9 in Fig. 3) of upper half of the shear box. The shear box was then set free by removing the holder restraining the upper half. All the tests were carried out at three normal stresses (σn) of 21, 37 and 58 kPa at a displacement rate of 1 mm/min. During the test the evolution of shear band formation in sand alone and sand–geomembrane is visualized and captured through the transparent side using a high definition camera of 12 megapixels.

Results and Discussion

Direct Shear Tests

A series of symmetric loading direct shear tests were conducted on sand. The shear stress–shear displacement responses of the sands under three different normal stresses are presented in Fig. 4. Corresponding volume change responses are presented in Fig. 5. These figures indicate that the stress–strain responses of the sands prior to the peak are almost same for all three sands, irrespective of the large variations in their particle sizes. All the three sands used in the present study have similar morphological characteristics, which is evident from Table 4. Santamarina and Cho [31] and Holtz and Kovacs [32] showed that the stress–strain response of the sands are influenced by the angularity and increasing angularity results in higher peak shear strength due to increase in interlocking forces. The angularity and roundness properties of the sands are almost same in the present study, indicating that the effect of morphology is masked and the variations in the shear behavior are only due to the particle size effects. According to Holtz and Kovacs [32], particle size does not have any effect on the peak friction angle (ϕp) if the void ratio is same. For the present samples, void ratios were calculated at 70 % relative density and the same are presented in Table 5. From the table, it is observed that the void ratios of coarse and medium sands are same, whereas the void ratio for FS is slightly higher. Table 1 also presents the range of void ratios for each type of sand. The difference in minimum and maximum void ratios is increasing with the decrease in particle size. However, for 70 % relative density, the difference in the void ratios is not significant, hence the effect of particle size on peak shear strength is negligible.

Fig. 4
figure 4

Shear stress versus shear displacement response of different sands. (Color figure online)

Fig. 5
figure 5

Vertical displacement versus shear displacement response of different sands. (Color figure online)

Table 5 Void ratios and unit weights of different sands at a relative density of 70 %

The post-peak stress–strain response of the coarse, medium and FSs is considerably different as observed from Fig. 4. In case of CS, the post peak drop in shear strength is less compared to medium and FSs. FS showed maximum post peak strain softening. The post peak strain softening is increasing with the decrease in particle size, though the relative density of the samples was same in all three cases. Herle and Gudehus [33] and Atkinson [34] showed that the ultimate friction angle (ϕr) in direct shear tests is equal to the angle of repose. To confirm this, angle of repose of the coarse, medium and FSs were determined as per ASTM C1444 [35]. In this method, a glass funnel filled with dry sand is held at a specific height to form a sand heap until the funnel is in contact with the heap. This method was used for measuring the angle of repose of sands as well as gravels by several earlier researchers with modifications [36, 37]. Typical measurement of angle of repose for the MS is shown in Fig. 6. After the heap is formed, high definition images were captured in four mutually perpendicular directions and the average angle of repose was estimated. Table 6 gives the values of measured angles of repose for all three types of sands along with their comparison with the ultimate friction angle obtained from direct shear tests. Comparisons showed that the ultimate friction angles were in good agreement with the angles of repose for sands of different grain sizes. Unlike the peak friction angles, the ultimate friction angles were influenced by the size of the particles.

Fig. 6
figure 6

Determination of angle of repose (ϕa) for medium sand

Table 6 Peak friction angle, ultimate friction angle and angle of repose values for different sands

One of the very important aspects of plane strain shearing in sand is the formation of shear bands because they can be related to the failure and strain localization. Well defined shear bands were formed in the direct shear tests with all three types of sands. The shear zone thickness is a significant factor which is correlated with the strain necessary to obtain the failure [3840]. Thickness of shear bands were determined for the direct shear tests through image segmentation technique explained by Vangla and Latha [2]. Videos of deforming colored sand columns were captured and the images at regular displacement intervals (about 0.1 mm) were extracted and used to find the shear band thickness. Figure 7 shows the deformed (20 mm) and initial (0 mm) profiles of colored sand columns for CS traced and overlapped one over the other so as to see the net deformation and shear band thickness (t) before and after the test. Similar exercise was done for all other sands and it was observed that the thickness of shear bands is increasing with the particle size. Literature [4143] suggests that the average thickness of the shear zone varies between 5 and 25 times the particle size. It has been also observed that the ratio between shear zone thickness and soil particle size depends on the scale of the test (ratio between sample size and particle size), with values considerably larger than 25 having been observed in large scale tests [44]. The shear band thickness measured in the present study is agreeing with these observations. Variation of normalized shear band thickness (t/D50) with the median particle size (D50) is shown in Fig. 8. As observed from Fig. 8, normalized thickness of shear bands is decreasing with the increase in D50, which is in agreement with the observations of earlier researchers [45].

Fig. 7
figure 7

Estimation of shear band thickness from initial and final profiles of colored sand columns. (Color figure online)

Fig. 8
figure 8

Variation of normalized shear band thickness with D50

Interface Direct Shear Tests

Interface shear tests were carried out on sand–geomembrane interfaces in SDST setup. All three types of sands were used in these tests and the tests were carried out under three normal stresses of 21, 37 and 53 kPa. The shear stress-shear displacement responses of the sand–geomembrane interfaces under three different normal stresses are presented in Fig. 9. The shear strength of all the sand–geomembrane interfaces is less than the corresponding sand alone shear strength. As observed from Fig. 9, MS with D50 of 0.87 mm is showing highest interfacial frictional resistance with the geomembrane compared to the coarse and fine sands. Also, there is considerable post peak drop in the shear strength in case of FS, unlike coarse and medium sands where the peak strength is sustained for large deformations. Figure 10 presents the comparison on peak and post peak interfacial shear strengths of different sand–geomembrane interfaces, which clearly shows that the drop in peak strength is minimal for coarse and medium sands and significant for FS. The calculated interface friction angle is 0.4 times the ϕp for CS, 0.49 times ϕp for MS and 0.37 times ϕp for FS.

Fig. 9
figure 9

Shear stress versus shear displacement response of different sands with GM from interface direct shear tests. (Color figure online)

Fig. 10
figure 10

Peak and post peak interface shear stress versus normal stress from interface direct shear tests. (Color figure online)

To understand this behavior, microscopic images of geomembrane before and after the test were examined for their surface roughness. Figure 11 presents the microscopic images of geomembrane before the test and the images after the tests for coarse, medium and FSs carried out at 37 kPa. These images clearly show that the CS could make strong impressions leading to plowing and wearing of the surface. Plowing was also evident in case of MS. However, FS could only make lighter impressions on the geomembrane. To quantify the effect of particles on the interfacial friction, both the number of contacts per unit area and the groove depth need to be considered. Groove depth depends on the particle size and angularity as well as the contact stress per particle, which in turn depends on the roughness and hardness of the particles [17]. In the present study, the roughness and hardness of the sand particles is almost same. The only factor that can influence the interfacial friction in this case is the particle size. A combined effect of number of grooves per unit area and groove depth is realized in the present study. MS showed highest interfacial friction because of increased number of contacts per unit area along with groove depths comparable to CS. In case of CS, though the groove depths are more, number of contacts per unit area were less because of bigger particle sizes. Though the number of contacts per unit area were highest in case of FS, the particles were able to slide along the interface without scarring the surface much as coarse and medium sands do, this can be seen in Fig. 11. Thus the effective number of contacts per unit area (number of particles in contact with the surface and causing grooves of considerable depth in the interfacing material) is the parameter that governs the interfacial shear strength of sands.

Fig. 11
figure 11

Microscopic images of geomembrane before and after the interface shear tests with different sands at a normal stress of 37 kPa

Further, surface roughness studies were carried out on the geomembranes before and after the test. A skidless stylus profilometer with maximum traverse length of 50 mm and tip diameter of 2 μm was used for this purpose. The vertical range of this profilometer is 500 μm. Representative specimens of size 100 × 20 mm were cut from the geomembranes and the surface roughness was determined in the direction perpendicular to the direction of shearing to get the closest estimates of the profile roughness. The measurements were made at four distinct locations for a length of 40 mm in each specimen and the average value was taken. Figure 12 shows the photograph of the stylus profilometer with a specimen in place. Roughness of the specimens is represented by root mean square (RMS) roughness value (Rq), which is defined as the geometric average height of roughness component irregularities from the mean line measured within the sampling length. This value is obtained directly from the measurement using profilometer. Variation of Rq with the normal stress for geomembranes used in interface shear tests with all three types of sands are plotted in Fig. 13. As observed from figure, the surface roughness values increased with the increase in normal stress for all cases because of higher stress applied per unit area. Surface roughness of the geomembrane used in tests with MS was measured to be highest compared to the other two cases because of higher number of contacts per unit area complemented by the particle size, causing grooves of considerable depth. Geomembrane specimen used in the test with FS showed least surface roughness because of the smaller particle size, leading to the sliding of particles without forming deeper grooves. These measurements are in agreement with the visual observations from the microscopic images of these membranes.

Fig. 12
figure 12

Stylus profilometer with geomembrane specimen

Fig. 13
figure 13

Variation of Rq for geomembrane after interface tests with sands of different particle size

Videos captured during the interface shear tests showed that there was no formation of shear bands for sand–geomembrane interfaces because the particles were sliding or being dragged along the interface because the surface of the geomembrane is very smooth. The interface shearing is non-dilative for smooth geomembranes with extremely low vertical deformations as explained in earlier studies of interface shear tests with different geosynthetics [2, 46]. In case of interface shear tests, the governing factor for interfacial shear strength is the relative roughness of the sand and geomembrane which leads to the plowing and formation of grooves rather than the formation of shear bands as observed in sand alone tests.

Conclusions

From the series of direct shear tests on three different sands with different particle sizes and similar morphological characteristics, interface shear tests on these sands with geomembrane and image analysis techniques, the following major conclusions are drawn from the present work.

  1. (1)

    Peak shear strength of sands with similar morphological characteristics is unaffected by the particle size of sand, when the tests are conducted at same void ratio.

  2. (2)

    Post peak shear stress response of sands varies with the change in particle size even when the tests are conducted on sands of similar morphology at same void ratio. Ultimate friction angle increases with the particle size, which is confirmed from the angle of repose tests.

  3. (3)

    Thickness of the shear bands in direct shear tests increases with the increase in particle size. However, the normalized shear band thickness (t/D50) increases with the decrease in particle median size (D50) of sands.

  4. (4)

    Interface shear strength does not depend on the particle size at a specific void ratio but is influenced by the number of effective contacts per unit area. Among the three sands tested with geomembrane in interface shear, MS showed higher shearing resistance due to increase in number effective contacts, which is verified through microscopic image analysis and surface roughness studies on the geomembranes before and after the tests.