Next Article in Journal
Surface Plasmon Resonance (SPR)-Based Biosensors as Instruments with High Versatility and Sensitivity
Next Article in Special Issue
Structural Operativity Evaluation of Strategic Buildings through Finite Element (FE) Models Validated by Operational Modal Analysis (OMA)
Previous Article in Journal
Inline Inspection with an Industrial Robot (IIIR) for Mass-Customization Production Line
Previous Article in Special Issue
Application and Algorithm of Ground-Penetrating Radar for Plant Root Detection: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A New Approach to Explore the Surface Profile of Clay Soil Using White Light Interferometry

1
School of Civil Engineering, Qingdao University of Technology, Qingdao 266033, China
2
Key Laboratory of Special Area Highway Engineering, Ministry of Education, Chang’an University, Xi’an 710064, China
3
School of Engineering and Built Environment, Griffith University, Gold Coast, QLD 4222, Australia
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(11), 3009; https://doi.org/10.3390/s20113009
Submission received: 7 April 2020 / Revised: 15 May 2020 / Accepted: 23 May 2020 / Published: 26 May 2020
(This article belongs to the Special Issue Sensors for Nondestructive Testing and Evaluation)

Abstract

:
In order to have a better understanding of the real contact area of granular materials, the white light interference method is applied to explore the real surface morphology of clay soils under high stress. Analysis of the surface profile indicates that there exists a support point height z0 with the highest distribution frequency. A concept of a real contact region (from z0 to z0 + d90; d90 represents the particle size corresponding to 90% of the volume fraction) is proposed by combining a surface profile with the particle size distribution of clay soil. It was found that under the compressive stress of 106 MPa–529 MPa, the actual contact area ratio of clay soil varies between 0.375 and 0.431. This demonstrates an increasing trend with the rise of stress. On the contrary, the apparent porosity decreases with an increasing stress, varying between 0.554 and 0.525. In addition, as the compressive stress increases, the cumulative frequency of apparent profile height (from z0d90 to z0 + d90) has a concentrated tendency with a limited value of 0.9.

1. Introduction

The contact characteristics of granular materials are important for prediction of friction between interfaces in geophysical, geological, petroleum, ocean and geotechnical engineering [1,2,3,4,5,6]. Clay particles are widely distributed in the mainland and ocean. In view of the interfacial friction of clay soil, the theory of “adhesion-ploughing” [7] can be used to explain the macroscopic friction force composition with meso-physical quantities. Friction is the resistance encountered when two mutually contacting surfaces are sliding or rolling. There are numerous factors affecting friction, and lots of researchers explain the phenomenon of friction by the early mechanical interlocking theory and molecular attraction theory within a certain range. However, interpreting friction as a single factor has a one-sided effect and does not work well for the hard-soft interface. Bowden and Tabor [7] proposed the adhesion-ploughing theory to lay the theoretical foundation for modern solid friction. Their theory is based on the joint action of mechanism and molecules. With the continuous improvement of the theory [8,9,10,11], it can be well applied to the hard-soft interface.
The study of tribological behavior of granular materials has been remarkably developed in recent years [12,13,14,15,16,17,18,19,20,21,22,23,24,25,26]. There have been theoretical advances, including discrete element method and modern tribology theory. The discrete element method can better simulate the application of granular materials, but the constitutive model relationship between granular materials needs to be further explored [27,28,29,30,31,32]. Interface friction has been studied for a long time and several modern tribology theories have been proposed simultaneously. These theories can well explain the composition of the interface friction formed by a continuum, but they still need further systematic work to better explain the composition of the interface friction formed by discrete particles.
On the other hand, the continuous improvement of meso-testing methods, such as SEM (Scanning Electron Microscope), TEM (Transmission Electron Microscope), AFM (Atomic Force Microscope), and other testing methods [33,34,35,36], can be used to conduct mesoscopic research on granular materials. Numerous scholars [36,37,38,39,40,41,42,43,44,45] applied advanced experimental methods to study granular materials and explore the friction characteristics of granular materials. It has been reported in the literature that the frictional force at the interface of granular materials is generated by ploughing [17], and the actual contact stress between particles is 0.8–2 Gpa. The magnitude of the ploughing force is closely related to the actual contact area. The friction at the interface depends on the number of contact points between the particles, that is, the total contact area of all contact points.
It is obvious that the actual contact area Ar is indispensable in the constitution of macro-friction [46,47,48,49]. In order to explore the actual contact area of granular materials [50,51,52,53,54,55], many researchers have applied the binarization method to the gray image produced by scanning electron microscopy, thereby establishing a 3D model. The real contact area ratio (CAR), which is defined as the ratio between the real contact area and the total interface area, is finally quantitatively calculated. The specified principle of the binarization method is illustrated in Figure 1. However, the arbitrary choice of gray value will cause significant variations of contact area ratio.
For non-contact measurement methods of thixotropic samples, the original appearance of the samples can be better maintained. It is difficult for an optical microscope to form a clear image through the sample, and the white light interference technique [56] can scan the profile height information of the sample surface well. Especially in the medical field, it is applied to in vivo measurement [57], which reflects the non-damaging advantage of non-contact measurement. White light interference technology was proposed in the 1990s [58] and has been widely applied in the fields of machinery [59], medical treatment, chemistry [60], civil engineering [61], and precision manufacturing. However, there are few relevant references in the field of geotechnical engineering.
In order to accurately capture the real contact area, white light interferometry was introduced to explore the surface profile of the clay soil interface under different compressive pressures, using non-contact measurement with nano-precision. To the authors’ best knowledge, this is the first study of its kind with an aim to understand the microscopic changes in the surface contact area of clay soil.
The aim of this study is to explore the real source of interface resistance through the constituent elements. The present paper is organized into the following sections: in Section 2 and Section 3, the adhesion–ploughing friction theory will be outlined and the composition of the interface friction resistance will be discussed. In Section 4 and Section 5, the surface profile of the clay soil under normal load will be tested by the white light interference technique, which also explains real contact area.

2. Adhesion–Ploughing Friction Theory

2.1. Basic Concept

Under loading, the contact point of the two asperities undergoes plastic deformation due to excessive pressures, and the contact points are firmly adhered into one, which is called adhesion. When adhesion occurs, two surfaces of the asperities slide relatively, the adhesion point is sheared, and the shearing force is the adhesion resistance in the frictional resistance.
When the surface of asperities with different hardness contacts and relative sliding occurs, hard asperities will generate ploughing on the surface of soft asperities, and the ploughing resistance at the asperities is another part of the frictional resistance.
The source of frictional resistance [62] is divided into two parts based on the adhesion–ploughing friction theory, namely, adhesion resistance and ploughing resistance. The generation of adhesion resistance mainly depends on the number of contact areas, and the number of contact areas is highly related to the normal load and surface profile. The ploughing resistance mainly depends on the surface profile height of the soft and hard interface. In the final analysis, the height of the surface profile and the magnitude of the normal load are the dominant factors that determine the resistance.

2.2. Surface Interaction

2.2.1. Adhesion Effects

Since the actual contact area only occupies a small percent of the apparent contact area, the stress at the actual contact point reaches the yield strength σs and plastic deformation occurs, and the actual contact point is firmly adhered by the compressive stress. Because the stress of the actual contact points has reached the yield strength, the loading increment can only be borne by increasing the actual contact area (as asperity A,C in Figure 2).
New adhesion points will be continuously formed during the relative sliding of two surfaces regardless of the adhesion points cutting. With the alternation of the old and new adhesion points, total adhesion area of the adhesion points will remain unchanged.

2.2.2. Ploughing Effects

When these two asperities relatively slide, a variety of material migration phenomena may occur at the contact points due to the ploughing, which is related to the hardness of the asperities material (as ploughing asperities A,B in Figure 2). The material migration phenomenon is mainly divided into the following:
  • When asperity B (or A) passes through asperity D, asperity D undergoes relatively large plastic deformation and becomes adhesive. If the adhesive strength between the hard and soft medium is greater than the adhesive strength of the soft medium, sliding shear occurs inside the soft medium, making the soft medium transfer from the lower surface to the upper surface (asperity D is ploughed and adhered to the right side of asperity B, even filling asperity B).
  • When asperity B (or A) passes through asperity D, it undergoes plastic deformation. If the adhesive strength between the hard and soft medium is smaller than the adhesive strength of the soft medium, the adhesion between the soft and hard medium is not firm (asperity D is ploughed and cannot be adhered to B). Asperity B (or A) ploughs along D, which means cutting along the interface. At this time, material migration and deformation (ploughing) instead of transfer occurs at asperity D.
  • Asperity D only produces elastic deformation, and asperity B (or A) slides through D more easily than D produces plastic deformation.

2.3. Adhesion–Ploughing Calculation Formula

Friction is the sum of the resistance created by adhesion effects and ploughing effects. Hence, the friction F is:
F = T + P e = A r · τ b + S · p e
where T and Pe are adhesion resistance and ploughing resistance, respectively; pe is ploughing resistance unit area; τb is shear strength at the adhesion point; S is ploughing contact area; Ar is adhesion contact area.

3. Discretization of Asperities in a Hard–Soft Interface

When considering the hard-soft interface with uniform material properties, it can be assumed that the height of the actual asperities on hard and soft surfaces performs normal distribution function. Therefore, the contacting asperities can be calculated based on the probability. The standard deviation of the asperity height of the hard and soft interface is σ1 and σ2, respectively, and h is the distance between the center lines of the two surfaces’ contact. The contact between two interface asperities can be transferred into the contact between a smooth rigid surface and a two-dimensional normal distribution surface for which standard deviation of the height distribution is σ = σ 1 2 + σ 2 2 . As shown in Figure 3, the outline of asperities is represented by the Greenwood and Williamson random surface model (GW model) [4] with a radius of curvature R and an elastic modulus E’. The problem of the actual contact area of the hard-soft interface is transferred into the height distribution problem of the particles which formed asperities.
The adhesion contact occurs in the portion where the asperities height z > h:
P ( z > h ) = h Ψ z d z
If the total number of asperities is n, then the contacting asperities number m is
m = n h Ψ z d z
The normal deformation of each contact point is z − h. According to the contact mechanics, the contact area of a single asperity is obtained and then multiplied by the contacting asperities number m, obtaining the actual contact area Ar.
A r = m π R z h = n π R h z h Ψ z d z
Converting the contact interface between two random surfaces into a smooth interface and a random surface lays the foundation for the experiment. Using this result, the experiment can use a smooth surface and a random surface to observe the surface profile, reducing a random variable and reducing the test difficulty. It can be seen from Equation (4) that the contact area is related to the surface profile distribution at the height of the interface (z − h).

4. White Light Interferometry Test of Clay Soil

The basic principle of white light interference is to use the same light source to emit two beams of coherent white light at a fixed distance to form interference fringes as the reference optical path. The interference fringes formed on the surface of the measured object due to the difference in profile are used as the detection optical path. By calculating the change of the optical path difference between the reference optical path and the detection light path, the height change of the surface profile is obtained. The white light interferometry system is shown in Figure 4.
The 3D-20 surface profiler from the well-known manufacturer Zeta was used in the test. The measurement accuracy in the XY direction was ≤0.1 μm, and the measurement accuracy in the Z-direction was <20 nm. In order to overcome the difficulty of sample handling and eliminate other interference factors, higher compressive stress was adopted in sample preparation.
The white light interference instrument needs to emit high-frequency white light at each specific height to scan the entire area. After the specific height scan is completed, the next height scan is performed until all the scans within the set scan height range are completed. Due to the fixed focal length, the image with a specific height matching the focal length is clear, and the images outside the specific height are blurred, as shown in Figure 5.

4.1. Sample Preparation

Clay soil used in this experiment was collected from the Jiaozhou Bay of Jiaodong Peninsula, China. The original clay soil possessed a flowing plastic to soft plastic condition with a slight gloss reaction. The specific parameters of the original clay soil are shown in Table 1.
The soil sample after drying and sieving was put into the mold and compressive stress was applied through new polymethyl methacrylate (PMMA) as a smooth and rigid plane. The role of PMMA is the same as that of the top rigid plane shown in Figure 3. The new PMMA used in the test ensured a smooth surface with a diameter of 8.5 mm. In order to make the granular materials have obvious contact, higher compressive stress was set. A total of five samples were made, and the applied normal forces were 2.4 t, 4.8 t, 7.2 t, 9.6 t, and 12 t. The compressive stresses were 106 MPa, 211 MPa, 317 MPa, 423 MPa, and 529 MPa. The whole sample preparation process is shown in Figure 6.

4.2. Experimental Measurement

Specimens were placed in the center of the operation platform (as shown in Figure 7). A proper horizontal and vertical position of the lens was selected to capture the viewing surface. By adjusting the focal length, the sharpness of the image changes. During the process of adjusting the focal length from far to near and then from near to far, the image changes from blurred to clear and then clear to blurred. The critical point where the image changes from blurred to clear is the starting height of the scan, and the critical point that changes from clear to blurred is the ending height of the scan. The number of scanning steps is set according to the required accuracy and the height difference between the upper and lower limits of the scan. The accuracy of each step in this experiment was about 1 μm.

4.3. Particle Size Distribution by Laser Measurement

In order to accurately obtain the soil particle distribution, the clay soil after drying and sieving were subjected to laser particle size analysis. The laser particle size analyzer used in the test was the model Mastersizer 3000, which is measured by the dry method. The analysis results are summarized in Table 2 and the curve of particle size distribution is shown in Figure 8.

4.4. Profile Data Processing

The profile data were processed as follows:
  • A partial section (1500 μm × 1500 μm) in the center of the scanning area (4687 μm × 3516 μm) was selected for analyzing, with a total 45,000 points in this area.
  • Horizontal correction was done by the linear fitting of profile data, thereby avoiding the inclination of contour induced from specimen fabrication and placement. The slope was restricted by 2%, which roughly met the experimental requirements.
  • Statistics analysis on the point cloud data of Z-direction was conducted by taking 5 μm as an interval. A support point height z0 with the highest distribution frequency was identified regardless of different vertical stress.

5. Discussion

5.1. Actual Contact Plane Location

Loose soil particles were compressed under high stress. With an increase of vertical stress, the growth rate of effective stress between soil particles was lower than that of total stress. In order to bear the additional compressive stress ΔP, the height of the outer surface contour decreased and the actual contact area increased ΔS. Since the contact surface was pressed by a smooth PMMA, the height of the actual contact plane contour was unified, as illustrated in Figure 9.
There existed a unified profile height for the bearing surface. In other words, it was the point with the highest distribution frequency, defined as the support point height z0. Particle size distribution was combined in the location of the actual contact surface. The average particle diameters d10, d25, d50, d75, d90 were chosen to calculate the cumulative frequency of apparent profile height under various compression stresses (as shown in Table 3).

5.2. The Relationship between Actual Contact Area and Overburden Stress

As defined above, the support center z0 should be the height of the actual bearing surface, hence, there should be no contour above this height theoretically. However, the scanning results indicate that 96% of profile height is distributed in the range from z0d90 to z0 + d90, rather than a simple cross-sectional area with a certain height. This can be explained by the rebound effect of the soil particles after unloading. Hereby, a concept of real contact region (from z0 to z0 + d90) is proposed in this study, where all the particles undertake the external pressure in a cooperative working mechanism. The actual contact area ratio should be the accumulated frequency from z0 to z0 + d90, in which d90 is the representative value of particle size. The real contact region located in the profile is illustrated in Figure 10, while Figure 11 shows its distribution in the frequency histogram.
Table 4 shows the accumulated frequency of surface profile near the support center. From Figure 12, it can be seen that the actual contact area ratio is relatively stable, changing from 0.375 to 0.431, performing a slightly increasing trend with the rise of overburden stress.
For the sake of making a comparison between the contact area ratio (CAR) produced by SEM and white light interferometry, in this study, the binarization method was also used in analyzing the original SEM image of clay soil under 106 MPa compressive pressure. As shown in Figure 13, the arbitrary choice of gray value can bring about a significant variation of the contact area ratio. More specifically, when the gray value changed from 40 to 80, the contact area ratio decreased drastically from 98.39% to 79.56%. Hereby, it is can be proven the contact area ratio calculated directly from the 3D image of white light interferometry is more reliable than that produced from the binarization of the SEM image.

5.3. Apparent Porosity and Compressive Stress

The pores near the support center change as the stress increases. According to the surface profile data of each point obtained by the white light interference, Equation (5) calculates the difference between the z0 + d90 height and the height of each point, and integrates over the entire range to calculate the porosity. Equation (5) is used to calculate the volume integration in the range of d90, above and below the support center z0. The schematic diagram of apparent porosity calculation is illustrated in Figure 14.
e = z 0 d 90 z 0 + d 90 z 0 + d 90 z d z d x d y z 0 d 90 z 0 + d 90 d z d x d y
where e is the apparent porosity.
The relationship between apparent porosity and compressive stress is illustrated in Figure 13. It shows that the apparent porosity decreased linearly with increasing stress, but at a slower rate as observed through the amplifying relationship.
Table 4 also shows that with an increase of compressive stress, the cumulative frequency of apparent profile height (z0d90z0 + d90) had a concentrated tendency with a limited value, which resulted from the stabilization of porosity. A logistic function was used to derive this limit asymptotic value under the different initial accumulated frequency of apparent profile; the limit asymptotic value was about 0.9. The fitted curve is shown in Figure 15.
The data in Table 5 show that, as the initial accumulated frequency of apparent profile (zero stress or natural accumulation state) increased near the support center, the limit asymptotic value was reduced and maintained at a high level of 0.907–0.921. The fitted curve and corresponding fitted equation are illustrated in Figure 16 and Equation (6), respectively. This observation provides an experimental basis for the PFC (Particle Flow Code) simulation.
y = 0.92126 0.00561 x 0.04357 x 2

6. Conclusions

An innovative mesoscopic study on the surface profile of clay soil was conducted in this study by white light interference. The concept of real contact region (from z0 to z0 + d90) was proposed when considering the rebound of soil particles. It was proven that the proposed new method can accurately reveal the real contact area ratio (CAR), which is indispensable in the modern tribology theory of adhesion-ploughing. The stresses used in the present experiments were all over 100 MPa, and stresses of this magnitude mostly appear in the contact between the tips of the soil particles, and rarely occur in the large-area contact in engineering. The particle contact under stress commonly used in engineering still needs further study. Based on the experiments presented, the following conclusions can be drawn:
  • Under high compressive stress, the real contact area ideally should be the total area of all soil particles in the real contact region (from z0 to z0 + d90), when considering the rebound of the soil particles after unloading, rather than a simple cross-sectional area with a certain height. Hence, compared with the traditional SEM, white light interferometry has an obvious advantage to reveal the real contact area, which would be quite useful in the mesoscopic research of soil particles friction.
  • The actual contact area ratio of clay soil is in the range of 0.375–0.431, and the apparent porosity is approximately 0.525–0.554, demonstrating a linear increase with the rise of overburden pressure.
  • With an increase of compressive stress, the cumulative frequency of apparent profile height (from z0d90 to z0 + d90) has a concentrated tendency with limited value, and a logistic function is used to derive this limit asymptotic value.
  • A regression formula for estimating the limit asymptotic value under various initial values was established for clay soil, which provides an experimental basis for the DEM (Discrete Element Method) simulation.

Author Contributions

The design of the experiments was conducted by M.Z. The execution of the experiments, data analysis, and preparation of the manuscript were conducted by S.Y. under the supervision of L.X., J.L., and D.-S.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 41772318, 51778312), Shandong Key Research and Development Plan (2017GSF20107), Cooperative Innovation Center of Engineering Construction and Safety in Shandong Blue Economic Zone, Open Fund of State key Laboratory of Coastal and Offshore Engineering (LP1712), Ministry of Housing and Urban-Rural Development of China (2014-K3-026).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yang, S.C.; Liu, J.L.; Zhang, M.Y. Analytical solution and field test of critical bearing capacity and settlement of pile tip. Civ. Eng. J. Staveb. Obz. 2020, 1, 61–73. [Google Scholar] [CrossRef]
  2. Yang, S.C.; Zhang, M.Y.; Wang, Y.H.; Sang, S.K.; Miao, D.Z. Field test on pile tip resistance of closed-end jacked pipe pile penetrating into layered foundation. Rock Soil. Mech. 2018, 39, 91–99. [Google Scholar] [CrossRef]
  3. Liu, J.W.; Cui, L.; Zhu, N.; Han, B.; Liu, J. Investigation of cyclic pile-sand interface weakening mechanism based on large-scale CNS cyclic direct shear tests. Ocean Eng. 2019, 194, 11. [Google Scholar] [CrossRef]
  4. Zhu, J.F.; Zhao, H.Y.; Jeng, D.S. Effect of principal stress rotation on dynamic characteristics of a sandy seabed under a partially reflected standing wave. Ocean Eng. 2020, 196, 17. [Google Scholar] [CrossRef]
  5. Zhu, J.F.; Zhao, H.Y.; Jeng, D.S. Effects of principal stress rotation on wave-induced soil response in a poro-elastoplastic sandy seabed. Acta Geotech. 2019, 14, 1717–1739. [Google Scholar] [CrossRef]
  6. Zhu, J.F.; Zhao, H.Y.; Jeng, D.S. Dynamic characteristics of a sandy seabed under storm wave loading considering the effect of principal stress rotation. Eng. Geol. 2019, 259, 15. [Google Scholar] [CrossRef]
  7. Bowden, F.P.; Tabor, D.F. The friction and lubrication of solids. Am. J. Phys. 1951, 19, 428–429. [Google Scholar] [CrossRef]
  8. Greenwood, J.; Williamson, J. Contact of Nominally Flat Surfaces. Proc. R. Soc. Lond. A 1966, 295, 300–319. [Google Scholar] [CrossRef]
  9. Prokopovich, P.; Perni, S. Multiasperity Contact Adhesion Model for Universal Asperity Height and Radius of Curvature Distributions. Langmuir 2010, 26, 17028–17036. [Google Scholar] [CrossRef]
  10. Prokopovich, P.; Perni, S. Comparison of JKR- and DMT-based multi-asperity adhesion model: Theory and experiment. Colloids Surf. A Physicochem. Eng. Asp. 2011, 383, 95–101. [Google Scholar] [CrossRef]
  11. Petean, P.G.C.; Aguiar, M.L. Determining the adhesion force between particles and rough surfaces. Powder Technol. 2015, 274, 67–76. [Google Scholar] [CrossRef]
  12. Hanaor, D.A.H.; Gan, Y.; Einav, I. Static friction at fractal interfaces. Tribol. Int. 2016, 93, 229–238. [Google Scholar] [CrossRef]
  13. Yang, L.; Wang, D.; Guo, Y.; Liu, S. Tribological behaviors of quartz sand particles for hydraulic fracturing. Tribol. Int. 2016, 102, 485–496. [Google Scholar] [CrossRef]
  14. Nardelli, V.; Coop, M.R.; Andrade, J.E.; Paccagnella, F. An experimental investigation of the micromechanics of Eglin sand. Powder Technol. 2017, 312, 166–174. [Google Scholar] [CrossRef]
  15. Sandeep, C.S.; He, H.; Senetakis, K. An experimental micromechanical study of sand grain contacts behavior from different geological environments. Eng. Geol. 2018, 246, 176–186. [Google Scholar] [CrossRef]
  16. Sandeep, C.S.; Senetakis, K. Effect of Young’s Modulus and Surface Roughness on the Inter-Particle Friction of Granular Materials. Materials 2018, 11, 217. [Google Scholar] [CrossRef] [Green Version]
  17. Sandeep, S.; Senetakis, K. Grain-scale mechanics of quartz sand under normal and tangential loading. Tribol. Int. 2018, 117. [Google Scholar] [CrossRef]
  18. He, H.; Sandeep, C.S.; Senetakis, K. The interface behavior of recycled concrete aggregate: A micromechanical grain-scale experimental study. Constr. Build. Mater. 2019, 210, 627–638. [Google Scholar] [CrossRef]
  19. He, H.; Senetakis, K.; Coop, M.R. An investigation of the effect of shearing velocity on the inter-particle behavior of granular and composite materials with a new micromechanical dynamic testing apparatus. Tribol. Int. 2019, 134, 252–263. [Google Scholar] [CrossRef] [Green Version]
  20. Jeong, S.W.; Park, S.S. Effect of the Surface Roughness on the Shear Strength of Granular Materials in Ring Shear Tests. Appl. Sci. 2019, 9, 10. [Google Scholar] [CrossRef] [Green Version]
  21. Nardelli, V.; Coop, M.R. The experimental contact behaviour of natural sands: Normal and tangential loading. Géotechnique 2019, 69, 672–686. [Google Scholar] [CrossRef]
  22. Sandeep, C.S.; Senetakis, K. An experimental investigation of the microslip displacement of geological materials. Comput. Geotech. 2019, 107, 55–67. [Google Scholar] [CrossRef]
  23. Cundall, P.A.; Strack, O.D.L. A discrete numerical model for granular assemblies. Géotechnique 1979, 29, 47–65. [Google Scholar] [CrossRef]
  24. Mo, Y.F.; Turner, K.T.; Szlufarska, I. Friction laws at the nanoscale. Nature 2009, 457, 1116–1119. [Google Scholar] [CrossRef] [PubMed]
  25. Forterre, Y.; Pouliquen, O. Flows of dense granular media. In Annual Review of Fluid Mechanics; Annual Reviews: Palo Alto, CA, USA, 2008; Volume 40, pp. 1–24. [Google Scholar]
  26. Vakis, A.I.; Yastrebov, V.A.; Scheibert, J.; Nicola, L.; Dini, D.; Minfray, C.; Almqvist, A.; Paggi, M.; Lee, S.; Limbert, G.; et al. Modeling and simulation in tribology across scales: An overview. Tribol. Int. 2018, 125, 169–199. [Google Scholar] [CrossRef]
  27. Luding, S. Cohesive, frictional powders: Contact models for tension. Granul. Matter 2008, 10, 235–246. [Google Scholar] [CrossRef] [Green Version]
  28. Wensrich, C.M.; Katterfeld, A. Rolling friction as a technique for modelling particle shape in DEM. Powder Technol. 2012, 217, 409–417. [Google Scholar] [CrossRef]
  29. Lu, G.; Third, J.R.; Muller, C.R. Discrete element models for non-spherical particle systems: From theoretical developments to applications. Chem. Eng. Sci. 2015, 127, 425–465. [Google Scholar] [CrossRef]
  30. Kruggel-Emden, H.; Rickelt, S.; Wirtz, S.; Scherer, V. A study on the validity of the multi-sphere Discrete Element Method. Powder Technol. 2008, 188, 153–165. [Google Scholar] [CrossRef]
  31. Mcdowell, G.; Li, H.; Lowndes, I. The importance of particle shape in discrete-element modelling of particle flow in a chute. Géotechnique Lett. 2011, 1, 59–64. [Google Scholar] [CrossRef]
  32. Liu, J.; Duan, N.; Cui, L.; Zhu, N. DEM investigation of installation responses of jacked open-ended piles. Acta Geotech. 2019, 14, 1805–1819. [Google Scholar] [CrossRef]
  33. Thakur, M.M.; Penumadu, D.; Bauer, C. Capillary Suction Measurements in Granular Materials and Direct Numerical Simulations Using X-Ray Computed Tomography Microstructure. J. Geotech. Geoenvironmental Eng. 2020, 146. [Google Scholar] [CrossRef]
  34. Zhu, J.; Zhu, L.; Zhu, R.; Tian, S.; Li, J. Surface microtopography of surfactant modified montmorillonite. Appl. Clay Sci. 2009, 45, 70–75. [Google Scholar] [CrossRef]
  35. Desbois, G.; Urai, J.L.; Hemes, S.; Brassinnes, S.; De Craen, M.; Sillen, X. Nanometer-scale pore fluid distribution and drying damage in preserved clay cores from Belgian clay formations inferred by BIB-cryo-SEM. Eng. Geol. 2014, 179, 117–131. [Google Scholar] [CrossRef]
  36. Omrani, E.; Tafti, A.P.; Fathi, M.F.; Moghadam, A.D.; Rohatgi, P.; D’Souza, R.M.; Yu, Z.Y. Tribological study in microscale using 3D SEM surface reconstruction. Tribol. Int. 2016, 103, 309–315. [Google Scholar] [CrossRef]
  37. Zhou, W.H.; Jing, X.Y.; Yin, Z.Y.; Geng, X.Y. Effects of particle sphericity and initial fabric on the shearing behavior of soil-rough structural interface. Acta Geotech. 2019, 14, 1699–1716. [Google Scholar] [CrossRef]
  38. Knight, C.; O’Sullivan, C.; van Wachem, B.; Dini, D. Computing drag and interactions between fluid and polydisperse particles in saturated granular materials. Comput. Geotech. 2020, 117. [Google Scholar] [CrossRef]
  39. Hanaor, D.A.H.; Gan, Y.; Einav, I. Effects of surface structure deformation on static friction at fractal interfaces. Géotechnique Lett. 2013, 3, 52–58. [Google Scholar] [CrossRef] [Green Version]
  40. Han, Y.; Abousleiman, Y.N.; Hull, K.L.; Al-Muntasheri, G.A. Numerical Modeling of Elastic Spherical Contact for Mohr-Coulomb Type Failures in Micro-Geomaterials. Exp. Mech. 2017, 57, 1091–1105. [Google Scholar] [CrossRef]
  41. Cavarretta, I.; Rocchi, I.; Coop, M.R. A new interparticle friction apparatus for granular materials. Can. Geotech. J. 2011, 48, 1829–1840. [Google Scholar] [CrossRef]
  42. Aparicio, P.; Galán, E.; Valdrè, G.; Moro, D. Effect of pressure on kaolinite nanomorphology under wet and dry conditions: Correlation with other kaolinite properties. Appl. Clay Sci. 2009, 46, 202–208. [Google Scholar] [CrossRef]
  43. Cavarretta, I.; Coop, M.; O’Sullivan, C. The influence of particle characteristics on the behaviour of coarse grained soils. Géotechnique 2010, 60, 413–423. [Google Scholar] [CrossRef] [Green Version]
  44. Liu, J.; Guo, Z.; Zhu, N.; Zhao, H.; Garg, A.; Xu, L.; Liu, T.; Fu, C. Dynamic Response of Offshore Open-Ended Pile under Lateral Cyclic Loadings. J. Mar. Sci. Eng. 2019, 7, 128. [Google Scholar] [CrossRef] [Green Version]
  45. Liu, J.; Guo, Z.; Han, B. Load Transfer of Offshore Open-Ended Pipe Piles Considering the Effect of Soil Plugging. J. Mar. Sci. Eng. 2019, 7, 313. [Google Scholar] [CrossRef] [Green Version]
  46. Gaboreau, S.; Robinet, J.C.; Pret, D. Optimization of pore-network characterization of a compacted clay material by TEM and FIB/SEM imaging. Microporous Mesoporous Mater. 2016, 224, 116–128. [Google Scholar] [CrossRef]
  47. Matsuda, K.; Hashimoto, D.; Nakamura, K. Real contact area and friction property of rubber with two-dimensional regular wavy surface. Tribol. Int. 2016, 93, 523–529. [Google Scholar] [CrossRef]
  48. Pedrotti, M.; Tarantino, A. Effective stresses for unsaturated states stemming from clay microstructure. Geomech. Energy Environ. 2018, 15, 74–84. [Google Scholar] [CrossRef] [Green Version]
  49. Eyo, E.U.; Ng’ambi, S.; Abbey, S.J. Effect of intrinsic microscopic properties and suction on swell characteristics of compacted expansive clays. Transp. Geotech. 2019, 18, 124–131. [Google Scholar] [CrossRef]
  50. Otsu, N. Threshold Selection Method from Gray-Level Histograms. IEEETrans. Syst. Man Cybern. 1979, 9, 62–66. [Google Scholar] [CrossRef] [Green Version]
  51. Pal, N.R.; Pal, S.K. A Review on Image Segmentation Techniques. Pattern Recognit. 1993, 26, 1277–1294. [Google Scholar] [CrossRef]
  52. Liu, Y.; Srihari, S.N. Document image binarization based on texture features. IEEE Trans. Pattern Anal. Mach. Intell. 1997, 19, 540–544. [Google Scholar] [CrossRef]
  53. Blayvas, I.; Bruckstein, A.; Kimmel, R. Efficient computation of adaptive threshold surfaces for image binarization. Pattern Recognit. 2006, 39, 89–101. [Google Scholar] [CrossRef]
  54. Ramirez-Ortegon, M.A.; Tapia, E.; Ramirez-Ramirez, L.L.; Rojas, R.; Cuevas, E. Transition pixel: A concept for binarization based on edge detection and gray-intensity histograms. Pattern Recognit. 2010, 43, 1233–1243. [Google Scholar] [CrossRef]
  55. Andersen, D.H.; Zhang, Z.L. Contact area on rough surface of nonlinear isotropic brittle materials. Wear 2011, 271, 1017–1028. [Google Scholar] [CrossRef]
  56. Fercher, A.F.; Hu, H.Z.; Vry, U. Rough-Surface Interferometry with A 2-wavelength heterodyne Speckle Interferometer. Appl. Opt. 1985, 24, 2181–2188. [Google Scholar] [CrossRef]
  57. Shekhawat, V.K.; Laurent, M.P.; Muehleman, C.; Wimmer, M.A. Surface topography of viable articular cartilage measured with scanning white light interferometry. Osteoarthr. Cartil. 2009, 17, 1197–1203. [Google Scholar] [CrossRef] [Green Version]
  58. Dresel, T.; Hausler, G.; Venzke, H. 3-Dimensional Sensing of Rough Surfaces by Coherence Radar. Appl. Opt. 1992, 31, 919–925. [Google Scholar] [CrossRef]
  59. Zhang, Y.; Kovalev, A.; Meng, Y. Combined effect of boundary layer formation and surface smoothing on friction and wear rate of lubricated point contacts during normal running-in processes. Friction 2018, 6, 274–288. [Google Scholar] [CrossRef] [Green Version]
  60. Maniscalco, B.; Kaminski, P.M.; Walls, J.M. Thin film thickness measurements using Scanning White Light Interferometry. Thin Solid Film. 2014, 550, 10–16. [Google Scholar] [CrossRef] [Green Version]
  61. Yuan, L.; Zhou, L.; Jin, W.; Chan, C.C. Recent progress of white light interferometric fiberoptic strain sensing techniques. Rev. Sci. Instrum. 2000, 71. [Google Scholar] [CrossRef] [Green Version]
  62. Bowden, F.P.; Moore, A.J.W.; Tabor, D. The Ploughing and Adhesion of Sliding Metals. J. Appl. Phys. 1943, 14, 80–91. [Google Scholar] [CrossRef]
Figure 1. Principle diagram of the binarization method.
Figure 1. Principle diagram of the binarization method.
Sensors 20 03009 g001
Figure 2. Adhesion and ploughing.
Figure 2. Adhesion and ploughing.
Sensors 20 03009 g002
Figure 3. Greenwood and Williamson random surface model.
Figure 3. Greenwood and Williamson random surface model.
Sensors 20 03009 g003
Figure 4. White light interferometry system.
Figure 4. White light interferometry system.
Sensors 20 03009 g004
Figure 5. The scanned image at a specific height.
Figure 5. The scanned image at a specific height.
Sensors 20 03009 g005
Figure 6. Fabrication of clay soil specimen.
Figure 6. Fabrication of clay soil specimen.
Sensors 20 03009 g006
Figure 7. White light interferometry using Zeta 3D surface topography instrument. (a) White light interferometer; (b) Sample scanning.
Figure 7. White light interferometry using Zeta 3D surface topography instrument. (a) White light interferometer; (b) Sample scanning.
Sensors 20 03009 g007
Figure 8. The curve of particle size distribution.
Figure 8. The curve of particle size distribution.
Sensors 20 03009 g008
Figure 9. Illustration of the contact plane under compression. (a) 2D profile; (b) 3D profile.
Figure 9. Illustration of the contact plane under compression. (a) 2D profile; (b) 3D profile.
Sensors 20 03009 g009
Figure 10. Real contact region located in the 2D profile.
Figure 10. Real contact region located in the 2D profile.
Sensors 20 03009 g010
Figure 11. Frequency histogram corresponding to various height in the Z direction. (a) 106 MPa; (b) 211 MPa; (c) 317 MPa; (d) 423 MPa; (e) 529 MPa.
Figure 11. Frequency histogram corresponding to various height in the Z direction. (a) 106 MPa; (b) 211 MPa; (c) 317 MPa; (d) 423 MPa; (e) 529 MPa.
Sensors 20 03009 g011
Figure 12. Relationship between the actual contact area ratio and compression stress.
Figure 12. Relationship between the actual contact area ratio and compression stress.
Sensors 20 03009 g012
Figure 13. Analysis of SEM image by binarization method for the clay soil compacted under 106 MPa vertical pressure. (a) Original SEM image; (b) CAR = 98.39% (Gray value = 40); (c) CAR = 92.57% (Gray value = 60); (d) CAR = 79.56% (Gray value = 80).
Figure 13. Analysis of SEM image by binarization method for the clay soil compacted under 106 MPa vertical pressure. (a) Original SEM image; (b) CAR = 98.39% (Gray value = 40); (c) CAR = 92.57% (Gray value = 60); (d) CAR = 79.56% (Gray value = 80).
Sensors 20 03009 g013
Figure 14. The relationship between apparent porosity and compressive stress.
Figure 14. The relationship between apparent porosity and compressive stress.
Sensors 20 03009 g014
Figure 15. Regression curve of contour distribution near the support center.
Figure 15. Regression curve of contour distribution near the support center.
Sensors 20 03009 g015
Figure 16. Relationship between the initial accumulated frequency of apparent profile and its corresponding limit asymptotic value.
Figure 16. Relationship between the initial accumulated frequency of apparent profile and its corresponding limit asymptotic value.
Sensors 20 03009 g016
Table 1. Physical and mechanical parameters of clay soil sample.
Table 1. Physical and mechanical parameters of clay soil sample.
ds
(g/cm3)
Γ
(kN/cm3)
W
(%)
wL
(%)
wp
(%)
Ip
(%)
Cohesion Force
(kPa)
IFA
(o)
Es1-2
(MPa)
2.7318.034.834.821.213.514.48.63.3
Table 2. Particle size representative value.
Table 2. Particle size representative value.
Leveld10d25d50d75d90
Particle size (μm)10.5522.2142.8265.9586.69
d10 represents the particle size corresponding to 10% of the volume fraction.
Table 3. Cumulative frequency of contact surfaces of each height under high stress of clay soil.
Table 3. Cumulative frequency of contact surfaces of each height under high stress of clay soil.
Pressure (MPa)Cumulative Frequency (CF)
z0d90z0d75z0d50z0d25z0d10z0z0 + d10z0 + d25z0 + d50z0 + d75z0 + d90
1060.2800.3320.4140.4960.5560.6080.6620.75400.8810.9630.983
2110.1390.1930.2960.4020.4770.5420.5910.6710.8070.9160.965
3170.1380.1910.2840.4130.4160.5900.6570.7850.9160.9740.986
4230.1100.1630.2580.3870.4950.5570.6430.7450.8940.9700.988
5290.0880.1360.2340.3710.4790.5520.6520.7540.8910.9640.980
Table 4. Cumulative frequency of surface profile near the support center.
Table 4. Cumulative frequency of surface profile near the support center.
Pressure (MPa)Cumulative Frequency (CF)
(z0 + d90) − z0z0 − (z0d90)(z0 + d90) − (z0d90)
1060.3750.3280.703
2110.4230.4030.826
3170.3960.4520.848
4230.4310.4470.878
5290.4280.4640.892
Table 5. Relationship between initial contour cumulative frequency and limit contour CF.
Table 5. Relationship between initial contour cumulative frequency and limit contour CF.
Initial Contour
(z0 + d90) − (z0d90) CF
R2Limit Profile
(z0 + d90) − (z0d90) CF
0.110.920 ± 0.026
0.20.9990.918 ± 0.025
0.30.9990.916 ± 0.024
0.40.9990.912 ± 0.022
0.50.9980.907 ± 0.020

Share and Cite

MDPI and ACS Style

Yang, S.; Liu, J.; Xu, L.; Zhang, M.; Jeng, D.-S. A New Approach to Explore the Surface Profile of Clay Soil Using White Light Interferometry. Sensors 2020, 20, 3009. https://doi.org/10.3390/s20113009

AMA Style

Yang S, Liu J, Xu L, Zhang M, Jeng D-S. A New Approach to Explore the Surface Profile of Clay Soil Using White Light Interferometry. Sensors. 2020; 20(11):3009. https://doi.org/10.3390/s20113009

Chicago/Turabian Style

Yang, Suchun, Junwei Liu, Longfei Xu, Mingyi Zhang, and Dong-Sheng Jeng. 2020. "A New Approach to Explore the Surface Profile of Clay Soil Using White Light Interferometry" Sensors 20, no. 11: 3009. https://doi.org/10.3390/s20113009

APA Style

Yang, S., Liu, J., Xu, L., Zhang, M., & Jeng, D. -S. (2020). A New Approach to Explore the Surface Profile of Clay Soil Using White Light Interferometry. Sensors, 20(11), 3009. https://doi.org/10.3390/s20113009

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

Article Metrics

Back to TopTop