Abstract
Based on the granular-solid-hydrodynamic theory, the constitutive model considering the thermo-hydro-mechanical (THM) coupled action is established, and the dilatancy property of sandy soil under coupled high mechanical pressure and temperature is simulated. The relationship between the energy dissipation and the macroscopic stress-strain changes at the grain level of saturated sandy soil is connected by defining the transfer coefficient and the energy function, without considering the concepts of yield surface and hardening parameters in classical plastic mechanics. Additionally, the changes in temperature, relative density and confining pressure during the shearing process cause particle rolling, slipping and friction. The energy dissipation in this process is described by defining the concept of particle entropy and particle temperature. In the calculation, the isotropic compression test, drained and undrained shear test of sandy soil under high stress are simulated respectively. The validity of the model is proved by comparing with the test results. Meanwhile, the stress-strain relationship and pore pressure variation law of sandy soil under different temperatures are predicted. The results show that the effect of temperature on shear strength is limited, and the pore pressure will gradually increase and become stable with the increase of temperature. Thus, this work establishes the soil THM coupled model from the perspective of micro energy dissipation, which can provide new theoretical support for the prediction of natural disasters such as landslides and debris flow.
Introduction
In mountainous areas, under the action of high ground stress, water flow and temperature, sandy soil will be disturbed by the construction of gravity dams, deep-buried tunnels and other large-scale projects, which will generate pore water pressure in the sandy soil, and cause particles to roll or even break (; ). This process intensifies the occurrence of natural disasters such as debris flow and landslide, and seriously affects the environmental change and sustainable development in mountainous areas (; Zhou et al., 2021). Through the triaxial drained and undrained shear tests of sandy soil under high confining pressure and temperature, researchers found that the higher the confining pressure, the more obvious the strain softening trend of sandy soil samples. Meanwhile, with the increase of the confining pressure, the axial strain at failure reaches the peak value, then decreases, and the volumetric strain gradually changes from dilatancy to shear shrinkage (; ). Attempts have been made to consider the factors that influence the mechanical properties of sandy soil including initial density, temperature and stress level. Researchers (; ; ) introduced that the sand particles were broken obviously under high stress, and the isotropic compression curves of different initial density samples finally coincided with each other. At the same time, when the confining pressure is large enough, the initial relative density has little effect on the stress-strain relationship and pore pressure, and the sandy soil sample shows the nature of normally consolidated clay. divided the compression characteristics of sandy soil into three stages through one-dimensional compression test: under low pressure, the volume of sample decreases due to the rearrangement of particles; under high pressure, the crushing and rearrangement of sand particles produce stronger compression performance; under higher pressure, the contact between particles increases greatly, and the material properties are controlled by particle crushing. The results of the above experiments provide guidance for the establishment of constitutive relationship of sandy soil under high stress.
In the earlier studies considering the impact of particle breakage on the mechanical properties of sandy soil, a variety of mathematical models were established by the modification and expansion of the existing soil model. Scholars (; ; ) have improved the yield surface or hardening law in the Cambridge model and established the constitutive model of sandy soil by using the unrelated flow law, which can well predict the characteristics of sandy soil. Meanwhile, the theory of sub plasticity based on thermodynamics has been increasingly applied in the establishment of sand’s constitutive model. The theory uses tensor function as a tool, abandons the concepts of strain decomposition into elastic strain and plastic strain, yield surface and loading and unloading in traditional elastic-plastic theory, and directly establishes the relationship between stress rate and strain rate (; ; ). In order to describe the material state change of sandy soil during loading, some scholars introduced state parameters into the constitutive relation of sandy soil. successfully carried out a series of triaxial test against Kogyuk sand and proposed the concept of state parameters, which have been used by many researchers to study the dilatancy of sandy soil. For better reflection of the deformation behavior of saturated sandy soil in general stress space, established a constitutive model, which was based on the basic framework of critical state soil mechanics, and combined the theory of bilateral plastic model with state parameters. verified that the dilatancy of granular soil was controlled by the distribution of contact stress between particles through experiments on different materials, and deduced the stress dilatancy relationship under plane strain condition on the assumption that the minimum energy ratio is reached when the soil is damaged. proposed to use dilatancy index to reflect the current state of sandy soil, which can estimate the dilatancy of sandy soil samples. However, the constitutive equations used in these models cannot effectively describe the coupled characteristic of soil under multi field action, and are only suitable for soil samples under a certain condition.
Different from the above methods, some scholars used the granular-solid-hydrodynamic (GSH) theory, which has been applied in solid materials such as polymers and crystals, to establish the constitutive model of soil. The GSH theory is based on the method of fluid dynamics, and constructs the constitutive relationship of materials on the basis of non-equilibrium theory and physical conservation law. Using the definitions of particle entropy and particle temperature, established the theoretical constitutive model of single granular material, which can comprehensively consider the macro and micro behaviour of granular material. On this basis, some researchers (Zhang and Cheng, 2017; ) extended the GSH theory to saturated soil and unsaturated soil, respectively, to simulate the mechanical characteristics of related soil, which confirmed the effectiveness of the model. further extended the theory and focused on the factors that influence the characteristics of unsaturated soil including over consolidation ratio, temperature and stress path. The traditional plastic theory is to choose the yield function, plastic potential surface and hardening function independently, which is not rigorous enough in theory, and GSH theory can effectively avoid this point (). In addition, the constitutive equation of soil based on the GSH theory can automatically meet the first and second laws of thermodynamics. However, taking into account the complexity of the interaction between the various phases in the soil, the application of the theory in various types of soil needs to be further promoted.
In the present study, based on the GSH theory and the dilatancy equation considering the state variables, a THM coupled model is established, and the shear test of saturated sandy soil under the coupled action of temperature and stress is simulated. The new THM coupled model considers the energy dissipation of sandy soil material from the micro mechanism by using the linear non-equilibrium theory and the critical state theory, and further defines the entropy and temperature of the particle to link the energy dissipation at particle level with the change of system energy. The effectiveness of the proposed model is verified by the simulation of sandy soil samples under different initial void ratios and confining pressures. The simulation tests include isotropic compression test, drained and undrained shear test. At the same time, based on undrained test, the stress-strain relationship and pore pressure change of sandy soil under different temperatures are predicted, which provides a reference for the development of follow-up experiments.
Methods
Granular-Solid-Hydrodynamic Theory in Sandy Soil
It is assumed that both the solid and liquid phases of geomaterials are incompressible and have the same temperature at the same location. In the closed system, the change of internal energy of granular materials caused by unit time and unit volume can be divided into two categories: recoverable energy and irrecoverable energy. For granular materials, the irrecoverable energy dissipation is mainly due to the sliding, collision, rolling and friction between granular materials, while the recoverable energy is mainly due to the mutual extrusion between particles, which transforms the external input energy into elastic potential energy for storage. Thus, the composition of energy density of granular materials in the closed system can be obtained, as shown in Eq. 1:where is the elastic potential energy stored by the granular material under the action of external force; is the irrecoverable energy produced in the mesoscopic scale; is the kinetic energy of the granular material itselfp; w is the energy density in isolated system; is the heat energy of the material itself in the closed system.
In addition, when characterizing the independent state variables of granular materials, there are inelastic collisions between particles in the mesoscopic scale of granular materials. This further results in the dissipation of system energy, which shows as the change of system entropy. linked the entropy generation caused by the dissipation of granular materials at mesoscopic level with the change of macroscopic entropy, and determined the thermodynamic independent state variable and their conjugate quantities of granular materials. By constructing thermodynamic identities, the energy density expression is obtained:
In combination with Eq. 1, it is concluded that:where is the particle temperature; is the elastic stress; is the elastic strain; is the material temperature; is material entropy; is the chemical potential; is the material density; is particle entropy; is the velocity; is momentum.
According to the second law of thermodynamics, the entropy transformation rate R caused by the irreversible process in the closed system is always greater than zero. The higher the entropy of the system, the closer the isolated system will be to the equilibrium state, and the entropy will reach the maximum. The linear non-equilibrium theory defines the factor that drives the closed system to deviate from the equilibrium state as dissipative force, and the energy dissipation caused by the unit dissipative force in unit time is called dissipative flow. At the same time, when the system is close to the equilibrium state, it can be considered that there is a linear relationship between the dissipative force and the dissipative flow. The dissipative force and the corresponding dissipative flow of granular materials in the closed system can be determined by using the conservation equations of mass, momentum and energy combined with the entropy increasing equation of the closed system. It shows that the elastic stress and the plastic strain rate are the corresponding dissipative force and flow. The dissipative force and dissipative flow at the particle level are and , respectively, which can be used to represent the entropy generation rate of particles Rg at the mesoscopic level (). The following results can be obtained:where is the transfer coefficient corresponding to the elastic stress , is the transfer coefficient at the particle level.
According to the basic function form of the elastic potential energy density function of cohesive soil materials given by and the basic concept of the specific elastic energy proposed by , for the saturated sandy soil without considering the influence of the cohesive force, the expression of the elastic potential energy stored in the sandy soil under the action of temperature is obtained as follows:where B is the parameter describing the degree of hardness or softness of the material. is the thermal expansion coefficient of soil, is the elastic volume strain , is the elastic shear strain , is the friction angle coefficient, e is the void ratio, reflects the material property of sandy soil, andis the slope of e-logp curve in consolidation test of sandy soil sample. Combined with the isotropic compression consolidation curve, B is assumed to be the form.
By using the Expressions (3), (4) and (5) and referring to the expressions in reference (; ), the elastic stress and plastic strain rate are obtained:where elastic stress is equal to the effective stress between particles (), is the elastic deviator strain , and are the corresponding transfer coefficients.
When the sandy soil is in the critical state, there exists liquefaction phenomenon. At this time, the elastic volume strain is 0, the average effective stress is 0, and will be taken as 1.5. Combined with Eq. 8, the specific expressions of average effective stress and shear stress q are obtained as follows:
In reference , the energy carried by particle motion is expressed as = . According to the linear non-equilibrium theory, the entropy generation rate at the particle level is . The motion equation of particle entropy at the mesoscopic level is obtained by combining Eqs 4–6
For the convenience of parameter calculation, Eq. 12 is simplified as follows:where is the deviator strain; is the volume strain; and are the transport coefficients at the particle level, respectively; is the dissipation at particle level caused by temperature. b is a material properties parameter, ,,,,. According to the hypothesis in references (; ), a is a material constant and is taken as 0.5 ().
Dilatancy in Sandy Soil
When the sandy soil sample reaches the critical state, the volume strain increment is zero and the deviator strain increment increases continuously in the last stage of shear deformation. The critical state of sandy soil is not a straight line in the e-logp plane (; ). Li and Dafalias () expressed the critical state of sandy soil as follows:where is the void ratio corresponding to the critical state line when p = 0; ec is the critical void ratio; is the slope in-plane; Pa is the standard atmospheric pressure, usually 100kPa; θ is the critical state line parameter.
By comparing the test data of Toyoura sandy soil by with the critical state curve proposed by , it is found that Eq. 14 can better reflect the critical state condition of sandy soil.
The state parameter is usually expressed by relative density and confining pressure. regard the dilatancy of sandy soil as an independent state variable, and proposed the difference between critical void ratio and current void ratio e to represent the state of sandy soil samples:
In this way, the current state (e, p) of sandy soil can be determined. When >0, the sandy soil is in a loose state, and shear shrinkage occurs during shearing. On the other side (<0), the sandy soil is in a relatively dense state, and shear dilatation occurs during shearing.
Referring to , the dilatancy ratio is the ratio of plastic volumetric strain to plastic shear strain , and the state parameter is introduced into the dilatancy equation to establish the dilatancy expression of sandy soil.whereis the critical state stress ratio; = is the stress ratio in the current state; is the state parameter, which is determined by combining Eq. 15; and are the fitting parameters.
The above is the constitutive model of sandy soil under high stress and temperature coupled based on GSH theory. The main governing equations are Eqs 9–11, Eqs 13–16. The model establishes the constitutive relationship from the point of microscopic energy dissipation and can well reflect the mechanical properties of sandy soil. The state variables required for the constitutive equation of GSH theoretical model are shown in Table 1.
TABLE 1
| State variables | Solution method |
|---|---|
| Dry density of sand sample (kg/m3). The initial value is determined by the following formula: where is the density of liquid; initial void ratio; particle specific gravity | |
| Particle entropy temperature.. Dimensionless. Its initial value is 0 | |
| Elastic volume strain, dimensionless, initial value 0 | |
| Elastic shear strain, dimensionless, initial value 0 | |
| T | The test temperature (°C). The initial value T0 is measured according to the initial conditions of the test, generally taking the ambient temperature of 25 (°C) |
| Effective stress of soil (kPa) | |
| Shear stress of soil (kPa) |
State variables in GSH theoretical model.
Determination of Model Parameters
There are 15 parameters in the model, including thermodynamic parameter
and parameter
in the dilatancy equation with state variable. Most of the parameters can be obtained through the conventional soil mechanics test, and the available specific parameter calibration method is provided here.
1) B0 andare calibrated by the curve of isotropic compression and rebound. B0 is the elastic stiffness parameter, and a set of elastic strain under confining pressure are obtained by using the rebound curve of soil. Because the shear strain is ignored in the rebound process, the value of B0 can be obtained by combining Eq. 10. is the slope of e-logp curve in one-dimensional consolidation test of sandy soil sample.
2) ξ is the parameter related to the friction angle of sandy soil (). In the critical state, the volume and stress of sandy soil remain unchanged, while the shear strain develops continuously. In the simulation, it is shown as, . Eq. 9 is simplified as:
where the parameter
is defined.
3) The dilatancy ratio d is the ratio of plastic volumetric strain rate to plastic shear strain rate. Eq. 18 is obtained by combining Eq. 17. Simultaneous Eqs 16–18 determine the value of c1.
4) The value of c2 can be obtained through Eq. 17 and Eq. 13, Eqs 10–11 as well as the critical state ratio.
5) The parameter c3 is determined by isotropic compression curve; c4 is determined by the thermal response caused by heating, and c4 = 0.003 kg m−3 kPa−1 by referring to ; ); c5 is determined by stress relaxation test, and the value in reference is 800 kg m−3 s−1.
6) The method of parameter selection in the dilatancy equation is explained in detail in reference .
Model Verification
Lade and Yamamuro (; ; ; ) carried out a series of experimental studies on the mechanical properties of three different initial densities of Cambria sandy soil under different pressures using a standard triaxial apparatus. The particle size distribution of Cambria sandy soil ranges from 0.83 to 2 mm, the particle specific gravity Gs is 2.69, and the maximum and minimum void ratio are 0.792 and 0.503, respectively. The temperature of sample test is set at 25°C, and the experimental results are estimated by the sandy soil model. Through the parameter analysis of part 2.3, the values of relevant control symbols including thermodynamic model parameters and state parameters are listed in Table 2.
TABLE 2
| Thermodynamic model parameters | Values | State parameters | Values |
|---|---|---|---|
| B0 | 2.65 × 106 MPa | M | 1.370 |
| Gs | 2.69 | ek | 0.074 |
| ξ | 13 | λc | 0.019 |
| λ | 0.367 | θ | 0.537 |
| βT | 4.5 × 10–5°C−1 | d0 | 15.322 |
| c3 | 9,000 s2 | m | 0.328 |
| c4 | 0.003 kg m−3 kPa−1 | ||
| c5 | 800 kg m−3 s−1 |
Thermodynamic parameters and state parameters in the model.
B0 is the elastic stiffness parameter; Gs is the specific gravity of soil; ξ is the parameter related to the friction angle of sandy soil; λ is the slope of e-logp curve in one-dimensional consolidation test of sandy soil sample; βT is the thermal expansion coefficient of soil; c3, c4 and c5 are the simplified parameters related to the transfer coefficients; M is the slope of the critical stress ratio; ek is the void ratio corresponding to the critical state line when p = 0; λc is the slope lne-(p’/pa)θ plane; θ is the critical state line parameter; d0 and m is the fitting parameters.
Isotropic Compression Test
According to Eqs 10, 11, the corresponding differential expression is obtained as follows:
In the isotropic compression experiment, the elastic shear strain is 0, and the effect of temperature is not considered. According to Eq. 19, the stress expression is described as:where is the density of soil particles, obtained by ; is the water density, = 1,000 kg/m3.
Figure 1 illustrates the simulation of sandy soil isotropic compression test using the newly proposed model, where the initial relative densities are 30, 60, and 90%, respectively. The results show that when the confining pressure exceeded 12 MPa, the sandy soil samples with different initial void ratios have similar void ratios under the same confining pressure, indicating that the initial relative density has no significant effect on the deformation and strength of sandy soil when the stress level exceeds a certain value. This is due to the fact that under the action of larger confining pressure, the rolling and extrusion deformation of soil particles in the sandy soil sample gradually turns into the crushing of soil particles. Meanwhile, with the increase of stress level, the amount of broken sandy soil particles increases gradually, and the broken particles fill the gap, resulting in the compression phenomenon of clay like materials. Through comparison, the proposed model can well reflect the compression characteristics of sandy soil under high stress.
FIGURE 1
Uudrained Shear Test
In the undrained shear test of sandy soil, ,, it can be deduced that:
By introducing Eqs 19–20 into Eqs 22, 23, the pore pressure growth rate can be written as:
It can be found from Figure 2 that, at lower confining pressure (6.4 MPa), the deviator stress increases and tends to be stable with the increase of axial strain, and there exists no peak point, showing strain hardening characteristics. However, when the confining pressure is more than 34 MPa, under the condition of small strain, the deviator stress increases rapidly, then the growth rate slows down and the peak value appears. In the final stage, as the increase of axial strain, the deviator stress decreases and tends to be stable. It is worth noting that the higher the confining pressure, the higher the peak strength, and the greater the difference between the peak strength and the bias stress under the stable state. The confining pressure increased from 34 to 52 MPa, the difference between them was 3.67 and 6.65 MPa, and the relative increase rate was 37.7 and 48.6% respectively. This phenomenon can be explained as the breakage of particles in the sample. When the confining pressure is low, at the initial stage of shear deformation, the shear strength of sandy soil samples is mainly affected by sliding friction, dilatancy effect and particle rearrangement, and the particle breakage rate is low. The higher the confining pressure is, the higher the particle breakage rate is under the same axial strain condition. When the shear strength reaches the maximum value, the particle breakage inside the sample becomes the main factor affecting the shear strength, and its shear strength gradually decreases and tends to be stable, which is similar to the softening characteristics of clay soil.
FIGURE 2
Moreover, in conjunction with Figure 3, it can be seen that under the confining pressure of 6.4 MPa, the average effective stress decreases rapidly and then increases. This is mainly due to the fact that in the initial stage of shear deformation, the volume of the sandy soil sample has a faster trend of shrinking, resulting in a rapid increase in pore pressure. With the further development of shear deformation, pore pressure gradually decreases, and the sample shows the tendency of volume expansion and finally the sample is damaged. The simulation results are in good agreement with the experimental results, which verifies the validity of the proposed model.
FIGURE 3
Figure 4 indicated that the stress-strain curves of sandy soil with different initial densities have similar changing trends. Under the condition of small axial strain (<3%), loose sandy soil and medium dense sandy soil reach the same peak strength. However, after the peak value, the shear strength of medium dense sandy soil is always higher than that of loose sandy soil under the same axial strain. This is mainly due to the influence of two aspects on the mechanical properties of sandy soil under larger confining pressure, namely, the rearrangement of particles and particle breakage. The former is affected by initial relative density, while the latter is affected by confining pressure.
FIGURE 4
Three temperature gradients (25°C, 40°C and 80°C) are set according to the method in the reference (; ; ). The undrained shear tests of sandy soil samples at different temperatures were simulated. The stress-strain relationship and pore pressure relation with strain under different temperature conditions are presented in Figure 5 and Figure 6. Figure 5 illustrates that the stress variation law of sandy soil samples under different temperatures is consistent, showing that the higher temperature lead to an increase of the peak strength. It is indicating that within a certain strain range, the increase in temperature can affect the undrained shear strength of sandy soil samples, which is consistent with the conclusion in literature (; ). It is worth noting that when the strain exceeds a certain value (15%), with the further development of shear deformation, the deviatoric stress values under different temperatures are finally close to each other. At this time, the ambient temperature during shear test has no significant effect on the strength of sandy soil. The variation of pore pressure under different temperatures is presented in Figure 6. The pore pressure increased gradually with increasing shear deformation and finally reached a constant value until axial strain reached 10%, which is caused by the decrease trend of volume in the initial stage of shear. At the same time, the comparison of pore pressure changes at different temperatures indicate that the higher the temperature is, the higher the pore pressure is, which is mainly caused by the difference of thermal expansion coefficient between water and sandy soil particles.
FIGURE 5
FIGURE 6
Drained Shear Test
In the sandy soil drainage shear test, the relationship between shear stress and average effective stress are as follows:
By substituting Eqs 19, 20 into Eq. 25, the specific expression of volume strain rate without considering temperature effect is obtained as follows:
Specifying that the volume strain of the specimen is negative with volume shrinkage and positive with volume expansion. It can be seen in Figure 7 that when the confining pressure is 2.1 MPa, with the axial strain increases, the sandy soil sample experiences two trends: shear shrinkage first and then shear expansion, until the sample is destroyed. This is in line with the volume change trend of medium density sandy soil under the action of low confining pressure. Figure 8 demonstrates that, at the lower confining pressure (2.1 MPa), the relationship between axial strain and shear stress ratio can be divided into two stages: surge section and platform section. In the initial stage, when the axial strain is in 0–5%, the stress ratio increases obviously. With the increase of axial strain, an obvious platform section appears, and the stress increase trend is gradually gentle until the sandy soil sample is destroyed. In contrast, when the confining pressure gradually increases (5.8 and 11.5 MPa), the volumetric strain presents the deformation dominated by shear shrinkage. Meanwhile, with the increase of confining pressure, the initial slope of stress-strain curve decreases, the platform section gradually decreases, and the maximum principal stress ratio of the specimen decreases. The main reason for the formation of the platform section is that under the action of small confining pressure (2.1 MPa), the load-bearing particles in the sandy soil sample change from the initial mutual extrusion and rolling to the particle breakage and particle reorientation, which is also the main reason for the specimen deformation from shear shrinkage to shear dilatation. This conclusion is consistent with the point of reference (Zhou et al., 2021). When the confining pressure is large (5.8 and 11.5 MPa), the plateau section in the stress-strain curve is no longer obvious, indicating that a large number of particles are broken in the initial stage and distributed throughout the test process, which can be verified by the fact that the slope of volume strain in Figure 7 remains basically constant.
FIGURE 7
FIGURE 8
Conclusion
By introducing the state variable into the dilatancy equation and using the theory of GSH and critical state, an inelastic coupled model of thermo-hydro-mechanical is established without considering the concepts of yield surface and flow rule, and the characteristics of sandy soil under the action of high stress and temperature are simulated. The effectiveness of the model is proved by the isotropic compression tests and the shear tests of drained and undrained. Moreover, based on the undrained shear test of sandy soil, the stress-strain relationship and pore pressure of sandy soil under different temperature conditions are predicted. The following conclusions are drawn:
Under the action of high stress, the microscopic mechanism of sandy soil sample deformation is gradually transformed from the rolling and extrusion of the particles to the particles broken and filled with pores. The particle breakage rate and the degree of shear shrinkage deformation of sandy soil samples are positively correlated with the pressure on the samples.
The stress-strain relationship of sandy soil under atmospheric pressure is strain hardening type, which gradually changes to strain softening type with the increase of confining pressure, and the peak strength value increases with the increase of confining pressure.
The higher the temperature is, the greater the pore pressure of sandy soil under undrained shear condition are. However, with the increase of shear strain, the effect of temperature on shear strength is gradually weakened, while the pore pressure is gradually increased and tends to be flat.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
The contribution of RZ is methodology, and the others are investigation and analysis.
Funding
Hebei Province Postdoctoral Research Projects Merit-based Funding Program (B2020005008).
Conflict of interest
WW was employed by Hebei Research Institute of Construction and Geotechnical Investigation Co. Ltd.
The remaining 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.
Publisher’s note
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.
References
1
Abuel-NagaH. M.BergadoD. T.LimB. F. (2007). Effect of Temperature on Shear Strength and Yielding Behavior of Soft Bangkok clay. Soils Found.47 (3), 423–436. 10.3208/sandf.47.423
2
BaiB.LiT. (2013). Irreversible Consolidation Problem of a Saturated Porothermoelastic Spherical Body with a Spherical Cavity. Appl. Math. Model.37 (4), 1973–1982. 10.1016/j.apm.2012.05.003
3
BaiB.ShiX. (2017). Experimental Study on the Consolidation of Saturated Silty clay Subjected to Cyclic thermal Loading. Geomech. Eng.12 (4), 707–721. 10.12989/gae.2017.12.4.707
4
BaiB.SuZ. (2012). Thermal Responses of Saturated Silty Clay during Repeated Heating-Cooling Processes. Transp. Porous Med.93 (1), 1–11. 10.1007/s11242-012-9939-6
5
BaiB.GuoL.HanS. (2014). Pore Pressure and Consolidation of Saturated Silty clay Induced by Progressively Heating/cooling. Mech. Mater.75, 84–94. 10.1016/j.mechmat.2014.04.005
6
BaiB.LongF.RaoD.XuT. (2017). The Effect of Temperature on the Seepage Transport of Suspended Particles in a Porous Medium. Hydrol. Process.31 (2), 382–393. 10.1002/hyp.11034
7
BaiB.RaoD.XuT.ChenP. (2018). SPH-FDM Boundary for the Analysis of thermal Process in Homogeneous media with a Discontinuous Interface. Int. J. Heat Mass Transfer117, 517–526. 10.1016/j.ijheatmasstransfer.2017.10.004
8
BaiB.RaoD.ChangT.GuoZ. (2019). A Nonlinear Attachment-Detachment Model with Adsorption Hysteresis for Suspension-Colloidal Transport in Porous media. J. Hydrol.578, 124080. 10.1016/j.jhydrol.2019.124080
9
BaiB.YangG.-c.LiT.YangG.-s. (2019). A Thermodynamic Constitutive Model with Temperature Effect Based on Particle Rearrangement for Geomaterials. Mech. Mater.139, 103180. 10.1016/j.mechmat.2019.103180
10
BaiB.ZhangJ.LiuL.JiY. (2020). The Deposition Characteristics of Coupled lead Ions and Suspended Silicon Powders along the Migration Distance in Water Seepage. Transp. Porous Med.134 (3), 707–724. 10.1007/s11242-020-01464-3
11
BaiB.XuT.NieQ.LiP. (2020). Temperature-driven Migration of Heavy Metal Pb2+ along with Moisture Movement in Unsaturated Soils. Int. J. Heat Mass Transfer153, 119573. 10.1016/j.ijheatmasstransfer.2020.119573
12
BaiB.NieQ.ZhangY.WangX.HuW. (2021). Cotransport of Heavy Metals and SiO2 Particles at Different Temperatures by Seepage. J. Hydrol.597, 125771. 10.1016/j.jhydrol.2020.125771
13
BaiB.JiangS.LiuL.LiX.WuH. (2021). The Transport of Silica Powders and lead Ions under Unsteady Flow and Variable Injection Concentrations. Powder Techn.387, 22–30. 10.1016/j.powtec.2021.04.014
14
BaiB.NieQ.WuH.HouJ. (2021). The Attachment-Detachment Mechanism of Ionic/nanoscale/microscale Substances on Quartz Sand in Water. Powder Techn.394, 1158–1168. 10.1016/j.powtec.2021.09.051
15
BaiB. (2006). Fluctuation Responses of Saturated Porous media Subjected to Cyclic thermal Loading. Comput. Geotechn.33 (8), 396–403. 10.1016/j.compgeo.2006.08.005
16
BaiB. (2006). Thermal Consolidation of Layered Porous Half-Space to Variable thermal Loading. Appl. Math. Mech.27, 1531–1539. 10.1007/s10483-006-1111-1
17
BeenK.JefferiesM. G. (1985). A State Parameter for Sands. Géotechnique35 (2), 99–112. 10.1680/geot.1985.35.2.99
18
BillamJ. (1971). “Some Aspects of the Behaviour of Granular Materials at High Pressures,” in Proceedings of the Roscoe Memorial Symposium, Cambridge, March (Cambridge University), 69–80.
19
BoltonM. D. (1986). The Strength and Dilatancy of Sands. Géotechnique36 (1), 65–78. 10.1680/geot.1986.36.1.65
20
CuiC.MengK.XuC.LiangZ.LiH.PeiH. (2021). Analytical Solution for Longitudinal Vibration of a Floating Pile in Saturated Porous media Based on a Fictitious Saturated Soil Pile Model. Comput. Geotechn.131, 103942. 10.1016/j.compgeo.2020.103942
21
HagertyM. M.HiteD. R.UllrichC. R.HagertyD. J. (1993). One‐Dimensional High‐Pressure Compression of Granular Media. J. Geotechn. Eng.119 (1), 1–18. 10.1061/(asce)0733-9410(1993)119:1(1)
22
JiangY.LiuM. (2007). A Brief Review of “Granular Elasticity”. Eur. Phys. J. E22, 255–260. 10.1140/epje/e2007-00009-x
23
JiangY.LiuM. (2009). Granular Solid Hydrodynamics. Granular Matter11 (3), 139–156. 10.1007/s10035-009-0137-3
24
JiangY.EinavI.LiuM. (2017). A Thermodynamic Treatment of Partially Saturated Soils Revealing the Structure of Effective Stress. J. Mech. Phys. Sol.100, 131–146. 10.1016/j.jmps.2016.11.018
25
LadeP.YamamuroJ. (1997). Undrained Sand Behavior in Axisymmetric Tests at High Pressures. J. Geotechn. Eng.122 (2), 120–129.
26
LadeP.YamamuroJ.BoppP. (2005). “Relative Density Effects on Drained and Undrained Strengths of Sand at High Pressures,” in Proceedings of the 16th International Conference on Soil Mechanics and Geotechnical Engineering: Geotechnology in Harmony with the Global Environment, Osaka, Japan, 537–541. 10.3233/978-1-61499-656-9-537
27
LiX. S.DafaliasY. F. (2000). Dilatancy for Cohesionless Soils. Géotechnique50 (4), 449–460. 10.1680/geot.2000.50.4.449
28
ManzatiM.DafaliasY. (1997). A Critical State Two-Surface Plasticity Model for Sands. Geotechnique47 (2), 255–272.
29
MengK.CuiC.LiangZ.LiH.PeiH. (2020). A New Approach for Longitudinal Vibration of a Large-Diameter Floating Pipe Pile in Visco-Elastic Soil Considering the Three-Dimensional Wave Effects. Comput. Geotechn.128, 103840. 10.1016/j.compgeo.2020.103840
30
PoulV.LadeP.BoppP. (2020). Relative Density Effects on Drained Sand Behavior at High Pressures. Soils and Foundations45 (1), 1–13. 10.3208/sandf.45.1_1
31
RoweP. (1962). The Stress–Dilatancy Relation for Static Equilibrium of an Assembly of Particles in Contact. Proc. R. Soc. Lond.269 (1339), 500–527.
32
TaherdangkooR.TatomirA.SauterM. (2020a). Modeling of Methane Migration from Gas Wellbores into Shallow Groundwater at basin Scale. Environ. Earth Sci.79 (432), 3–16. 10.1007/s12665-020-09170-5
33
TaherdangkooR.TatomirA.TaherdangkooM.QiuP.SauterM. (2020b). Nonlinear Autoregressive Neural Networks to Predict Hydraulic Fracturing Fluid Leakage into Shallow Groundwater. Water12 (3), 841. 10.3390/w12030841
34
TaherdangkooR.LiuQ.XingY.YangH.CaoV.SauterM.et al (2021). Predicting Methane Solubility in Water and Seawater by Machine Learning Algorithms: Application to Methane Transport Modeling. J. Contam. Hydrol.242, 103844. 10.1016/j.jconhyd.2021.103844
35
VerdugoR.IshiharaK. (1996). The Steady State of sandy Soils. Soils Found.36 (2), 81–91. 10.3208/sandf.36.2_81
36
WangY.ChenF.LiX.YinX.LeiY. (2020). The Variable-Mass Seepage Law of Broken Porous Rock: an Experimental Study. Geomatics Nat. Hazards Risk11 (1), 1991–2005. 10.1080/19475705.2020.1821791
37
WangY.LiZ.JingH.LiY.WangM. (2021). Study on the Seepage Characteristics of Deep Buried Tunnels under Variable High-Pressure Water Heads. Bull. Eng. Geol. Environ.80 (2), 1477–1487. 10.1007/s10064-020-01986-6
38
WuW.KolymbasD. (1990). Numerical Testing of the Stability Criterion for Hypoplastic Constitutive Equations. Mech. Mater.9, 245–253. 10.1016/0167-6636(90)90006-2
39
YamamuroJ. (1996). Drained Sand Behavior in Axisymmetric Tests at High Pressures. J. Geotechn. Eng.122, 109. 10.1061/(asce)0733-9410(1996)122:2(109)
40
YangG.BaiB. (2019). Thermo-hydro-mechanical Model for Unsaturated clay Soils Based on Granular Solid Hydrodynamics Theory. Int. J. Geomech.19 (10), 04019115. 10.1061/(asce)gm.1943-5622.0001498
41
ZarifiM.KvammeB.KuznetsovaT. (2021). Modeling Heat Transport in Systems of Hydrate-Filled Sediments Using Residual Thermodynamics and Classical Nucleation Theory. Appl. Sci.11 (9), 4124. 10.3390/app11094124
42
ZhangZ.-C.ChengX.-H. (2016). A Thermo-Mechanical Coupled Constitutive Model for clay Based on Extended Granular Solid Hydrodynamics. Comput. Geotechn.80, 373–382. 10.1016/j.compgeo.2016.05.010
43
ZhangZ.ChengX. (2017). A Fully Coupled THM Model Based on a Non-equilibrium Thermodynamic Approach and its Application. Int. J. Numer. Anal. Meth. Geomech.41, 527–554. 10.1002/nag.2569
44
ZhouH.-f.LiuB.YeF.FuW.-x.TangW.-q.QinY.-d.et al (2021). Landslide Distribution and Sliding Mode Control along the Anninghe Fault Zone at the Eastern Edge of the Tibetan Plateau. J. Mt. Sci.18, 2094–2107. 10.1007/s11629-020-6573-6
Summary
Keywords
GSH theory, thermo-hydro-mechanical response, sandy soil, state variable, stress-strain changes
Citation
Zhou R, Gao S and Wang W (2021) A Granular Thermodynamic Model to Describe the Temperature/Mechanical Characteristics of Sandy Soil. Front. Earth Sci. 9:796523. doi: 10.3389/feart.2021.796523
Received
17 October 2021
Accepted
02 November 2021
Published
09 December 2021
Volume
9 - 2021
Edited by
Xiaodong Fu, Institute of Rock and Soil Mechanics (CAS), China
Reviewed by
Chunyi Cui, Dalian Maritime University, China
Yue Gui, Kunming University of Science and Technology, China
Xiaobing Xu, Zhejiang University of Technology, China
Updates
Copyright
© 2021 Zhou, Gao and Wang.
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.
*Correspondence: Rui Zhou, 18121179@bjtu.edu.cn
This article was submitted to Environmental Informatics and Remote Sensing, a section of the journal Frontiers in Earth Science
Disclaimer
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.