<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Mater.</journal-id>
<journal-title>Frontiers in Materials</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Mater.</abbrev-journal-title>
<issn pub-type="epub">2296-8016</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1137758</article-id>
<article-id pub-id-type="doi">10.3389/fmats.2023.1137758</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Materials</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Energy multi-scale method to analyze the scale effect of soil particles</article-title>
<alt-title alt-title-type="left-running-head">Chen et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fmats.2023.1137758">10.3389/fmats.2023.1137758</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Chen</surname>
<given-names>Jian</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1983888/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Tong</surname>
<given-names>Huawei</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Yuan</surname>
<given-names>Jie</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Fang</surname>
<given-names>Yingguang</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Huang</surname>
<given-names>Xiaofeng</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>School of Civil Engineering, Guangzhou University</institution>, <addr-line>Guangzhou</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>School of Civil Engineering &#x0026; Transportation, South China University of Technology</institution>, <addr-line>Guangzhou</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1453183/overview">Xianze Cui</ext-link>, China Three Gorges University, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2165926/overview">Mingxun Hou</ext-link>, Guangdong Construction Vocational and Technical College, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2177174/overview">Xiaodan Cai</ext-link>, Huizhou University, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Jie Yuan, <email>Jiey8904@163.com</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Structural Materials, a section of the journal Frontiers in Materials</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>16</day>
<month>03</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>10</volume>
<elocation-id>1137758</elocation-id>
<history>
<date date-type="received">
<day>04</day>
<month>01</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>27</day>
<month>02</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Chen, Tong, Yuan, Fang and Huang.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Chen, Tong, Yuan, Fang and Huang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>As a common geological material, soil is a key aspect of construction engineering. Soil has typical multi-scale characteristics, but current multi-scale methods analyze these characteristics only in regard to geometric space. More exploration of the coupling influence mechanism of the basic properties of particles on the microstructure and macroscopic properties of soil is needed. This study analyzed the influence of geometric scale and mineral composition on the surface energy of particles at the microscopic level for development of the energy multi-scale method. Experiments were performed to determine the influence of mineral composition and particle size on the plasticity index (<italic>I</italic>
<sub>p</sub>) of the soil, and experimental results are discussed and interpreted quantitatively using the energy multi-scale method. The conclusions derived from this work are as follows: 1) the mineral composition and particle size of the soil can cause interface and surface effects; 2) the comprehensive ratio of micro-force to weight (CRFW) of the particles can be determined using the energy multi-scale method and quantitatively reflects the influence of particle size and mineral composition on the microscopic properties of the soil; and 3) the energy multi-scale method explains the mechanism of the plasticity index of soil and has allowed identification of a new division of soil plasticity. When the CRFW was used as the control index, the plasticity index of the three materials was practically the same, even if the mineral composition and particle sizes of the three materials were different.</p>
</abstract>
<kwd-group>
<kwd>scale effect</kwd>
<kwd>microstructure</kwd>
<kwd>energy multi-scale method</kwd>
<kwd>adsorbed water</kwd>
<kwd>plasticity index</kwd>
</kwd-group>
<contract-num rid="cn001">51908151</contract-num>
<contract-num rid="cn002">2022A1515110793</contract-num>
<contract-sponsor id="cn001">National Natural Science Foundation of China<named-content content-type="fundref-id">10.13039/501100001809</named-content>
</contract-sponsor>
<contract-sponsor id="cn002">Basic and Applied Basic Research Foundation of Guangdong Province<named-content content-type="fundref-id">10.13039/501100021171</named-content>
</contract-sponsor>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Soil is one of the most common geological materials and it has clear effects on human life; slope slides, debris flow, and the collapse of building foundations are all due to soil failure (<xref ref-type="bibr" rid="B16">Guzzetti and Peruccacci, 2008</xref>; <xref ref-type="bibr" rid="B36">Pham and Kim, 2018</xref>; <xref ref-type="bibr" rid="B22">Li and Zhao, 2022</xref>; <xref ref-type="bibr" rid="B40">Samprada Pradhan and Toll, 2022</xref>; <xref ref-type="bibr" rid="B23">Liang and Ge, 2023</xref>). Particles in soils span about five orders of magnitude in size, from rubble to clay, and show obvious cross-scale characteristics in geometric space. Soil particles also show mineral composition diversity, with varying ratios that indicate differences in primary and secondary minerals. These particles, with irregular shapes, different sizes, and varying mineral content, accumulate to form the complex spatial structure of soil (<xref ref-type="bibr" rid="B29">Mitchell and Soga, 2005</xref>). During the soil deformation process, particles of different sizes and different mineral composition in different scale-spaces vary in mechanical and deformation performance. Therefore, there is an obvious size effect on soil behavior (<xref ref-type="bibr" rid="B12">Fang and Bo, 2016</xref>; <xref ref-type="bibr" rid="B47">Yang and Liu, 2017</xref>).</p>
<p>A series of studies showed that particles with different size and mineral composition influenced the force transmission and coordination of deformation between particles at the contact surface between particles at the micro-level and demonstrated the properties of soils at the macro-level. <xref ref-type="bibr" rid="B30">Monkul and Ozden (2007)</xref> found that, in the microstructure of fine mixed sand, the interactions between fine and coarse particles affected the overall stress&#x2013;strain behavior of soils. <xref ref-type="bibr" rid="B4">Chang and Yin (2010)</xref> also suggested that the soil was an assembly of particles and that the stress&#x2013;strain relationship of the assembly was determined by integrating the behavior of the interparticle contacts in all orientations. <xref ref-type="bibr" rid="B39">Sabbar and Chegenizadeh (2017)</xref> found that fine soil particles formed a certain structure to reduce the pores between the particles, thereby changing the soil&#x2019;s resistance to liquefaction. <xref ref-type="bibr" rid="B10">Dai and Yang (2015)</xref> revealed that the occurrence of liquefaction and temporary reduction of soil strength was mainly ascribed to loss of support from the contacts between fine particles, as well as those between fine and large particles. <xref ref-type="bibr" rid="B19">Li and Zhang (2016)</xref> used CT images and FEM to demonstrate that the first cracks usually appeared at soil&#x2013;gravel interfaces, which is the weak part of a soil&#x2013;rock mixture. <xref ref-type="bibr" rid="B5">Chen and Fang (2019)</xref> found that the diffuse double layers formed on particle surfaces at the microscopic level affected the effective pore for seepage, thereby leading to non-Darcy flow in some clay soils. <xref ref-type="bibr" rid="B13">Gong and Nie (2019)</xref> analyzed the type of contact between particles, including coordination number and normal contact force, at the microscopic level, to reveal the contribution of each contact type to the residual shear strength of the soil. <xref ref-type="bibr" rid="B34">Ni and Huang (2020)</xref> indicated that, according to the electric double layer theory, mineral particles have an important impact on the properties of soil due to their strong surface electrical characteristics. <xref ref-type="bibr" rid="B37">Phan and Bui (2021)</xref> suggested that the shear resistance of the soil is governed mainly by particle contact type and the current fractions of those types. <xref ref-type="bibr" rid="B52">Zhu and Zhang (2020)</xref> believed that larger particles were more likely to form an occlusal structure and provide higher frictional force, and that sand was more likely to slide when the particle size was smaller. <xref ref-type="bibr" rid="B2">Bai and Zhou (2021)</xref> built a coupled thermo-hydro-mechanical mechanism to explore the soil particle rearrangement of granular thermodynamics. <xref ref-type="bibr" rid="B6">Chen and Tong (2022)</xref> found that the charges on particle surfaces adsorb water into the pores, thereby affecting surface characteristics between the particles and thus the permeability characteristics of the soil. Therefore, it is necessary to establish a multi-scale research method linking macro and micro scales to explore the microcosmic mechanism behind the complex scale phenomenon of soil; this endeavor has become the frontier of current geotechnical field research (<xref ref-type="bibr" rid="B27">Mingjing, 2019</xref>).</p>
<p>Scholars have proposed different multi-scale methods to research the scale phenomenon of soil. First, the strain gradient theory and the couple stress theory (<xref ref-type="bibr" rid="B9">Cosserat, 1909</xref>; <xref ref-type="bibr" rid="B1">Ashby, 1970</xref>) were proposed to solve the problem of the granularity of materials and rotation of particles, as well as the discontinuity in the stress and strain transformation between the contact interface of the particles when entering a certain microscopic space scale. Researchers continued to improve these two theories for use in soil. For example, <xref ref-type="bibr" rid="B31">M&#xfc;hlhaus (1989)</xref> derived a continuous model of regular structures within the framework of the Cosserat theory. Based on this model, the ultimate load problem of tunnel statics was solved by using the finite element method. After homogenization, <xref ref-type="bibr" rid="B32">M&#xfc;hlhaus and Oka (1996)</xref> introduced continuous concepts, such as stress and couple stress, to directly homogenize discrete motion equations to obtain a general continuous model for granular media. <xref ref-type="bibr" rid="B35">Oda and Iwashita (2000)</xref>, through tests on sand, found that couple stresses existed in a shear band in a manner consistent with the change of the particle rotation gradient from negative to positive. The presence of small couple stress still plays an important role in the development of microstructure in shear bands. <xref ref-type="bibr" rid="B45">Wang and Zhang (2003)</xref> theoretically analyzed the stability and instability sliding of an elastic rock mass in a shear zone and introduced the strain gradient plasticity theory, which allowed determination of the theoretical relationship between shear stress and the relative displacement of the upper and lower walls of a shear zone. <xref ref-type="bibr" rid="B26">Lu and Cui (2003)</xref> processed numerical simulation to investigate the influences of shear strain gradient and other factors on the evolution of a shear band. In the framework of a Cosserat continuum, <xref ref-type="bibr" rid="B14">Gudehus and Nuebel (2004)</xref> explored the evolution of a sand shear band trough using a numerical method, and the results showed that shear bands can evolve spontaneously in the interior of a granular body, even under homogeneous stress or displacement boundary conditions. <xref ref-type="bibr" rid="B38">Riahi and Curran (2009)</xref> used the finite element method to implement the adopted Cosserat formulation in buckling analysis of three-dimensional layered continua. Research in this area has continued to increase (<xref ref-type="bibr" rid="B8">Chu and Yu, 2012</xref>; <xref ref-type="bibr" rid="B41">Tang and Song, 2015</xref>; <xref ref-type="bibr" rid="B42">Tang and Wei, 2020</xref>; <xref ref-type="bibr" rid="B43">Tang and Wei, 2021</xref>; <xref ref-type="bibr" rid="B3">Borst and Alizadeh, 2022</xref>). The couple stress and strain gradient theories may be used to solve some problems with scale effects of the soil caused by the change in research scale, but these theories are based on the continuum; therefore, the anisotropy problem caused by granularity and mineral composition are not solvable using the same methods.</p>
<p>The discrete element method (DEM) takes the soil particle as the basic research unit, clearly describes the particle properties of the soil, and describes the macroscopic properties of the soil through the interaction of particle groups (<xref ref-type="bibr" rid="B44">Thornton, 2000</xref>). This method has become popular for studying multi-scale phenomena of soil. <xref ref-type="bibr" rid="B49">Zeghal and Shamy (2008)</xref> analyzed the liquefaction of saturated soil particles in loosely cemented sediments through the coupling of discrete elements and continuous elements. <xref ref-type="bibr" rid="B21">Li and Wan (2011)</xref> combined the DEM with the finite element method (FEM) and proposed a method for multi-scale analysis of granular materials at two macroscopic levels. <xref ref-type="bibr" rid="B33">Nguyen and Combe (2014)</xref> proposed a multi-scale modeling method for cohesive granular materials; at the micro-level, the DEM was used to simulate dense particle accumulation, and at the macro level, the numerical solution was obtained through the FEM. <xref ref-type="bibr" rid="B51">Zhao and Feng (2017)</xref> conducted a 3D numerical study on the formation of landslide dams in open fluid channels by using the DEM combined with computational fluid dynamics. Additionally, <xref ref-type="bibr" rid="B24">Liu and Zhixiang (2021)</xref> constructed a numerical model capable of reproducing the estimation of the impact of debris flow by coupling FEM and DEM, <xref ref-type="bibr" rid="B50">Zhang and Zhang (2022)</xref> conducted DEM simulations to quantify the relationship between fine particle loss and upscaling from particle to cell and then revealed the multi-scale response of seepage erosion in sandy strata, and <xref ref-type="bibr" rid="B11">Fang and Yang (2023)</xref> analyzed the progressive failure process of the stratified rock mass through the analysis of contact fracture events and the evolution of contact forces using the DEM-FEM numerical method. DEM and the DEM-FEM coupling method can intuitively display the particle properties of soil. However, limited by current computing capability, the scale of DEM or DEM-FEM models is generally small (<xref ref-type="bibr" rid="B48">Yang and Wang, 2023</xref>) and the particles in the models are greatly simplified. For example, for fine-grained sand&#x2013;gravel mixed soil, the DEM method only models large particles and simplifies fine particles into contact parameters between particles to reduce the amount of calculation required. Clearly, this method cannot reflect the effect of gradation of particles on soil properties. In addition, due to the unclear mechanism of the coupling effects of factors such as particle size and mineral composition on the particle contact surface, the various mineral compositions and the complex contact properties between the particles result in cumbersome calibration of contact parameters in the DEM that ultimately lacks credibility.</p>
<p>The multi-scale phenomenon of soil has long been known, and a variety of multi-scale analysis methods have been developed. However, there are still defects in these methods when used for solving some multi-scale problems with soil. Regardless of the kind of model that is constructed to connect the macro and micro, the transmission of force between particles, the coordination of strain, and even the construction of soil microstructure are realized through the contact surface. The impact of various characteristics of soil particles at the microscopic level on soil engineering properties should be directly reflected on the contact surface between particles, and the surface energy of soil particles determines the surface of particles and the contact surface between particles. Therefore, in this study, the comprehensive ratio of micro-force to weight was proposed to quantitatively analyze the coupled influence of size and mineral composition of particles on the surface energy of soil particles. We used this energy multi-scale method of analysis to reveal the mechanism of the soil plasticity index. Our results are significant as they support accurate calibration of contact parameters of the DEM method and in promoting the application of the DEM method in soil.</p>
</sec>
<sec id="s2">
<title>2 Energy multi-scale method</title>
<sec id="s2-1">
<title>2.1 Energy multi-scale property of soil</title>
<p>Materials are not infinitely separable, and the intrinsic properties of a particle in the internal length range will affect the mechanical properties of the microstructure. Interface effects and surface effects are ignored at the macroscopic level, but they play important roles at the microscopic level (<xref ref-type="bibr" rid="B7">Chiu and Chen, 2011</xref>; <xref ref-type="bibr" rid="B20">Li and Mi, 2019</xref>; <xref ref-type="bibr" rid="B25">Lu and Guo, 2019</xref>; <xref ref-type="bibr" rid="B46">Wang and Zhang, 2019</xref>). Soil shows a typical multi-scale property; that is, the stress and deformation of soil are different in each scale space. A continuous soil slice at the macroscopic level is shown in <xref ref-type="fig" rid="F1">Figure 1A</xref>, which demonstrates that the stress distribution in the soil slice is nearly uniform under the load. However, at the microscopic level, where the research focus is a single soil particle, contact surface discontinuity causes the stress to appear abruptly on the contact surfaces of the particles. In <xref ref-type="fig" rid="F1">Figure 1</xref>B, the cross-scale deformation behavior of the foundation from the macro-level to the micro-level is illustrated. At the macroscopic level, the foundation produced continuous compressive deformation under uniform load, whereas at the microscopic level, the soil particles showed an obvious rotational effect. Microscopically, the particles appeared in individual forms and the contact between particles was highlighted. However, macroscopically, the number of particles in a soil slice is huge and the individuality of particles is weakened; therefore, particle groups are treated as a continuous medium and the contact between particles is no longer considered. In general, the effects of interface properties between particles on soil properties change at different scales.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Scale effect of soil.</p>
</caption>
<graphic xlink:href="fmats-10-1137758-g001.tif"/>
</fig>
<p>For soil particles, when particle size decreases, the specific surface area of particles and the proportion of charges on the surface of particles increase gradually; thus, the ability to adsorb ions and water molecules increases. As a result, the adsorbed water and the electric double layer formed on the surface of the particles are more obvious, which in turn affects the contact surface between the particles. At the same time, a typical soil mixture contains clay minerals and non-clay minerals. Most clay minerals have excess negative charges on the surface of particles due to isomorphous replacement. The negative charges adsorb the polar water molecules and cations to form an electric double layer (<xref ref-type="fig" rid="F2">Figure 2</xref>). Clay particles tend to be smaller in size and have a larger specific surface area; therefore, the volume of adsorbed water and the electric double layer on the surface of particles is considerably larger than the volume of the particles. When clay particles are in contact, the contacts occur mainly between the water films on the particle&#x2019;s surface. For non-clay minerals, the ability to adsorb ions and water molecules is very weak because the excess charges on the surface of the non-clay minerals are less; the adsorbed water on the surface of such mineral particles is so thin that the contact between particles is mainly direct contact between mineral crystals. For soil materials, both particle size and mineral composition influence the properties of the particle surface, thereby affecting the contact surface between particles, as shown in <xref ref-type="fig" rid="F3">Figure 3</xref>. Therefore, as a granular medium, soil particles also exhibit the surface effects and interface effects.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Diffuse electric double layer.</p>
</caption>
<graphic xlink:href="fmats-10-1137758-g002.tif"/>
</fig>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Contact states of different particles.</p>
</caption>
<graphic xlink:href="fmats-10-1137758-g003.tif"/>
</fig>
<p>The surface and interface effects of soil particles affect the soil microstructure (<xref ref-type="bibr" rid="B15">Guo and Yu, 2019</xref>) and thereby affect the soil engineering properties at the macroscopic level, which is the basis for the soil-scale phenomenon. For example, particle size affects the occlusion between sand particles; increased coarse particle content leads to increased microstructure rolling and increased internal friction angle of the soil at the macroscopic level; the electric double layer (DEL) and the adsorbed water on the surface of clay particles result in a cohesive effect among the particles with increased cohesive force and increased clay mineral content; the balance between the microscopic force and the gravity of the particles affects the formation of the microscopic and macroscopic structure of the soil; at the same time, the mineral components have different adsorption water states on the surface, which will also directly affect the state of the liquid flowing over its surface, thereby affecting the permeability characteristics of the soil.</p>
</sec>
<sec id="s2-2">
<title>2.2 Energy multi-scale equation for soil particles</title>
<p>The surface state of a particle is largely determined by the particle&#x2019;s ability to adsorb water or other ions, and this ability is called the particle&#x2019;s surface energy. This adsorption energy is mainly manifested as the van der Waals force between particles and the electric field force generated by the excess charges on the particle surface. At the micro-level, the forces with high size factor power, such as gravity and inertial force, are relatively small; while the forces with low size factor power, such as surface tension, van der Waals force, and electric field force, are relatively large, as shown in <xref ref-type="fig" rid="F4">Figure 4</xref>. The surface and contact interface features of the particles are produced mainly by the direct action of these forces. Among these forces, gravity, van der Waals force, and Coulomb force can be used to characterize both the particle size and the mineral composition of soil particles.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Microscopic force between particles.</p>
</caption>
<graphic xlink:href="fmats-10-1137758-g004.tif"/>
</fig>
<p>In this study, a relationship between gravity, van der Waals force, and Coulomb force was established. As introduced previously, these forces can be used to characterize the adsorption energy of the particle surface. Therefore, our energy multi-scale method involves calculation of the ratio of the micro-Coulomb force to gravity and the ratio of the micro-van der Waals force to gravity. The equations used for these calculations are expressed as follows:<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>G</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
<mml:msub>
<mml:mi>&#x3b5;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mfrac>
<mml:mo>/</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3c1;</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>&#x3c0;</mml:mi>
<mml:msup>
<mml:mi>d</mml:mi>
<mml:mn>3</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mn>6</mml:mn>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>3</mml:mn>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msup>
<mml:mi>&#x3c0;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msub>
<mml:mi>&#x3b5;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:msup>
<mml:mi>r</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mi>&#x3c1;</mml:mi>
<mml:mi>g</mml:mi>
<mml:msup>
<mml:mi>d</mml:mi>
<mml:mn>3</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
<disp-formula id="e2">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>G</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>24</mml:mn>
<mml:msup>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mo>/</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3c1;</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>&#x3c0;</mml:mi>
<mml:msup>
<mml:mi>d</mml:mi>
<mml:mn>3</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mn>6</mml:mn>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>&#x3c1;</mml:mi>
<mml:mi>g</mml:mi>
<mml:mi>&#x3c0;</mml:mi>
<mml:msup>
<mml:mi>d</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:msup>
<mml:mi>H</mml:mi>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
</p>
<p>where <italic>&#x3bb;</italic>
<sub>e</sub> and <italic>&#x3bb;</italic>
<sub>w</sub> are the ratio of micro-Coulomb force to gravity and the ratio of micro-van der Waals force to gravity, respectively, <italic>F</italic>
<sub>e</sub> is Coulomb force, <italic>F</italic>
<sub>w</sub> is van der Waals force, <italic>G</italic> is gravity, <italic>&#x3b5;</italic>
<sub>0</sub> is vacuum permittivity, <italic>q</italic>
<sub>1</sub> and <italic>q</italic>
<sub>2</sub> are particle charge, <italic>r</italic> is particle center spacing, <italic>&#x3c1;</italic> is the particle density, <italic>d</italic> is the particle diameter, <italic>A</italic> is the Hamaker constant, and <italic>H</italic> is the clear particle spacing. Soil particles exist mainly in a water environment, and the amount of charge associated with the particles can be obtained according to colloid and interface chemistry and electrostatic theory (<xref ref-type="bibr" rid="B17">Helmy, 1998</xref>; <xref ref-type="bibr" rid="B28">Missana and Adell, 2000</xref>):<disp-formula id="e3">
<mml:math id="m3">
<mml:mrow>
<mml:mi>q</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
<mml:msub>
<mml:mi>&#x3b5;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:msub>
<mml:mi>&#x3d5;</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
</p>
<p>where <italic>&#x3c6;</italic>
<sub>0</sub> is the particle surface potential and the clay particle surface potential can be determined by the cation exchange capacity (CEC) method. For the convenience of calculation and the unification of measurement indicators, the adjustment coefficient <italic>&#x3b1;</italic> is used to combine <italic>&#x3bb;</italic>
<sub>e</sub> and <italic>&#x3bb;</italic>
<sub>w</sub> as one parameter, and this parameter is called the comprehensive ratio of micro-force to weight (CRFW). The expression of <italic>&#x3bb;</italic> is as follows:<disp-formula id="e4">
<mml:math id="m4">
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msub>
<mml:mi>&#x3bb;</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>e</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>G</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>w</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi>G</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
<p>The adjustment coefficient <italic>&#x3b1;</italic> must be determined according to the properties of the material. Taking the common soil minerals montmorillonite and quartz as examples, the CRFW of these two mineral particles were calculated. When the value of the adjustment coefficient <italic>&#x3b1;</italic> was varied, the CRFW of montmorillonite and quartz particles also changed according to particle size; these data are displayed in <xref ref-type="fig" rid="F5">Figure 5</xref>. As shown in the figure, when <italic>&#x03B1;</italic> is 0 the <italic>&#x03BB;</italic> degenerates as <italic>&#x03BB;<sub>w</sub>
</italic> and when <italic>&#x03B1;</italic> is 1 the <italic>&#x03BB;</italic> degenerates as <italic>&#x03BB;<sub>e</sub>
</italic>, respectively.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Relationship of CRFW with particle diameter.</p>
</caption>
<graphic xlink:href="fmats-10-1137758-g005.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F5">Figure 5</xref> demonstrates that the relationships between <italic>&#x3bb;</italic> and the particle sizes of quartz and montmorillonite were basically the same, and when particle size decreased, <italic>&#x3bb;</italic> increased gradually. The specific surface area of particles is relative to particle size. As a result, the amount of negative charge on the particle surface increased and the surface energy increased; therefore, the thickness of the water film on the particle surface became larger and the interface between particles was mainly soft contact between water layers. In contrast, with increased particle size, the effect of adsorbed water on the particle surface decreased, and the contact between particles was mainly hard contact between mineral crystals. It is reasonable to use the CRFW to describe the variation of soil properties with particle size. At the same time, it can be seen from <xref ref-type="fig" rid="F5">Figure 5</xref> that, although the particle size of two different mineral particles is the same, their CRFW values can be very different. As demonstrated previously, compared with non-clay minerals, the isomorphic replacement of clay minerals is more obvious, and the negative charge density of particle surfaces is higher, so the surface effect is more obvious. This also shows the rationality of the formula for CRFW. In summary, the CRFW in the energy multi-scale method reflects the mechanism of property changes based on particle size and mineral composition at the micro-level.</p>
</sec>
</sec>
<sec id="s3">
<title>3 The plasticity index of soil explained by the energy multi-scale method</title>
<sec id="s3-1">
<title>3.1 Experimental plan</title>
<p>The plasticity index of soil is important for engineering endeavors involving soils such as clay and silty clay. A large number of studies have found that the shear strength, permeability coefficient, and liquefaction index of plastic soils were directly or indirectly related to the plasticity index of soil. Furthermore, the plasticity index of plastic soil reflects the water adsorption capacity of the surface of the particles, which is a direct reflection of the surface energy of the particles. At the same time, the mineral composition and the diameter of the soil particle are the main factors in the calculation of the CRFW by the energy multi-scale method. Therefore, in this study, the surface potential and the diameter of the particles were tested first. Next, the samples with different mineral compositions and relative content, as well as different particle sizes, were prepared to assess the plasticity index of the soil. In terms of soil mineral composition, three typical clay minerals, montmorillonite, kaolinite, and illite, and two non-clay minerals, quartz and feldspar, were selected for these experiments. Due to the high purity of commercial soil minerals, customized commercial soil mineral powders were selected as test material. Samples with clay mineral content of 100%, 60%, 40%, 20%, and 0% were assessed, and the diameter of all mineral particles was 2&#xa0;<italic>&#x3bc;</italic>m. Sample characteristics are summarized in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Sample characteristics.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Sample No.</th>
<th align="center">Composition</th>
<th align="center">Sample No.</th>
<th align="center">Composition</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">YSX-1</td>
<td align="center">100%, montmorillonite (M)</td>
<td align="center">YSX-8</td>
<td align="center">M (60%) &#x2b; Q (20%) &#x2b; F (20%)</td>
</tr>
<tr>
<td align="center">YSX-2</td>
<td align="center">100%, illite (I)</td>
<td align="center">YSX-9</td>
<td align="center">I (20%) &#x2b; Q (40%) &#x2b; F (40%)</td>
</tr>
<tr>
<td align="center">YSX-3</td>
<td align="center">100%, kaolinite (K)</td>
<td align="center">YSX-10</td>
<td align="center">I (40%) &#x2b; Q (30%) &#x2b; F (30%)</td>
</tr>
<tr>
<td align="center">YSX-4</td>
<td align="center">100%, quartz (Q)</td>
<td align="center">YSX-11</td>
<td align="center">M (60%) &#x2b; Q (20%) &#x2b; F (20%)</td>
</tr>
<tr>
<td align="center">YSX-5</td>
<td align="center">100%, feldspar (F)</td>
<td align="center">YSX-12</td>
<td align="center">K (20%) &#x2b; Q (40%) &#x2b; F (40%)</td>
</tr>
<tr>
<td align="center">YSX-6</td>
<td align="center">M (20%) &#x2b; Q (40%) &#x2b; F (40%)</td>
<td align="center">YSX-13</td>
<td align="center">K (40%) &#x2b; Q (30%) &#x2b; F (30%)</td>
</tr>
<tr>
<td align="center">YSX-7</td>
<td align="center">M (40%) &#x2b; Q (30%) &#x2b; F (30%)</td>
<td align="center">YSX-14</td>
<td align="center">K (60%) &#x2b; Q (20%) &#x2b; F (20%)</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The dried mineral powders were mixed evenly according to the designed mass ratio, then a small amount of water was added to wet the mineral particles of the soil; the samples were left to stand overnight. Finally, the liquid and plastic limits of the aforementioned soil samples were measured using a liquid&#x2013;plastic limit tester, and the plasticity index of each group of samples was calculated using the liquid and plastic limit data.</p>
</sec>
<sec id="s3-2">
<title>3.2 Results and discussion</title>
<p>After the liquid&#x2013;plastic test, the plasticity index (<italic>I</italic>
<sub>p</sub>) of each sample was assessed, as shown in <xref ref-type="fig" rid="F6">Figure 6</xref>; the <italic>I</italic>
<sub>p</sub> of all soil samples increased with increased clay mineral content, but the three clay mineral samples had substantially different <italic>I</italic>
<sub>p</sub> values despite having the same particle size and content.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Influence of soil mineral content on sample plasticity index.</p>
</caption>
<graphic xlink:href="fmats-10-1137758-g006.tif"/>
</fig>
<p>According to the energy multi-scale method, the ratio of micro-Coulomb force to gravity and the ratio of micro-van der Waals force to gravity can be calculated using equations <xref ref-type="disp-formula" rid="e1">1</xref> and <xref ref-type="disp-formula" rid="e2">2</xref> and the elements from the literature (<xref ref-type="bibr" rid="B18">Huang, 2018</xref>). When <italic>d</italic> &#x3d; 2 <italic>&#x3bc;</italic>m, <italic>&#x3bb;</italic>
<sub>e,</sub> and <italic>&#x3bb;</italic>
<sub>w</sub> of the three materials could be calculated and is shown in <xref ref-type="table" rid="T2">Table 2</xref>; the <italic>&#x3bb;</italic>
<sub>e</sub> values of montmorillonite and illite were 6.91 and 2.1 times the value of kaolinite, respectively, whereas the <italic>&#x3bb;</italic>
<sub>w</sub> values of montmorillonite and illite were 0.73 and 0.79 times that of kaolinite, respectively, and the <italic>I</italic>
<sub>P</sub> values of montmorillonite and illite were 4.8 and 1.68 times that of kaolinite, respectively. The difference in <italic>&#x3bb;</italic>
<sub>w</sub> between these minerals was not large. All samples were composed of clay minerals; therefore, the adjustment coefficient <italic>&#x3b1;</italic> could be set to 1 for the test. As showed in <xref ref-type="fig" rid="F5">Figure 5</xref>, the <italic>&#x03BB;<sub>e</sub>
</italic> is bigger than <italic>&#x03BB;<sub>w</sub>
</italic> for montmorillonite , while <italic>&#x03BB;<sub>e</sub>
</italic> is samller than <italic>&#x03BB;<sub>w</sub>
</italic> for Quartz ,within a certain particle size range. These data indicate that, in clay, the Coulomb force of particles has a greater influence on the particle surface adsorption energy, and the contribution value of <italic>&#x3bb;</italic>
<sub>e</sub> to <italic>&#x3bb;</italic> is larger for clay.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>CRFW and <italic>I</italic>
<sub>p</sub> of three types of mineral particles (<italic>d</italic> &#x3d; 2&#xa0;<italic>&#x3bc;</italic>m).</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Type</th>
<th align="center">Diameter (<italic>&#x3bc;</italic>m)</th>
<th align="center">
<italic>&#x3bb;</italic>
<sub>e</sub>
</th>
<th align="center">
<italic>&#x3bb;w</italic>
</th>
<th align="center">
<italic>&#x3bb;</italic>
<sub>e</sub>/<italic>&#x3bb;</italic>
<sub>e</sub> (K)</th>
<th align="center">
<italic>&#x3bb;w</italic>/<italic>&#x3bb;w</italic> (K)</th>
<th align="center">
<italic>I</italic>p/<italic>I</italic>p (K)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Montmorillonite</td>
<td align="center">2</td>
<td align="center">53.8</td>
<td align="center">6.89</td>
<td align="center">6.91</td>
<td align="center">0.73</td>
<td align="center">4.8</td>
</tr>
<tr>
<td align="center">Illite</td>
<td align="center">2</td>
<td align="center">17.2</td>
<td align="center">7.54</td>
<td align="center">2.21</td>
<td align="center">0.79</td>
<td align="center">1.68</td>
</tr>
<tr>
<td align="center">Kaolinite</td>
<td align="center">2</td>
<td align="center">7.79</td>
<td align="center">9.49</td>
<td align="center">1</td>
<td align="center">1</td>
<td align="center">1</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>As described previously, the surface state and interface state are the basis of stress transfer and deformation coordination, and the surface state and contact state of the particles at the microscopic level determine the internal structure and engineering properties of the soil. The adsorbed water and the electric double layer are mainly caused by the negative charges on the surface of the particles. The more charges on the particle surface, the stronger the surface adsorption energy of the particle, and the relatively thicker the adsorbed water layer. For clay minerals, the charge density of the particles is high, and the particle size is small, so the thickness of the adsorbed water layer on the surface of the particles is relatively large for a particle with such a tiny diameter. Currently, the properties of the interface between particles are generally determined by the adsorbed water layer. Therefore, many engineering properties of clay soils are affected by the adsorbed water layer. The plasticity index of clay minerals in this study was typical, and the plasticity index refers to the relative water content of plastic soil from the plastic state to the flowing state. It is an important index for measurement of surface adsorption capacity of mineral particles to free water in the soil, and it is also a typical measure of performance of the microscopic adsorption energy of mineral particles at the macroscopic level. The surface charge density of clay mineral particles such as montmorillonite is largest, illite has the second largest, and kaolinite has the smallest, so the plasticity index of each of the three minerals decreases accordingly, as shown in <xref ref-type="fig" rid="F6">Figure 6</xref>. As shown in <xref ref-type="table" rid="T2">Table 2</xref>, <italic>&#x3bb;</italic>
<sub>e</sub> of montmorillonite is the largest, that of illite is the second largest, and that of kaolinite is the smallest. These data indicate that the energy multi-scale method quantitatively reflects the adsorption capacity of the microscopic surface of particles. The adsorbed water layer discussed in this paper is relatively thick compared to the particle diameter. When the particle size is large, even if the adsorbed water layer on the surface of the particle is of the same thickness as that of a smaller particle, the relative thickness can be easily ignored compared to the huge particle size. Therefore, using the energy multi-scale method, both <italic>&#x3bb;</italic>
<sub>e</sub> and <italic>&#x3bb;</italic>
<sub>w</sub> decreased with increased particle size, as shown in <xref ref-type="fig" rid="F5">Figure 5</xref>. In addition, we assessed the plasticity index of each of the three clay minerals when the <italic>&#x3bb;</italic>
<sub>e</sub> parameter was controlled. First, according to Eq. <xref ref-type="disp-formula" rid="e1">1</xref> and Eq. <xref ref-type="disp-formula" rid="e2">2</xref>, the relationship curves between the <italic>&#x3bb;</italic>
<sub>e</sub> and particle size of each of the three materials were calculated, as shown in <xref ref-type="fig" rid="F7">Figure 7</xref>.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Relationship of <italic>&#x3bb;</italic>
<sub>e</sub> with particle size.</p>
</caption>
<graphic xlink:href="fmats-10-1137758-g007.tif"/>
</fig>
<p>The particle diameters of montmorillonite, illite, and kaolinite were set to 20&#xa0;<italic>&#x3bc;</italic>m, 14 <italic>&#x3bc;</italic>m, and 8&#xa0;<italic>&#x3bc;</italic>m, respectively, and the <italic>&#x3bb;</italic>
<sub>e</sub> of each of the three materials were caculated from the fitted equation in <xref ref-type="fig" rid="F7">Figure 7</xref> and the results were shoud in <xref ref-type="table" rid="T3">Table 3</xref>.</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>
<italic>&#x3bb;</italic>
<sub>e</sub> of mineral particles with different diameters.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Type</th>
<th align="center">Montmorillonite</th>
<th align="center">Illite</th>
<th align="center">Kaolinite</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Diameter (<italic>&#x3bc;</italic>m)</td>
<td align="center">20</td>
<td align="center">14</td>
<td align="center">8</td>
</tr>
<tr>
<td align="center">
<italic>&#x3bb;</italic>
<sub>e</sub>
</td>
<td align="center">0.0557</td>
<td align="center">0.0525</td>
<td align="center">0.0558</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>
<xref ref-type="fig" rid="F8">Figure 8</xref> indicates that the plasticity indices of these three clay minerals, calculated using the energy multi-scale method, are very close to each other and nearly on the same straight line. Therefore, there are benefits to use of this method to classify clay soil samples, as the parameter <italic>&#x3bb;</italic> in the energy multi-scale method accounts for the influence of mineral composition and particle size on the microscopic surface and interface states of the particles.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Plasticity index of soil samples with different mineral contents using the energy multi-scale method.</p>
</caption>
<graphic xlink:href="fmats-10-1137758-g008.tif"/>
</fig>
</sec>
</sec>
<sec sec-type="conclusion" id="s4">
<title>4 Conclusion</title>
<p>Soil has an obvious scale effect, and the microscopic characteristics of particles affect the macroscopic properties of soil. This study proposed an energy multi-scale method based on the microscopic properties of particles and established a formula for calculating the CRFW. Additionally, the mechanism of clay plasticity index change was quantitatively explained by the CRFW. The conclusions of this study are as follows:<list list-type="simple">
<list-item>
<p>(1) Both mineral composition and particle size can influence the surface energy of particles and cause interface and surface effects. Higher clay content is associated with greater surface charge of the soil, which results in greater surface adsorption energy and greater relative thickness of the adsorbed water layer on the particle surface. At the same time, when soil particle size is very small, the specific surface area of the particles is large and the thickness of the adsorbed water layer is large. At the microscopic level, the interaction between particles occurs primarily at the contact interface, so these two factors affect the surface characteristics of the particle and thereby affect macroscopic soil properties.</p>
</list-item>
<list-item>
<p>(2) The CRFW of the particles calculated using the energy multi-scale method quantitatively reflects the influence of particle size and mineral composition on the microscopic properties of the soil. The energy multi-scale method takes the microscopic force of the particle as the research object and establishes the CRFW. When the particle size increases, the CRFW of the particle decreases. When the type of the particle changes, the CRFW increases with an increase in particle surface charge.</p>
</list-item>
<list-item>
<p>(3) The energy multi-scale method explains the mechanism of the plasticity index of plastic soil and allowed identification of a new division of soil plasticity. The CRFW quantitatively predicts the influence of mineral composition on the soil plasticity index. When the CRFW is used as the control index, even if mineral composition and particle size of materials are different, the plasticity index of the materials is the same. The multi-scale method may also solve the scale effect of soils in terms of energy scale.</p>
</list-item>
</list>
</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s5">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/supplementary material; further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s6">
<title>Author contributions</title>
<p>Conceptualization, JC and FY; methodology, JC; resources, HT; data curation, XH; writing&#x2014;original draft preparation, JC; writing&#x2014;review and editing, JY; supervision, YF; funding acquisition, JY.</p>
</sec>
<sec id="s7">
<title>Funding</title>
<p>The study was supported by the following fundations: Guangdong Basic and Applied Basic Research Fundation No.20223A1515030051 Funder:JY. Guangdong Basic and Applied Basic Research Fundation No.2022A1515110793 Funder: JC.</p>
</sec>
<sec sec-type="COI-statement" id="s8">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s9">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ashby</surname>
<given-names>M. F.</given-names>
</name>
</person-group> (<year>1970</year>). <article-title>The deformation of plastically non-homogeneous materials</article-title>. <source>Philos. Mag.</source> <volume>21</volume> (<issue>170</issue>), <fpage>399</fpage>&#x2013;<lpage>424</lpage>. <pub-id pub-id-type="doi">10.1080/14786437008238426</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bai</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Cai</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Coupled thermo-hydro-mechanical mechanism in view of the soil particle rearrangement of granular thermodynamics</article-title>. <source>Comput. Geotechnics</source> <volume>137</volume> (<issue>8</issue>), <fpage>104272</fpage>. <pub-id pub-id-type="doi">10.1016/j.compgeo.2021.104272</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Borst</surname>
<given-names>R. D.</given-names>
</name>
<name>
<surname>Alizadeh</surname>
<given-names>S. S.</given-names>
</name>
<name>
<surname>Hageman</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Non-associated Cosserat plasticity</article-title>. <source>Int. J. Mech. Sci.</source> <volume>230</volume>, <fpage>107535</fpage>. <pub-id pub-id-type="doi">10.1016/j.ijmecsci.2022.107535</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chang</surname>
<given-names>C. S.</given-names>
</name>
<name>
<surname>Yin</surname>
<given-names>Z. Y.</given-names>
</name>
<name>
<surname>Hicher</surname>
<given-names>P. Y.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Micromechanical analysis for interparticle and assembly instability of sand</article-title>. <source>J. Eng. Mech.</source> <volume>137</volume> (<issue>3</issue>), <fpage>155</fpage>&#x2013;<lpage>168</lpage>. <pub-id pub-id-type="doi">10.1061/(asce)em.1943-7889.0000204</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Fang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Gu</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Shu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Ba</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Study on pore size effect of low permeability clay seepage</article-title>. <source>ARABIAN J. GEOSCIENCES</source> <volume>12</volume> (<issue>7</issue>), <fpage>238</fpage>. <pub-id pub-id-type="doi">10.1007/s12517-019-4375-3</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Tong</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Tong</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Yuan</surname>
<given-names>J.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>), <article-title>Permeability prediction model modified on kozeny-carman for building foundation of clay soil</article-title>. <source>Buildings</source> <volume>12</volume>(<issue>11</issue>), <fpage>1798</fpage>. <pub-id pub-id-type="doi">10.3390/buildings12111798</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chiu</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Higher-order surface stress effects on buckling of nanowires under uniaxial compression</article-title>. <source>Procedia Eng.</source> <volume>10</volume>, <fpage>397</fpage>&#x2013;<lpage>402</lpage>. <pub-id pub-id-type="doi">10.1016/j.proeng.2011.04.067</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>The dilatancy and numerical simulation of failure behaviorof granular materials based on Cosserat model</article-title>. <source>Interact. Multiscale Mech. Int. J.</source> <volume>5</volume> (<issue>2</issue>), <fpage>157</fpage>&#x2013;<lpage>168</lpage>. <pub-id pub-id-type="doi">10.12989/imm.2012.5.2.157</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Cosserate</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>1909</year>). <source>Th&#xe9;orie de corps deformable</source>. <publisher-loc>Pairs</publisher-loc>: <publisher-name>Herman</publisher-name>.</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dai</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>A numerical analysis of the shear behavior of granular soil with fines</article-title>. <source>Particuolog</source> <volume>2015</volume>(<issue>4</issue>), <fpage>13</fpage>. <pub-id pub-id-type="doi">10.1016/j.partic.2014.08.010</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Xie</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Numerical modeling of open TBM tunneling in stratified rock masses using a coupled FDM-DEM method</article-title>. <source>Comput. Geotechnics</source> <volume>156</volume>, <fpage>105251</fpage>. <pub-id pub-id-type="doi">10.1016/j.compgeo.2023.105251</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fang</surname>
<given-names>Y. G.</given-names>
</name>
<name>
<surname>Bo</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Multiscale problems and analysis of soil mechanics</article-title>. <source>Mech. Mater.</source> <volume>103</volume>, <fpage>55</fpage>&#x2013;<lpage>67</lpage>. <pub-id pub-id-type="doi">10.1016/j.mechmat.2016.09.003</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gong</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Nie</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Exploring the effects of particle shape and content of fines on the shear behavior of sand-fines mixtures via the DEM</article-title>. <source>Comput. Geotechnics</source> <volume>10</volume>, <fpage>161</fpage>&#x2013;<lpage>176</lpage>. <pub-id pub-id-type="doi">10.1016/j.compgeo.2018.10.021</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gudehus</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Nuebel</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>Evolution of shear bands in sand</article-title>. <source>Geotechnique</source> <volume>54</volume> (<issue>3</issue>), <fpage>187</fpage>&#x2013;<lpage>201</lpage>. <pub-id pub-id-type="doi">10.1680/geot.54.3.187.36346</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Guo</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2019</year>). "<article-title>A holistic computational model for prediction of clay suspension structure</article-title>." <source>Int. J. Sediment Res.</source> <volume>34</volume> (<issue>34</issue>), <fpage>345</fpage>. <pub-id pub-id-type="doi">10.1016/j.ijsrc.2018.12.002</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Guzzetti</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Peruccacci</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Rossi</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Stark</surname>
<given-names>C. P.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>The rainfall intensity&#x2013;duration control of shallow landslides and debris flows: An update</article-title>. <source>Landslides</source> <volume>5</volume> (<issue>1</issue>), <fpage>3</fpage>&#x2013;<lpage>17</lpage>. <pub-id pub-id-type="doi">10.1007/s10346-007-0112-1</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Helmy</surname>
<given-names>A. K.</given-names>
</name>
</person-group> (<year>1998</year>). <article-title>The limited swelling of montmorillonite</article-title>. <source>J. Colloid &#x26; Interface Sci.</source> <volume>207</volume> (<issue>1</issue>), <fpage>128</fpage>&#x2013;<lpage>129</lpage>. <pub-id pub-id-type="doi">10.1006/jcis.1998.5739</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Huang</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2018</year>). <source>Division of energy scale and analysis of engineering properties for soil granular group</source>. <publisher-loc>Guangzhou</publisher-loc>: <publisher-name>South China University of Technology</publisher-name>.</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>C. S.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>S. S.</given-names>
</name>
<name>
<surname>Shi</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Computed tomography based numerical simulation for triaxial test of soil&#x2013;rock mixture</article-title>. <source>Comput. Geotechnics</source> <volume>73</volume>, <fpage>179</fpage>&#x2013;<lpage>188</lpage>. <pub-id pub-id-type="doi">10.1016/j.compgeo.2015.12.005</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Mi</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Effects of surface tension and Steigmann&#x2013;Ogden surface elasticity on Hertzian contact properties</article-title>. <source>Int. J. Eng. Sci.</source> <volume>145</volume>, <fpage>103165</fpage>. <pub-id pub-id-type="doi">10.1016/j.ijengsci.2019.103165</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wan</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>A bridging scale method for granular materials with discrete particle assembly &#x2013; Cosserat continuum modeling</article-title>. <source>Comput. Geotechnics</source> <volume>38</volume> (<issue>8</issue>), <fpage>1052</fpage>&#x2013;<lpage>1068</lpage>. <pub-id pub-id-type="doi">10.1016/j.compgeo.2011.07.001</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Deng</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Further development of distinct lattice spring model for stability and collapse analysis of deep foundation pit excavation</article-title>. <source>Comput. geotechnics</source> <volume>144</volume>, <fpage>104619</fpage>. <pub-id pub-id-type="doi">10.1016/j.compgeo.2021.104619</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Ge</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zeng</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Lyu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>Y.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Debris flow susceptibility based on the connectivity of potential material sources in the Dadu River Basin</article-title>. <source>Eng. Geol.</source> <volume>312</volume>, <fpage>106947</fpage>. <pub-id pub-id-type="doi">10.1016/j.enggeo.2022.106947</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhixiang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>A coupled SPH-DEM-FEM model for fluid-particlestructure interaction and a case study of Wenjia gully debris flow impact estimation</article-title>. <source>Landslides</source> <volume>18</volume> (<issue>7</issue>), <fpage>2403</fpage>&#x2013;<lpage>2425</lpage>. <pub-id pub-id-type="doi">10.1007/s10346-021-01640-6</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>A unified size-dependent plate model based on nonlocal strain gradient theory including surface effects</article-title>. <source>Appl. Math. Model.</source> <volume>68</volume>, <fpage>583</fpage>&#x2013;<lpage>602</lpage>. <pub-id pub-id-type="doi">10.1016/j.apm.2018.11.023</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Cui</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>The influence of strain gradient on the shear band in a saturated soil</article-title>. <source>Iran. J. Sci. Technol. Transaction B Eng.</source> <volume>27</volume> (<issue>1</issue>), <fpage>57</fpage>&#x2013;<lpage>62</lpage>.</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mingjing</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>New paradigm for modern soil mechanics: Geomechanics from micro to macro</article-title>. <source>Chin. J. Geotechnical Eng.</source> <volume>41</volume> (<issue>2</issue>), <fpage>195</fpage>&#x2013;<lpage>254</lpage>. <pub-id pub-id-type="doi">10.11779/CJGE201902001</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Missana</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Adell</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>On the applicability of DLVO theory to the prediction of clay colloids stability</article-title>. <source>J. Colloid Interface Sci.</source> <volume>230</volume> (<issue>1</issue>), <fpage>150</fpage>&#x2013;<lpage>156</lpage>. <pub-id pub-id-type="doi">10.1006/jcis.2000.7003</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Mitchell</surname>
<given-names>J. K.</given-names>
</name>
<name>
<surname>Soga</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2005</year>). <source>Fundamentals of soil behavior</source>. <publisher-loc>Hoboken</publisher-loc>: <publisher-name>John Wiley &#x26; Sons</publisher-name>.</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Monkul</surname>
<given-names>M. M.</given-names>
</name>
<name>
<surname>Ozden</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Compressional behavior of clayey sand and transition fines content</article-title>. <source>Eng. Geol.</source> <volume>89</volume> (<issue>3-4</issue>), <fpage>195</fpage>&#x2013;<lpage>205</lpage>. <pub-id pub-id-type="doi">10.1016/j.enggeo.2006.10.001</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>M&#xfc;hlhaus</surname>
<given-names>Dr.</given-names>
</name>
</person-group> (<year>1989</year>). <article-title>Application of Cosserat theory in numerical solutions of limit load problems</article-title>. <source>Ingenieur-Archiv</source> <volume>59</volume> (<issue>2</issue>), <fpage>124</fpage>&#x2013;<lpage>137</lpage>. <pub-id pub-id-type="doi">10.1007/bf00538366</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>M&#xfc;hlhaus</surname>
<given-names>H. B.</given-names>
</name>
<name>
<surname>Oka</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>1996</year>). <article-title>Dispersion and wave propagation in discrete and continuous models for granular materials</article-title>. <source>Int. J. Solids Struct.</source> <volume>33</volume> (<issue>19</issue>), <fpage>2841</fpage>&#x2013;<lpage>2858</lpage>. <pub-id pub-id-type="doi">10.1016/0020-7683(95)00178-6</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nguyen</surname>
<given-names>T. K.</given-names>
</name>
<name>
<surname>Combe</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Caillerie</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Desrues</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>FEM &#xd7; DEM modelling of cohesive granular materials: Numerical homogenisation and multi-scale simulations</article-title>. <source>Acta Geophys.</source> <volume>62</volume> (<issue>5</issue>), <fpage>1109</fpage>&#x2013;<lpage>1126</lpage>. <pub-id pub-id-type="doi">10.2478/s11600-014-0228-3</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ni</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Rheological study on influence of mineral composition on viscoelastic properties of clay</article-title>. <source>Appl. Clay Sci.</source> <volume>187</volume>, <fpage>105493</fpage>. <pub-id pub-id-type="doi">10.1016/j.clay.2020.105493</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Oda</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Iwashita</surname>
<given-names>K.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Study on couple stress and shear band development in granular media based on numerical simulation analyses</article-title>. <source>Int. J. Eng. Sci.</source> <volume>38</volume> (<issue>15</issue>), <fpage>1713</fpage>&#x2013;<lpage>1740</lpage>. <pub-id pub-id-type="doi">10.1016/s0020-7225(99)00132-9</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pham</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Choi</surname>
<given-names>H. J.</given-names>
</name>
<name>
<surname>Lee</surname>
<given-names>I. M.</given-names>
</name>
<name>
<surname>Choi</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>A numerical framework for infinite slope stability analysis under transient unsaturated seepage conditions</article-title>. <source>Eng. Geol.</source> <volume>243</volume>, <fpage>36</fpage>&#x2013;<lpage>49</lpage>. <pub-id pub-id-type="doi">10.1016/j.enggeo.2018.05.021</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Phan</surname>
<given-names>Q. T.</given-names>
</name>
<name>
<surname>Bui</surname>
<given-names>H. H.</given-names>
</name>
<name>
<surname>Nguyen</surname>
<given-names>G. D.</given-names>
</name>
<name>
<surname>Bouazza</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Effect of particle rolling resistance on drained and undrained behaviour of silty sand</article-title>. <source>Acta Geotech.</source> <volume>16</volume> (<issue>3</issue>), <fpage>2657</fpage>&#x2013;<lpage>2682</lpage>. <pub-id pub-id-type="doi">10.1007/s11440-020-01128-y</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Riahi</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Curran</surname>
<given-names>J. H.</given-names>
</name>
<name>
<surname>Bidhendi</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Buckling analysis of 3D layered structures using a Cosserat continuum approach</article-title>. <source>Comput. Geotechnics</source> <volume>36</volume> (<issue>7</issue>), <fpage>1101</fpage>&#x2013;<lpage>1112</lpage>. <pub-id pub-id-type="doi">10.1016/j.compgeo.2009.03.012</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sabbar</surname>
<given-names>A. S.</given-names>
</name>
<name>
<surname>Chegenizadeh</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Nikraz</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Static liquefaction of very loose sand-slag-bentonite mixtures</article-title>. <source>Soils Found. -Tokyo-</source> <volume>57</volume> (<issue>3</issue>), <fpage>341</fpage>&#x2013;<lpage>356</lpage>. <pub-id pub-id-type="doi">10.1016/j.sandf.2017.05.003</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>SampradaPradhan</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Toll</surname>
<given-names>D. G.</given-names>
</name>
<name>
<surname>Rosser</surname>
<given-names>N. J.</given-names>
</name>
<name>
<surname>Brain</surname>
<given-names>M. J.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>An investigation of the combined effect of rainfall and road cut on landsliding</article-title>. <source>Eng. Geol.</source> <volume>307</volume>, <fpage>106787</fpage>. <pub-id pub-id-type="doi">10.1016/j.enggeo.2022.106787</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Cosserat continuum model and its application to the studies of progressive failure</article-title>. <source>Jpn. Geotech. Soc. Spec. Publ.</source> <volume>2</volume> (<issue>18</issue>), <fpage>703</fpage>&#x2013;<lpage>708</lpage>. <pub-id pub-id-type="doi">10.3208/jgssp.chn-57</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wei</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Elastoplastic Cosserat continuum model considering strength anisotropy and its application to the analysis of slope stability</article-title>. <source>Comput. Geotechnics</source> <volume>117</volume>, <fpage>103235</fpage>. <pub-id pub-id-type="doi">10.1016/j.compgeo.2019.103235</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wei</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>An anisotropic elastoplastic Cosserat continuum model for shear failure in stratified geomaterials</article-title>. <source>Eng. Geol.</source> <volume>293</volume> (<issue>1</issue>), <fpage>106304</fpage>. <pub-id pub-id-type="doi">10.1016/j.enggeo.2021.106304</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Thornton</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Numerical simulations of deviatoric shear deformation of granular media</article-title>. <source>Geotechnique</source> <volume>50</volume> (<issue>1</issue>), <fpage>43</fpage>&#x2013;<lpage>53</lpage>. <pub-id pub-id-type="doi">10.1680/geot.2000.50.1.43</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Unified criterion for stability failure of rock specimens subject to uniaxial tension, direct shear and uniaxial compression based on gradient-dependent plasticity</article-title>. <source>Rock Soil Mech.</source> <volume>24S</volume>, <fpage>138</fpage>&#x2013;<lpage>142</lpage>. <pub-id pub-id-type="doi">10.16285/j.rsm.2003.s2.032</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Shen</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Two-dimensional fretting contact analysis considering surface effects</article-title>. <source>Int. J. Solids Struct.</source> <volume>170</volume>, <fpage>68</fpage>&#x2013;<lpage>81</lpage>. <pub-id pub-id-type="doi">10.1016/j.ijsolstr.2019.04.027</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>H. Y.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Scale effect in macroscopic permeability of jointed rock mass using a coupled stress&#x2013;damage&#x2013;flow method</article-title>. <source>Eng. Geol.</source> <volume>228</volume>, <fpage>121</fpage>&#x2013;<lpage>136</lpage>. <pub-id pub-id-type="doi">10.1016/j.enggeo.2017.07.009</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Kang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>A unified high-order model for size-dependent vibration of nanobeam based on nonlocal strain/stress gradient elasticity with surface effect</article-title>. <source>Int. J. Eng. Sci.</source> <volume>182</volume>, <fpage>103785</fpage>. <pub-id pub-id-type="doi">10.1016/j.ijengsci.2022.103785</pub-id>
</citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zeghal</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Shamy</surname>
<given-names>U. E.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Liquefaction of saturated loose and cemented granular soils</article-title>, <source>Powder Technol.</source>, <volume>184</volume> <fpage>254</fpage>&#x2013;<lpage>265</lpage>.<pub-id pub-id-type="doi">10.1016/j.powtec.2007.11.032</pub-id>
</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>DEM-FEM based numerical analysis on mechanical responses of sandy soil and pipeline to seepage erosion</article-title>. <source>Eng. Geol.</source> <volume>310</volume> (<issue>5</issue>), <fpage>106868</fpage>. <pub-id pub-id-type="doi">10.1016/j.enggeo.2022.106868</pub-id>
</citation>
</ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Feng</surname>
<given-names>D.</given-names>
</name>
</person-group>, (<year>2017</year>). <article-title>Coupled DEM-CFD investigation on the formation of landslide dams in narrow rivers</article-title>. <source>Landslides</source> <volume>14</volume> (<issue>1</issue>), <fpage>189</fpage>&#x2013;<lpage>201</lpage>. <pub-id pub-id-type="doi">10.1007/s10346-015-0675-1</pub-id>
</citation>
</ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Dupla</surname>
<given-names>J. C.</given-names>
</name>
<name>
<surname>Canou</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Foerster</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Investigation on the undrained shear strength of loose sand with added materials at various mean diameter ratios</article-title>. <source>Soil Dyn. Earthq. Eng.</source> <volume>137</volume>, <fpage>106276</fpage>. <pub-id pub-id-type="doi">10.1016/j.soildyn.2020.106276</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>