International Journal of Geosciences
Vol.4 No.10A(2013), Article ID:41562,10 pages DOI:10.4236/ijg.2013.410A004

Advanced Statistical Analysis for Relationships between Particle Morphology (Size and Shape) and Shear (Static and Dynamic) Characteristics of Sands

Devapriya Chitral Wijeyesekera1, Alvin John Lim Meng Siang1, Ahmad Shukri Bin Yahaya2

1Department of Civil and Environmental Engineering, Universiti Tun Hussein Onn Malaysia, Batu Pahat, Malaysia

2School of Civil Engineering, Universiti Sains Malaysia, Penang, Malaysia

Email: devapriya@uthm.edu.my, alvinjlms@gmail.com, shukri@eng.usm.my

Copyright © 2013 Devapriya Chitral Wijeyesekera et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. In accordance of the Creative Commons Attribution License all Copyrights © 2013 are reserved for SCIRP and the owner of the intellectual property Devapriya Chitral Wijeyesekera et al. All Copyright © 2013 are guarded by law and by SCIRP as a guardian.

Received October 28, 2013; revised November 25, 2013; accepted December 11, 2013

Keywords: Morphology; Sand; Statistics; Static and Dynamic Shear

ABSTRACT

Morphology properties of naturally occurring sandy soil deposits are a consequence of the past geological transport history which gives rise to their different shape and size distributions. Five sand samples with different particle size and shape distributions were analysed statistically. Cluster analysis, performance indicators and probability distributions were used to find a representative shape parameter value. Static and dynamic tests were also contemporarily done to study the shear characteristics of the samples. Dimensionless relationships of the critical shear stress (τcr) normalised with respect to the normal stress (σN) against the representative values are presented. Dynamic testing on the samples showed that the modulus reduction curve (G/Gmax) and damping ratio (D) were also dependent on the morphology characteristics of the sand samples.

1. Introduction

Soil mechanics/science is a prime academic discipline of civil engineering importance, particularly in foundation design and hazardous ground analysis. The mechanical behaviour of sand is highly dependent on the interparticle reaction and therefore classification based on the morphological aspects is deemed necessary. Soils differ from other engineering materials because of its morphology characteristics which require a distributed particulate mechanics approach. For sandy soils (assumed cohesionless), the interparticle reactions of sands rely heavily on its initial density, initial confining pressure and most importantly the interlocking structure of the particles. The friction angle (Ø) is dependent on the relative density (Dr) and it is correlated with a linear relationship Ø = ADr + B, where A and B are constants [1].

Sand naturally occurs as granular material composed usually of broken rock fragments and mineral particles. The properties largely depend on the rocks from which they are derived. Different morphology (size and shapes)

characteristics are imposed on the sand particles as a consequence of the prevailing geological transport mechanism such as fluvial, aeolian or glacial process [2]. Table 1 explains this further in indicating different particle shapes based on its genetic transport mechanisms.

The study of the significance of particle morphology on the response to its mechanical behaviour needs further research. BSCS, ASTM, and AASHTO classifications do not consider particle shapes in their respective soil classification systems. This has resulted in only very few established theories being based on the influence of particle shape or “combined size and shape” relationships on the static and dynamic loading behaviour of granular soils [3-6].

The key aim of this study was to assess thoroughly the classification of particle morphology using state of the art experimental and statistical analysis techniques. The findings were put into a wider context on how it can have a significant influence on the static and dynamic parameters in geotechnical design.

Table 1. Description of the occurrence of particle shape based on its transportation mechanisms [2].

2. Design of Experiment, Materials and Testing Apparatus

The design of experiment is important to carefully plan the type of materials used, the experiments required, and the amount of data needed for the study. The following paragraphs explain the design incorporated in this study.

Five types of sand samples were used to study the morphology characteristics of the sand particles. The first sample is the commercially available natural river sands from Kahang, Johor, Malaysia which were used as the primary material classed as well graded sand (SW). This was strategically chosen and subsequent alteration through riffling of the SW sand was adopted to obtain the second and third sample of uniformly graded (SPuKahang) and gap graded (SPg) sands respectively which have different size distributions. The fourth sample is the windblown sands from Leighton Buzzard, UK (SPuL.Buzzard). It is documented as being uniformly graded sand and it has slightly rounded particle shapes. The fifth and last sample is the commercially manufactured glass beads available as “Glass Ballotini” were also used. These were also classed as uniformly graded but having near perfect rounded shapes.

The particle size distributions were analysed using the sieve analysis test as specified in BS 1377-2:1990. The shapes of the particles were observed and quantified using a digital microscope, shown in Figure 1(a), and where the magnified images of the sand particles are also given in Figure 1. A total of 100 randomly picked sand particles from each sample were analysed according to three established shape parameters suggested by Cho et al. [5], Goktepe and Sezer [4] and Tsomokos and Georgiannou [7]. The shape parameters are classified as sphericity (S), roundness (R) and regularity (ρ = [S + R]/2).

Sphericity is defined as a function of diameter of the largest inscribed sphere relative to the diameter of the smallest circumscribed sphere. Roundness is quantified as the average of the inscribed diameter of circles that are

(a)(b)(c)(d)(e)(f)

Figure 1. (a) The Digital microscope scanner. General view of the sand sample seen under the digital microscope for (b) SW; (c) SPuKahang; (d) SPg; (e) SPuL.Buzzard; sand and (f) Glass Ballotini samples.

tangent to the edges of a particle over the diameter of the largest inscribed circle. Figure 2 shows the classification of the particle shapes based on a scale of these parametric values. Each particle was weighed to obtain the mass of individual particles. The manufactured Glass Ballotini was not analysed because it has perfectly spherical (S ≈ 1) and very round (R ≈ 1) particles.

X-ray Diffraction (XRD) and X-ray Fluorescence (XRF) studies showed that all the test samples were of similar mineralogy. This ensured that the test results were only dependent on the particle size and shape, so as to minimise the number of independent variables. Direct shear box tests were used to simulate static shear loading conditions on the sand samples. The dynamic characteristics were determined using cyclic testing with the GDS dynamic triaxial apparatus. To avoid repetition, the relevant choice of experiment in the design of experiments is outlined in Sections 5 and 6 for the static and dynamic test respectively.

3. Particle Size and Shape Distribution Analysis

The conventional size distribution grading parameters are

Figure 2. Particle shape determination (after [5]).

summarised in Table 2. Advanced analysis was done on the size distribution data and cumulative mass percentage frequency distributions were computed at quarter phi intervals. Graphic measures were adopted as explained by Singh and Chaudhri [8], to determine the various distribution parameters. Figure 3 shows the frequency curves by mass of two extreme distributions of SW and Glass Ballotini samples on a phi scale from the data.

Figure 3 further shows that the peak modes are different between the two samples, as the manufactured Glass Ballotini shows a highly concentrated peak value compared to that of the SW sample which has a wider distribution of particle sizes.

Further statistical analysis of the size distributions are summarised in Table 3, including the graphic mean, standard deviation, skewness and kurtosis of the data. It is seen that the mean values are in phi scale. The negative mean phi values found on the SW and Glass Ballotini samples show that the average size of these test particles is gravel like in nature. Other samples show an average particle size corresponding to coarse to very coarse sand.

All the samples except for the Glass Ballotini show positive skews. This explains that most of the size distribution lies to the left side of the mean on the phi scale. The SW sand shows very slight skewness (−0.1 ≤ s < 0.1) indicating an almost symmetrical distribution, as compared to the rest of the samples which has slight skewness (0.1 ≤ s < 0.6). The low negative skew (−0.6 ≤ s < −0.1) values for the Glass Ballotini show that there is a slight distribution which lies to the right side of the mean. The kurtosis shows the pattern of the height and sharpness of the peak of the distribution relative to the rest of the data measured. All the samples except for the Glass Ballotini gave moderate kurtosis values (0.6 ≤ k < 1), indicating a flatter peak near the mode of the distribution.

Table 2. Conventional particle grading properties.

Figure 3. Frequency curves of SW and Glass Ballotini on the phi scale.

Table 3. Statistical analysis on the grading parameters.

Strong kurtosis value (k > 1) of the Glass Ballotini confirms a high peak on the distribution as seen in Figure 3.

Figure 4 illustrates the extreme roundness in the distributions for the SW and SPu(L.Buzzard) sand. The samples have distinctively different peak mode values. Table 4 gives the statistical analysis of the shape parameter data.

The sphericity (S) parameters for all the samples gave negative skew, indicating the size of the values lie to the right side of the mean. This shows that the majority of the particles have high sphericity values. The roundness (R) and regularity (ρ) parameter on the other hand showed positive skews. This indicates that most of the distribution lies to the left side of the mean.

The kurtosis for the S parameter shows that the SW and SPg sample have a strong kurtosis value (k > 1), indicating a sharper peak on the distribution. The small negative kurtosis values found on the SPu(Kahang) and SPu(L.Buzzard) indicates a less distinct peak in the sphericity distribution (−0.6 ≤ k < −0.1). For the R and ρ parameters, the highest value in the kurtosis (k > 1) is evident

Figure 4. Frequency curves of SW and SPu(L.Buzzard) on the roundness parameter.

Table 4. Statistical analysis on the shape parameters.

in the SPg sample. This shows a very distinct peak on the R and ρ distributions. The graphic mean and the moment mean were found to be similar as also observed by Singh and Chaudri [8] in their study.

4. Advanced Analysis of Particle Shapes

4.1. Cluster Analysis

Such a statistical study of particle shapes has not been done in the past. This study focussed on providing a better identification on the classification of the particle shapes. Previous researchers used only the average of the data to represent the shape parameters [5,9,10]. However, there is a possibility that such analysis will have some spatial discontinuities. In this study, statistical solutions of cluster analysis, performance indicators and probability distributions were used to find representative shape parameter values.

From the data collected, the k-means cluster analysis using the Minitab software was used to divide the data of the particle shape parameters and the mass into groups to examine similarities and dissimilarities in the characteristic data observed. k-means cluster analysis was used because of its suitability in catering a reasonably large amount of data and the ability to predetermine the number of clusters required [11]. Four clusters were used to separate the data which represented the extreme values of the particle shape and mass data. It can be seen from Table 5 that the data for each sample are divided into four groups of cluster each having its own centroid. The clusters represent the groups where it shows similarities in the data of the shape parameter and particle mass relationship. The cluster with the most number of data observations would be the most dominant group in the sample data.

4.2. Performance Indicators

Performance indicators (PI) were also used with the Matlab software to accurately judge which type of parent distribution (Normal, Weibull, Log-normal and Exponential) is the most appropriate to represent the shape parameters data. The performance indicators (PI) used in the research was as that suggested by Evans et al. [12] and Kottegoda and Rosso [13]. These consist of two measures that are error measures; mean absolute error (MAE), root mean squared error (RMSE), and the accuracy measures; coefficient of determination (R2) and the prediction accuracy (PA). These performance evaluators measure how close the predicted modelled values are to the observed values [14,15]. For error measures, the smallest error will indicate the best distribution while for accuracy measures, values approaching 1 will indicate the best distribution. Table 6 shows the distribution of the shape parameters for all the samples. The appropriate distribution based on the performance indicators are shaded in grey.

4.3. Probability Determination

Once the best fit distributions were identified by performance indicators, the probability density function (PDF) of the fitted distribution was used to calculate the exceedence, or the probability values of the shape parameters which were equalled or have exceeded the range of the most dominant cluster. Table 7 shows the probability percentage of the extreme values of each cluster according to its shape parameters. The probability percentages indicate that the range of the extreme values of the most dominant cluster appears to have the highest percentage, indicating a high likelihood that this range of values would occur in a set of samples. Hence, the most dominant cluster centroid was used as a representative

Table 5. Summary of the cluster analysis for all the samples of each shape parameter of sphericity (S), roundness (R) and regularity (ρ).

Table 6. The distributions of the shape parameters for all the samples based on the performance indicators.

Table 7. The probability percentage of the extreme values of each cluster according to its shape parameter.

value for the classification of the particle shape parameters. Table 8 summarises the centroids of the most dominant cluster for all the samples.

5. Static Shear Characteristics

The static shear behaviour of the sand samples was studied using the direct shear box apparatus. Three static normal stresses (σN) of 25, 50 and 100 kPa was imposed on all the samples to study the influence of particle morphology on the critical shear stress (τcr). The shearing rate used was kept the same at 0.5 mm/min for all the tests. Figure 5 shows the schematic diagram of the direct shear box test with annotations of the applied forces on the sample. Figure 6 illustrates the typical stress-strain relationship for all the samples where the critical stress occurs at large strains. Figure 7 shows the variation of the uniformity coefficient (Cu) and the critical shear stress parameter (τcrN). A dimensionless relationship was obtained from analysing the critical shear stress (τcr) normalised with respect to the normal stress (σN) against the uniformity coefficient (Cu) for all the samples.

A trendline was obtained that is presented by a single polynomial relationship encompassing a series of different normal stresses. The Cu and σN are known variables which is used to establish an independent variable which enables prediction of the critical shear stress (τcr) of the soil. The best fit curve is represented by a polynomial trendline with a fairly exceptional coefficient of determination, R2 value and is as shown in the Equation (1).

Table 8. The representative values of the shape parameters based on the cluster analysis and probability distribution.

Figure 5. Schematic diagram of the direct shear box test.

ΤcrN vs Cu, (R2 = 0.793):

(1)

Figure 6. Stress-strain relationship from the direct shear box test for all the samples.

Figure 7. The relationship of the normalised shear stress (τcrN) against uniformity coefficient (Cu).

Figure 8 shows the relationship of the effective diameter (D10) with the normalised critical shear stress (τcrN). The D10 relationship with τcrN is inversely proportional to the Cu versus the τcrN relationship where the decreasing values of D10 would also decrease the strength of the soil. Again a curve was obtained that showed a polynomial relationship (Equation (2)).

τcrN vs D10, (R2 = 0.556):

(2)

The representative shape parameters from the cluster analysis were also studied for its relationships with the critical shear stress of the soil. It is seen in Figure 9 that the roundness (R) parameter have better fitting relationships than the sphericity (S) parameter. However, the regularity (ρ) parameter incorporating both the S and R parameter to an average value shows R2 values that represent better relationships of the particle shapes and the critical shear strength. It seems that more spherical and rounder the particles are, the normalised shear stress would decrease to give a logarithmic trendline as presented in the set of equations below.

τcrN vs S, (R2 = 0.462):

(2)

τcrN vs R, (R2 = 0.907):

Figure 8. The relationship of the normalised shear stress (τcrN) against the effective diameter (D10).

Sphericity (S)

Roundness (R)

Regularity (ρ)

Figure 9. The relationship of the normalised shear stress (τcrN) against the (a) sphericity (S), (b) roundness (R) and (c) regularity (ρ).

(3)

τcrN vs ρ, (R2 = 0.755):

(4)

By studying the particle size and shape distribution parameters as well as the static stresses induced, relationships can be made to provide a significant understanding of predicting the critical shear strength behaviour of the soil.

6. Dynamic Loading Mechanical Behavior

Undrained cyclic triaxial tests on fully saturated samples were used to study the influence of particle morphology on the dynamic design parameters which is known as modulus reduction [16] or normalised modulus ratio (G/Gmax) and damping ratio (D). Anisotropic consolidation was used to specifically simulate stress conditions of a soil under asymmetric loading and at shallow depths [17]. This process was done by initially applying an all around confining pressure of 100 kPa, followed by an initial static axial stress (σi) which was always kept at half of the monotonic strength/peak deviator stress (σp) of the sample (σi = 0.5σp).

A cyclic loading frequency of 1.0 Hz and cyclic stress ratio (CSR) of 0.2 was used for the dynamic loading. The sand samples prepared were undercompacted as suggested by Ladd [18] and the size of the samples were 50 mm in diameter and 100 mm in length. Figure 10 shows the final set up of the tiraxial test apparatus where the cyclic loading was applied from the bottom of the sample.

The stress-strain behaviour observed in a typical cyclic triaxial test is represented by a hysteresis loop as shown in Figure 11. The hysteresis loop is an important feature which is used to determine important dynamic parameters such as the modulus reduction (G/Gmax) and damping ratio (D). These parameters were determined as explained by Kramer [16]. Figure 12 shows the effect of the particle size distribution parameter of effective size (D10) on the modulus reduction curve (G/Gmax).

It is seen that an increase in D10 would also show an

Figure 10. The dynamic triaxial test apparatus.

Figure 11. Typical stress-strain relation from a cyclic triaxial test.

Figure 12. The behaviour of the dynamic modulus reduction curve with increasing effective size (D10).

increase in the modulus reduction curve. The shape parameters of sphericity (S), roundness (R) and regularity (ρ) was also seen to affect the modulus reduction curve as pointed out in Figure 13. The decrease in the shape parameters on the other hand increased the modulus reduction curve.

The energy dissipated in the samples during the cyclic test was quantified using the damping ratio (D). As was observed with the G/Gmax behaviour, the morphology characteristics of the sand samples showed a considerable influence on the damping ratio (D). The results indicated that a decrease in the effective size (D10) resulted in an increase in the damping ratio (Figure 14). A decrease in the shape parameters (S, R and ρ) then showed an increased overall damping ratio (Figure 15).

The behaviour of the modulus reduction (G/Gmax) curve was found to be dependent on the effective size (D10) and the change in the shape parameters of sphericity (S), roundness (R) and regularity (ρ). Such variations were observed from past research by Vucetic and Dobry [19], Kramer [16], Okur and Ansal [20] to occur only on soils with high plasticity. However, despite the non existence of plasticity in the samples used in this

Figure 13. The behaviour of the dynamic modulus reduction curve with decreasing particle shape characteristics (S, R, ρ).

Figure 14. The relationships of the damping ratio (D) with increasing cyclic shear strain (γc) and effective size (D10).

Figure 15. The relationships of the damping ratio (D) with increasing cyclic shear strain (γc) and shape parameters (S, R, ρ).

research, there were still some variations in the dynamic shear modulus and modulus reduction curves indicating a high dependency on particle morphology.

7. Conclusions

The influence of particle morphology on the static and dynamic properties of sand was studied in this research. The importance of finding a representative shape parameter from the sand particles is required to have a better understanding on the mechanical properties of the soil based on the shapes of the particles.

In accordance with the statistical analysis done on the distributions of the particle size and shapes, the skewness of the particle size distribution showed that the sand samples are coarse. The shape distribution based on the skewness showed that the majority of the sphericity parameter data for all samples were higher than the mean sphericity. However, the majority of the data on roundness and regularity parameters for all the samples were found to have lower values than the mean.

The cluster analysis was able to determine the most dominant values of the shape parameters of sphericity (S), roundness (R) and regularity (ρ). The type of distribution of the shape parameter data was determined by performance indicators so that probability plots were able to search for the probability of the occurrence for the most dominant cluster. The analysis shows that there is a high percentage or a high likelihood that this range of shape parameter values is within the dominant cluster in a set of a sample. Hence, the most dominant cluster centroid was used as a representative value for the classification of the particle shape parameters.

The study also found that particle morphology has a significant influence on the static shear loading behaviour of the soil. The dimensionless relationships of the critical shear stress (τcr) normalised with respect to the normal stress (σN) against the particle size and representative shape parameters of the sand samples were obtained. These relationships provided a significant understanding in predicting the critical shear strength behaviour of the soil by studying morphological characteristics of the soil as well as the static stress induced.

The dynamic loading characteristics of the sand samples have shown that the modulus reduction curve (G/Gmax) and the damping ratio (D) were dependent on the effective size (D10) and also on the shape parameters of sphericity (R), roundness (R) and regularity (ρ). However observations reported in the past studies [18-20], found this to only occur on soils with high plasticity. Despite the non-existence of plasticity in the samples tested in this research, notable variations still existed in the dynamic shear modulus and modulus reduction curves indicating a significant dependency on particle morphology.

8. Acknowledgements

The authors wish to record their sincere gratitude to UTHM staff for their constant support and guidance in this research.

REFERENCES

  1. A. J. M. S. Lim, “The Influence of Particle Morphology and Structure on the Static and Dynamic Characteristics of Sand,” Ph.D. Manuscript (in preparation), 2013, University Tun Hussien Onn Malaysia.
  2. A.C. McLean and C. D. Gribble, “Geology for Civil Engineers,” E & FN Spon, an imprint of Chapman & Hall, London, 1979, p. 45.
  3. I. Cavarretta, “The Influence of Particle Characteristics on the Engineering Behaviour of Granular Materials,” Ph.D. Dissertation, Imperial College London, London, 2009.
  4. A. B. Goktepe and A. Sezer, “Effect of Particle Shape on Density and Permeability of Sands,” Proceedings of the Institution of Civil Engineers, Geotechnical Engineering, No. GE6, 2010, pp. 307-320.
  5. G. C. Cho, J. Dodds and J. C. Santamarina, “Particle Shape Effects on Packing Density, Stiffness and Strength: Natural and Crushed Sands,” Journal of Geotechnical and Geoenvironmental Engineering, Vol. 132, No. 5, 2006, pp. 591-602. http://dx.doi.org/10.1061/(ASCE)1090-0241(2006)132:5(591)
  6. P. Gou and X. Su, “Shear Strength, Interparticle Locking and Dilatancy of Granular Materials,” Canadian Geotechnical Journal, Vol. 44, 2006, pp. 570-591.
  7. A. Tsomokos and V. N. Georgiannou, “Effect of Grain Shape and Angularity on the Undrained Response of Fine Sands,” Canadian Geotechnical Journal, Vol. 47, No. 5, 2009.
  8. M. Singh and A. R. Chaudri, “Size analysis of the Late Pliocene-Early Pleistocene Upper Siwalak Sediments, Northwestern Himalaya, India,” International Journal of Geoscience, Vol. 4, No. 8, 2013, pp. 1120-1130. http://dx.doi.org/10.4236/ijg.2013.48106
  9. A. K. Alshibi and M. I. Alsaleh, “Characterizing Surface Roughness and Shape of Sands using Digital Microscopy,” Journal of Computing in Civil Engineering, Vol. 18, No. 1, 2004, pp. 36-45. http://dx.doi.org/10.1061/(ASCE)0887-3801(2004)18:1(36)
  10. M. R. Brooke, “The Influence of Grain Shape on Dilatancy,” Ph.D. Thesis, The University of Arizona, 2008.
  11. O. A. Abbas, “Comparisons Between Data Clustering Algorithms,” The International Arab Journal of Information Technology, Vol. 5, No. 3, 2008, pp. 320-325.
  12. M. Evans, N. Hasting and B. Peacock, “Statistical Distributions,” 3th Edition, Wiley, New York, 2000.
  13. N. T. Kottegoda and R. Rosso, “Statistic, Probability for Civil and Environmental Engineers,” McGraw-Hill, Singapore, 1998.
  14. J. A. Rice, “Mathematical Statistics and Data Analysis, Third Addition,” Cengage Learning, Belmont, Brooks/ Cole, 2007.
  15. D. H. Sanders, “Statistics, A First Course,” 5th Edition, Mcgraw-Hill, Inc., Fort Worth, 1995.
  16. S. L. Kramer, “Geotechnical Earthquake Engineering,” Prentice Hall, Upper Saddle River, 1996.
  17. M. T. Yilmaz, O. Pekcan and B. S. Bakir, “Undrained Cyclic Shear and Deformation Behaviour of Silt-Clay Mixtures of Adapazari, Turkey,” Soil Dynamic and Earthquake Engineering, Vol. 24, No. 7, 2004, pp. 497-507. http://dx.doi.org/10.1016/j.soildyn.2004.04.002
  18. R. S. Ladd, “Preparing Test Specimen using Undercompaction,” Geotechnical Testing Journal, Vol. 1, No. 1, 1978, pp. 16-23. http://dx.doi.org/10.1520/GTJ10364J
  19. M. Vucetic and R. Dobry, “Effect of Soil Plasticity in Cyclic Response,” ASCE Journal of Geotechnical Engineering, Vol. 117, No. 1, 1991, pp. 89-107. http://dx.doi.org/10.1061/(ASCE)0733-9410(1991)117:1(89)
  20. D. V. Okur and A. Ansal, “Stress-Strain Characteristics of Fine Grained Soils under Cyclic Loading,” Cyclic Behaviour of Soils and Liquefaction Phenomena, Proceedings of the International Conference, Bochum, 31 March- 2 April 2004, pp. 179-186.