Research Article  Open Access
Micromechanism Underlying Nonlinear StressDependent of Clays at a Wide Range of Pressures
Abstract
In order to investigate the mechanism underlying the reported nonlinear atrest coefficient of earth pressure, of clays at high pressure, a particlescale model which can be used to calculate vertical and horizontal repulsion between clay particles has been proposed. This model has two initial states which represent the clays at low pressure and high pressure, and the particles in this model can undergo rotation and vertical translation. The computation shows that the majority of particles in a clay sample at high pressure state would experience rotation during onedimensional compression. In addition, rotation of particles which tends to form a parallel structure causes an increase of the horizontal interparticle force, while vertical translation leads to a decrease in it. Finally, the link between interparticle force, microstructure, and macroscopic is analyzed and it can be used to interpret well the nonlinear changes in with both vertical consolidation stress and heightdiameter ratio.
1. Introduction
The coefficient of earth pressure at rest, , defined as the ratio of horizontal effective stress to vertical effective stress under the condition of zero horizontal deformation representing the in situ stress state of the ground, is a fundamental parameter in the analysis and design of geotechnical structures.
Numerous studies in the past have indicated that of a given soil is a constant depending on its strength parameter [1]. However, accumulated evidence over recent two decades demonstrated that is not necessarily a constant but generally a function of void ratio, stress level, and critical state friction angle even for given normally consolidated clay [2–5]. In particular, of clay increases nonlinearly with consolidation stresses over a wide range of pressures [6–9].
In addition, previous investigation has indicated that, during onedimensional compression, the distance between clay particles decreases continuously and the orientations of clay particles tend to be parallel to each other with their normal line pointing to the vertical direction [3, 10, 11].
It is well known that there exist noncontact forces such as repulsion between clay particles due to the electric charge on the surface of clay particle, and these interparticle forces usually dominate the mechanical behaviors of clayey material. In addition, interparticle forces, which are balanced with macroscopic stress in the clays, determine the arrangements and orientations (i.e., microscopic structure) of clay particles. The macroscopic stress and deformation can be readily measured during mechanical tests on the clay specimens, but the information related to microscopic structure such as pore size, arrangement, and orientation of clay particles only can be analyzed after stopping the mechanical tests. Since the evolution of microscopic structure of clay during tests is usually unknown, it is difficult to establish the relation between macroscopic mechanical behavior and microscopic structure. Moreover, it is almost impossible to measure the interaction forces between clay particles during macroscopic tests. Therefore, the studies on the links between microscopic structure, interparticle forces, and macroscopic behavior, which is important for the thorough understanding of the intrinsic mechanism of macroscopic mechanical properties of clayey material, are seldom reported. It is not surprising that the micromechanism relating above nonlinear of clay to its microscopic structure during high pressure onedimensional compression has not been exploited.
This study aims to reveal the links between microscopic structure, interparticle forces, and macroscopic using numerical method and get insight into the micromechanism underlying the reported nonlinear of clay at high pressure.
2. Numerical Study on the Links between Interparticle Forces and Microscopic Structure
2.1. Repulsive Forces between Clay Particles
Since it has been revealed that the calculated relations between void ratio and vertical pressure, based on the double layer repulsive forces between clay particles, agree well with the measured results of both low and high pressure onedimensional compression tests [12, 13], the noncontact forces between clay particles in this study are limited to the repulsive forces. The repulsion consists of osmotic pressure and electrical stress . and at arbitrary positon in a system of clayelectrolyte are given as follows [14]:where is the position vector, is the concentration of th ion in the electrolyte, J/(mol·K) is the universal gas constant, is the absolute temperature, is the permittivity of the electrolyte, and is electric field intensity. The repulsive pressure between two clay particles iswhere is placed at any position between two particles and is an arbitrary position outside the region between these two particles.
Both in (1) and in (2) are related to the electric potential:in which is the concentration of background electrolyte, is the electronic unit charge, is the valence of ion in the electrolyte, J/K is Boltzmann’s constant, denotes the gradient operator, and is the electric potential. It should be noted that (4) only applies to binary monovalent electrolytes which is the common case in the study of saturated clay.
In order to calculate the repulsion between clay particles, it is necessary to obtain the potential distribution in the clayelectrolyte system. Dimensionless potential can be determined by solving the following wellknown PoissonBoltzmann equation:whereis the reciprocal of doublelayer thickness in which ’s number.
Combined with the following boundary conditions, the particular solution of (7) in a given domain can be determined. Considerin which is the boundary enclosing the domain and , , , and are constant parameters used to determine the boundary conditions.
There are few analytical solutions of (6) due to its strong nonlinearity, and it is necessary to solve it using numerical method. Using a partial differential equation solver FreeFem++ which is open finite element software [14], (6) can be solved numerically. The first step to solve (6) in this software is to derive its weak form as follows: in which is an arbitrary weight function which is forced to be zero on the boundary . Equation (12) is nonlinear and can be solved using iterative method. Assuming that there is only Dirichlet boundary condition, the following linear iterative equation based on Newton method can be obtained:where the superscript denotes the number of iteration steps.
Figure 1 presents the finite element mesh used for computation of potential around two inclined charged plates with finite length and the corresponding potential nephogram. Substituting the calculated potential into (4), the repulsion between clay particles can be solved.
2.2. Numerical Model for Analysis
Actual microstructure of clay is complex and the amount of particles in a clay sample is huge as shown schematically in Figure 2(a). Because the very fine mesh near the particle surface as shown in Figure 1 is required to accurately calculate the electrical repulsion, it is impractical to simulate actual saturated clay samples. The numerical model shown in Figure 2(b) can be seen as a local point in the clay sample. This particlescale model consists of four clay particles. Each particle can feel both the vertical and horizontal forces from the neighbor particles. At the same time, the effects of macroscopic vertical compression and lateral constraint on the motion of particles can be considered in this model.
(a)
(b)
The particle is 100 nm long and 1 nm thick and the potential on the surface of particle is 248.2 mV, which are the properties of typical montmorillonite clays [15]. The pore space is occupied by binary monovalent electrolyte solution, and the background electrolyte concentration is 1 mol/m^{3}. According to (7)–(10), the corresponding dimensionless length of clay particle and surface potential are 10 and 9.94, respectively.
Two kinds of initial states of the numerical model will be used in this study. The first one represents the clays consolidated at low pressure, and the second is corresponding to the compacted clays at high pressure. According to the reported microscopic structure of clays during onedimensional compression [3, 11], edgetoedge distance dH1 and facetoface distance dV1 between clay particles shown in Figure 2(b) at high pressure are smaller than those at low pressure and the orientation of clay particle shown in Figure 2(b) at high pressure is also smaller than that at low pressure. corresponding to high and low pressure is assumed to be 10° and 30°, respectively. Based on the following relation between pore size and void ratio [16], dH1 and dV1 can be determined. Considerwhere is the specific surface area, is the relative density for deep clay, is the average distance between two neighbor parallel clay plates which represents the pore size, and is void ratio. It is assumed that m^{2}/g and g/m^{3} [13]. Initial void ratios related to high and low pressure are 0.49 and 0.97, respectively, and the calculated initial distances between clay particles using (14) are 3 nm and 6 nm.
The basic models of motion of rigid clay particles are translation and rotation. Due to the effects of macroscopic vertical compression and lateral constraint on the motion of particles, the initial configuration of numerical model shown in Figure 2(b) can change into the two new configurations. The first one is shown by the dotted lines in Figure 3(a), which results from the rotation of particles. The second as shown by the dotted lines in Figure 3(b) results from the vertical translation of particles.
(a)
(b)
2.3. Numerical Results
In order to investigate the microscopic mechanism underlying the difference in of clays at high and low pressure, it is necessary to study the incremental repulsion between particles resulting from the same amount of vertical displacement of clays at both low and high pressure during onedimensional compression. Table 1 presents such numerical results of both vertical and horizontal interparticle repulsion.

It can be seen from Table 1 that, experiencing the same amount of vertical compression in the case of low pressure, the required increase in vertical external forces for pure rotation is 0.3 which is comparable with that for pure translation (0.38). However, in the case of high pressure the required increase in vertical external forces for pure rotation (1.51) is much smaller than that for pure translation (3.89). Therefore, in the case of low pressure, the amount of rotating particles in clay specimens is comparable with that of pure translation. But, in the case of high pressure, the amount of rotation should be larger than that of pure translation.
Another observation from Table 1 is that pure rotation leads to an increase in horizontal interparticle force, but pure translation leads to a decrease in horizontal interparticle force. This holds true for low and high pressure. According to the definition of , it can be deduced that pure rotation would cause an increase in , while pure translation would lead to a decrease in . Apart from the numerical simulations shown above, we carried out a great number of calculations with variables dV1, dH1, and , and similar observations have been obtained.
3. Micromechanism Underlying Nonlinear of Clays over a Wide Range of Pressures
According to the numerical results described in Section 2.3, rotation and translation would cause an increase and a decrease in , respectively, and in the case of low pressure, the amount of clay particles rotating during onedimensional compression is comparable with that of translating particles. Accordingly, the increase in caused by rotation is almost counterbalanced by the decrease in caused by translation. Therefore, the measured during low pressure onedimensional compression of clays largely remains unchanged.
However, in the case of the high pressure, the amount of rotating particles is larger than that of translating particles, and so the increase in caused by rotation is more significant than the decrease in caused by translation. Consequently, the measured should increase with the vertical pressure during high pressure onedimensional compression of clays. This is consistent with the reported test results shown in Figure 4.
Apart from the nonlinear of clays as shown in Figure 4, it has been observed by Min [9] using triaxial testing system that of deep clay in East China at high pressure increases with the decreasing ratio of height to diameter of clay specimens as shown in Figure 5. This test result can also be interpreted well based on the above analysis.
Although it was usually thought that the ideal onedimensional compression tests do not have the problem of end friction, the small lateral deformation of specimens during such tests using triaxial apparatus and so a nonignorable end friction is unavoidable. If there were not such end frictions, of clay specimens in these tests with different heightdiameter ratio would be the same.
The end friction will exert a horizontal shear stress on the ends of the clay specimens as shown in Figure 6 and induces shear stresses inside the specimens. The shear stresses which would produce a torque on clay particles may make it easier to rotate the particles. Therefore, in contrast to those only subjected to vertical stress, the inclined clay particles subjected to both vertical stress and shear stress are more inclined to rotate themselves. According to the early analysis in this section, the increase in during onedimensional compression test with end friction will be more significant than that without end friction. As a result, a higher was observed in the tests with end restraints.
4. Conclusions
In order to explore the micromechanism underlying the reported nonlinear of clays at high pressure, a particlescale numerical model is established and the analysis on the links between particlescale interaction forces, structure, and macroscopic is carried out. The main conclusions of this study are as follows.(1)In the case of low pressure, the amount of purely rotating particles in a clay sample is comparable to that of purely translating particles. But the number of purely rotating particles is larger than that of purely translating particles in the case of high pressure.(2)Pure rotation which tends to make clays reach a parallel configuration would cause an increase in while pure translation which causes smaller facetoface distance between clay particles would lead to a decrease in .(3)Based on the above observations, the links between interparticle forces, microstructure, and macroscopic have been investigated. The links can be used to well interpret the nonlinear changes in of clays at high pressure with both vertical pressure and heightdiameter ratio of clay specimens.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This research is supported by the Major State Basic Research Development Program of China (no. 2012CB026103), the Natural Science Foundation of China (51009136), and the Natural Science Foundation of Jiangsu province (BK2011212).
References
 J. Jaky, “The coefficient of earth pressure at rest,” Journal of the Hungarian Society of Architects and Engineers, vol. 25, p. 355, 1944. View at: Google Scholar
 C. M. R. Ting, G. C. Sills, and D. C. Wijeyesekera, “Development of K_{0} in soft soils,” Geotechnique, vol. 44, no. 1, pp. 101–109, 1994. View at: Publisher Site  Google Scholar
 W.P. Li, Z.Y. Zhang, R.H. Sun, W.L. Wang, and X.Q. Li, “High pressure ${K}_{0}$ creep experiment and the anisotropy of microstructure of deep buried clay,” Chinese Journal of Geotechnical Engineering, vol. 28, no. 10, pp. 1185–1190, 2006. View at: Google Scholar
 N. O. Abdulhadi, J. T. Germaine, and A. J. Whittle, “Stressdependent behavior of saturated clay,” Canadian Geotechnical Journal, vol. 49, no. 8, pp. 907–916, 2012. View at: Publisher Site  Google Scholar
 X. Yao, J. Qi, and F. Yu, “Study on lateral earth pressure coefficient at rest for frozen soils,” Journal of Offshore Mechanics and Arctic Engineering, vol. 136, no. 1, Article ID 011301, 2013. View at: Publisher Site  Google Scholar
 X.Y. Wang, Y.Q. Tang, Y.Z. Zang, J. Chen, and S.P. Han, “Experimental studies and new ideas on the lateral stress in soil,” Chinese Journal of Geotechnical Engineering, vol. 29, no. 3, pp. 430–435, 2007. View at: Google Scholar
 Q. H. Tian, Z. W. Xu, G. Q. Zhou, X. D. Zhao, and K. Hu, “Coefficients of earth pressure at rest in thick and deep soils,” Mining Science and Technology, vol. 19, no. 2, pp. 252–255, 2009. View at: Publisher Site  Google Scholar
 Z.W. Xu, K.H. Zeng, Z. Wei, Z.Q. Liu, X.D. Zhao, and Q.H. Tian, “‘Nonlinear’ characteristics of the static earth pressure coefficient in thick alluvium,” Mining Science and Technology, vol. 19, no. 1, pp. 129–132, 2009. View at: Publisher Site  Google Scholar
 T. Min, Research on test method of earth pressure at rest for deep soil [M.S. thesis], China University of Mining and Technology, Xuzhou, China, 2010, (Chinese).
 A. Anandarajah, “Discreteelement method for simulating behavior of cohesive soil,” Journal of Geotechnical Engineering, vol. 120, no. 9, pp. 1593–1613, 1994. View at: Publisher Site  Google Scholar
 I. DjéranMaigre, D. Tessier, D. Grunberger, B. Velde, and G. Vasseur, “Evolution of microstructures and of macroscopic properties of some clays during experimental compaction,” Marine and Petroleum Geology, vol. 15, no. 2, pp. 109–128, 1998. View at: Publisher Site  Google Scholar
 S. Tripathy and T. Schanz, “Compressibility behaviour of clays at large pressures,” Canadian Geotechnical Journal, vol. 44, no. 3, pp. 355–362, 2007. View at: Publisher Site  Google Scholar
 X.Y. Shang, H.S. Yu, G.Q. Zhou, F. Wang, and Y. Lu, “Micro analysis of mechanical characteristics of deep clay under high stress level,” Chinese Journal of Geotechnical Engineering, vol. 34, no. 2, pp. 363–368, 2012. View at: Google Scholar
 A. J. Grodzinsky, Fields, Forces, and Flows in Biological Systems, Garland Science, Oxford, UK, 2011.
 D. W. Smith, G. A. Narsilio, and P. Pivonka, “Numerical particlescale study of swelling pressure in clays,” KSCE Journal of Civil Engineering, vol. 13, no. 4, pp. 273–279, 2009. View at: Publisher Site  Google Scholar
 X. Y. Shang, G. Q. Zhou, L. F. Kuang, and W. Cai, “Compressibility of deep clay in East China subjected to a wide range of consolidation stresses,” Canadian Geotechnical Journal, vol. 52, no. 2, pp. 244–250, 2015. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Xiangyu Shang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.