Abstract: Clay minerals can experience strong tensile and compressive forces in extreme environments such as the deep sea and subsurface. Moreover, the presence of water films greatly affects the mechanical properties of clay. To explore these properties, we use a molecular dynamics (MD) simulation method to study axial mechanical behavior and failure mechanisms of hydrated kaolinite. Two types of deformation are applied to kaolinite examples with varying water film thicknesses: stretching along the transverse (x) direction, and compression along the longitudinal (z) direction. The ultimate strengths of hydrated kaolinite with different water film thicknesses range from 8.12% to 27.53% (for stretching along the x-direction) and from 15.71% to 26.02% (for compression along the z-direction) less than those of dehydrated kaolinite. Additionally, we find that hydrated kaolinite is more prone to tensile than compressive failure under high stress. When stretched along the x-direction, the diffusion of water molecules results in unstable tensile properties. When compressed along the z-direction, water films weaken the compressive strength of the system and lead to greater compressive deformation, but also delay the time at which the system fails. Furthermore, we investigated the failure mechanisms of hydrated kaolinite through analysis of interaction energies. The tensile failure along the x-direction is caused by the breaking of the covalent bonds in the clay mineral sheet. On the other hand, the compressive failure along the z-direction is due to the crushing of the internal structure of the clay mineral sheet.

含水高岭石拉伸和压缩性能的分子动力学模拟研究

[1]ArmandG, NoiretA, ZghondiJ, et al., 2013. Short- and long-term behaviors of drifts in the Callovo-Oxfordian claystone at the Meuse/Haute-Marne Underground Research Laboratory. Journal of Rock Mechanics and Geotechnical Engineering, 5(3):221-230.

[2]BencoL, TunegaD, HafnerJ, et al., 2001. Upper limit of the O‒H···O hydrogen bond. Ab initio study of the kaolinite structure. The Journal of Physical Chemistry B, 105(44):10812-10817.

[3]ChenJ, MinFF, LiuLY, et al., 2019. Mechanism research on surface hydration of kaolinite, insights from DFT and MD simulations. Applied Surface Science, 476:6-15.

[4]CyganRT, LiangJJ, KalinichevAG, 2004. Molecular models of hydroxide, oxyhydroxide, and clay phases and the development of a general force field. The Journal of Physical Chemistry B, 108(4):1255-1266.

[5]Duque-RedondoE, ManzanoH, Epelde-ElezcanoN, et al., 2014. Molecular forces governing shear and tensile failure in clay-dye hybrid materials. Chemistry of Materials, 26(15):4338-4345.

[6]FanQC, WangZP, MengX, et al., 2022. Multi-scale analysis of the strengthening mechanism of functionalized graphene as reinforcement in cement composites. Colloids and Surfaces A: Physicochemical and Engineering Aspects, 651:129729.

[7]FanQC, ZhengYY, LiuYM, et al., 2023. Effect of modified cellulose nanocrystals on the structure of calcium silicate hydrate studied by molecular dynamics simulation and experiment. Langmuir, 39(46):16244-16260.

[8]GengZC, TangSW, WangY, et al., 2024. Stress relaxation properties of calcium silicate hydrate: a molecular dynamics study. Journal of Zhejiang University-SCIENCE A (Applied Physics & Engineering), 25(2):97-115.

[9]GonzálezNA, RouainiaM, ArroyoM, et al., 2012. Analysis of tunnel excavation in London clay incorporating soil structure. Géotechnique, 62(12):1095-1109.

[10]HooverWG, 1985. Canonical dynamics: equilibrium phase-space distributions. Physical Review A, 31(3):1695-1697.

[11]JorgensenWL, ChandrasekharJ, MaduraJD, et al., 1983. Comparison of simple potential functions for simulating liquid water. The Journal of Chemical Physics, 79(2):926-935.

[12]KalinichevAG, WangJW, KirkpatrickRJ, 2007. Molecular dynamics modeling of the structure, dynamics and energetics of mineral-water interfaces: application to cement materials. Cement and Concrete Research, 37(3):337-347.

[13]LengYS, CummingsPT, 2006. Hydration structure of water confined between mica surfaces. The Journal of Chemical Physics, 124(7):074711.

[14]LiuXD, LuXC, WangRC, et al., 2012. Atomic scale structures of interfaces between kaolinite edges and water. Geochimica et Cosmochimica Acta, 92:233-242.

[15]LuM, DiaoQF, ZhengYY, et al., 2024a. Influence of water on tensile behavior of Illite through the molecular dynamics method. International Journal of Geomechanics, 24(4):04024024.

[16]LuM, HuangW, ZhengYY, 2024b. Molecular dynamics simulation study of the effects of water content and wettability on the shear properties of kaolinite for the failure of clay soil. ACS Applied Nano Materials, 7(3):2843-2854.

[17]LuM, ZhengYY, YinZY, 2024c. From sedimentation to consolidation of kaolinite: a molecular dynamic study. Computers and Geotechnics, 170:106285.

[18]LuM, ZhengYY, YinZY, 2025. Investigating the consolidation of kaolinite with molecular dynamics: the micro effective stress principle. Applied Surface Science, 690:162653.

[19]MaZY, Pathegama GamageR, RathnaweeraT, et al., 2019. Review of application of molecular dynamic simulations in geological high-level radioactive waste disposal. Applied Clay Science, 168:436-449.

[20]MartinRT, 1960. Adsorbed water on clay: a review. Clays and Clay Minerals (National Conference on Clays and Clay Minerals), 9:28-70.

[21]MilheiroFAC, FreireMN, SilvaAGP, et al., 2005. Densification behaviour of a red firing Brazilian kaolinitic clay. Ceramics International, 31(5):757-763.

[22]Moreno-MarotoJM, Alonso-AzcárateJ, 2015. An accurate, quick and simple method to determine the plastic limit and consistency changes in all types of clay and soil: the thread bending test. Applied Clay Science, 114:497-508.

[23]NederRB, BurghammerM, GraslT, et al., 1999. Refinement of the kaolinite structure from single-crystal synchrotron data. Clays and Clay Minerals, 47(4):487-494.

[24]PlimptonS, 1995. Fast parallel algorithms for short-range molecular dynamics. Journal of Computational Physics, 117(1):1-19.

[25]PorionP, MichotLJ, FaugèreAM, et al., 2007. Structural and dynamical properties of the water molecules confined in dense clay sediments: a study combining ²H NMR spectroscopy and multiscale numerical modeling. The Journal of Physical Chemistry C, 111(14):5441-5453.

[26]PouvreauM, GreathouseJA, CyganRT, et al., 2019. Structure of hydrated kaolinite edge surfaces: DFT results and further development of the ClayFF classical force field with metal‍‒‍O‍‒‍H angle bending terms. The Journal of Physical Chemistry C, 123(18):11628-11638.

[27]SmirnovKS, BougeardD, 1999. A molecular dynamics study of structure and short-time dynamics of water in kaolinite. The Journal of Physical Chemistry B, 103(25):5266-5273.

[28]ŠolcR, GerzabekMH, LischkaH, et al., 2011. Wettability of kaolinite (001) surfaces—molecular dynamic study. Geoderma, 169:47-54.

[29]SpagnoliG, SridharanA, OresteP, et al., 2017. A probabilistic approach for the assessment of the influence of the dielectric constant of pore fluids on the liquid limit of smectite and kaolinite. Applied Clay Science, 145:37-43.

[30]SunHM, YangW, ChenRP, et al., 2021. Microfabric characteristics of kaolinite flocculates and aggregates—insights from large-scale molecular dynamics simulations. Applied Clay Science, 206:106073.

[31]SuterJL, CoveneyPV, GreenwellHC, et al., 2007. Large-scale molecular dynamics study of montmorillonite clay: emergence of undulatory fluctuations and determination of material properties. The Journal of Physical Chemistry C, 111(23):8248-8259.

[32]TangSW, A HB, ChenJT, et al., 2019. The interactions between water molecules and C-S-H surfaces in loads-induced nanopores: a molecular dynamics study. Applied Surface Science, 496:143744.

[33]TararushkinEV, PisarevVV, KalinichevAG, 2022. Atomistic simulations of ettringite and its aqueous interfaces: structure and properties revisited with the modified ClayFF force field. Cement and Concrete Research, 156:106759.

[34]TunegaD, GerzabekMH, LischkaH, 2004. Ab initio molecular dynamics study of a monomolecular water layer on octahedral and tetrahedral kaolinite surfaces. The Journal of Physical Chemistry B, 108(19):5930-5936.

[35]WeiPC, ZhengYY, ZaouiA, et al., 2023. Atomistic study on thermo-mechanical behavior and structural anisotropy of montmorillonite under triaxial tension and compression. Applied Clay Science, 233:106817.

[36]WongsarojJ, SogaK, MairRJ, 2013. Tunnelling-induced consolidation settlements in London clay. Géotechnique, 63(13):1103-1115.

[37]YangH, HeMC, LuCS, et al., 2019. Deformation and failure processes of kaolinite under tension: insights from molecular dynamics simulations. Science China Physics, Mechanics & Astronomy, 62(6):64612.

[38]ZhangCL, 2018. Thermo-hydro-mechanical behavior of clay rock for deep geological disposal of high-level radioactive waste. Journal of Rock Mechanics and Geotechnical Engineering, 10(5):992-1008.

[39]ZhangLH, LuXC, LiuXD, et al., 2014. Hydration and mobility of interlayer ions of (Na_x, Ca_y)-montmorillonite: a molecular dynamics study. The Journal of Physical Chemistry C, 118(51):29811-29821.

[40]ZhangLL, ZhengYY, WeiPC, et al., 2021. Nanoscale mechanical behavior of kaolinite under uniaxial strain conditions. Applied Clay Science, 201:105961.

[41]ZhangZ, SongXY, 2022. Molecular dynamics modeling of cracks in dry clay sheets at the nanoscale. Computers and Geotechnics, 152:105037.

[42]ZhaoMH, ZouXJ, ZouPX, 2007. Disintegration characteristics of red sandstone and its filling methods for highway roadbed and embankment. Journal of Materials in Civil Engineering, 19(5):404-410.