Numerical Solution of Desiccation Cracks in Clayey Soils: History
Please note this is an old version of this entry, which may differ significantly from the current revision.

This entry presents the theoretical fundamentals, the mathematical formulation, and the numerical solution for the problem of desiccation cracks in clayey soils. The formulation uses two stress state variables (total stress and suction) and results in a non-symmetric and nonlinear system of transient partial differential equations. A release node algorithm technique is proposed to simulate cracking, and the strategy to implement it in the hydromechanical framework is explained in detail. This general framework was validated with experimental results, and several numerical examples were published at international conferences and in journal papers.

  • desiccation cracking
  • hydromechanical coupling
  • unsaturated soil mechanics

Crack produced by desiccation in clayey soils is a problem that has implications in a wide range of ground-related fields. Geotechnical engineering, agriculture, mining, radioactive waste storage, tailings reservoirs, gravity dams, and public buildings are all affected by cracks due to desiccation [1][2][3][4][5].

The crack patterns that appear during desiccation are random, unique, and their development depends on multiple factors. Experimental approaches and numerical contributions have been made in the last 50 years [2][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22][23]. Comprehensive state of the art was published in the last decade [12][24].

The process of cracking due to desiccation in clayey soils is difficult to treat numerically because the features involved are not well understood and they are dependent on one another . The process of initiation and evolution of the cracks includes a desiccation process first (a hydraulic process) and a process of shrinkage after (a mechanical process). The hydraulic and the mechanical components of the problem are both nonlinear processes. Most of the soil’s properties that affect the cracking process depend on the suction (see definition of suction in Equation (3)) or moisture content values. The process is hydromechanical and coupled. Adding to the complexity is the not-yet-well-understood soil–atmosphere interaction . The mineral composition of the clay as well as the salts dissolved in the water could play a role in the process, but this influence is beyond the scope of this hydromechanical formulation.

In the past, several approaches have been presented for dealing with this problem, with models available for each phase. There are numerical approaches in the literature based on the finite element method (FEM), the discrete element method (DEM), the distinct element method (DiEM), the mesh fragmentation technique (MFT), the lattice spring model (DLSM) [2][11][25][26][27][28][29][30][31], etc.

There is a model [32] that presents a relatively simple application of linear elastic fracture mechanics (LEFM) to simulate tension cracks in soils. In that work, however, the soil is not drying, and the fracture toughness can be determined in a relatively easy way at the laboratory for the convenient water content of the specimen. Commercial codes, such as FLAC® [33], based on the finite difference method have been applied for simulating the curling of a high-plasticity clay subjected to desiccation with no restrictions and with no cracks [11]. Rodríguez et al. [2] developed a one-dimensional hydromechanical approach to study the initiation of a crack during a desiccation process applied to mine wastes. The model proposed by Trabelsi et al. [25] takes into account the heterogeneity of the soil and uses as the main variable the porosity to study the initiation and propagation of cracks in thin samples. The distinct element method proposed by Amarasiri et al. [27] (Amarasiri et al., 2011) is an adaptation of the commercial code UDEC® [34] for analysing desiccation cracks in soils. This approach uses an interface element to simulate the cracks but does not calculate the flow in the porous media (the variation of water content with time needs to be known). The discrete element method (DEM) has been used for the simulation of crack desiccation in thin cylindrical specimens [26] using the commercial code PFC3D® [35] in which the variation of the stiffness and the tensile strength was taken into account and drying shrinkage was imposed on each aggregate at the micro-scale. The mesh fragmentation technique [28] uses high-aspect-ratio elements between the standard elements of a mesh making the model easily adapted to standard finite element programs. The lattice spring model DLSM is a two-phase bond model [29]. Finally, there have been several attempts to use simple models of particles connected by springs or viscoelastic Maxwell elements [36].

The classical theories of continuum mechanics, unsaturated soil mechanics, and the strength of materials and concepts in fracture mechanics are the preferable framework when studying problems such as this. The mathematical problem that emerges from the classical theories in this approach is resolved by the finite element method using material parameters obtained from laboratory tests in a controlled environment. Then, the approach is theoretical, experimental, and numerical (Figure 1).

Figure 1. A theoretical, numerical, and experimental approach to the desiccation crack problem in clayey soils.

In this entry, the flow in deformable porous media is resolved by the finite element method with a  formulation [37]. A release-node technique is used to simulate the cracking process. The capabilities of the approach have been published recently [24][38]. In this approach, desiccation, shrinkage, and cracking are solved without imposing distribution or variation of water content. The system evolves from its initial conditions in suction and displacements first. After a while, the first crack appears when the tensile strength is reached at a point in the soil’s matrix. The release-node algorithm deals with the propagation of the cracks allowing changes in the suction at the boundaries at the new surfaces created by the cracks. Any other technique discussed above to tackle the cracking process can be implemented in the context of the general finite element hydromechanical model proposed here. The more complex the fracture mechanics approach is, the more parameters determined at the laboratory will be necessary and the more numerical instabilities will have to be resolved when simulating the process.

This entry is adapted from the peer-reviewed paper 10.3390/encyclopedia2020068


  1. Rodríguez, R.L. Estudio Experimental de Flujo y Transporte de Cromo, Níquel Y Manganeso en Residuos de la Zona Minera de Moa (Cuba): Influencia del Comportamiento Hidromecánico. Ph.D. Thesis, UPC-BarcelonaTech, Barcelona, Spain, 2002.
  2. Rodríguez, R.L.; Sánchez, M.J.; Ledesma, A.; Lloret, A. Experimental and numerical analysis of a mining waste desiccation. Can. Geotech. J. 2007, 44, 644–658.
  3. Tang, C.-S.; Shi, B.; Liu, C.; Zhao, L.; Wang, B. Influencing factors of geometrical structure of surface shrinkage cracks in clayey soils. Eng. Geol. 2008, 101, 204–217.
  4. Justo, J.L.; Vázquez, N.J.; Justo, E. Subsidence in saturated-unsaturated soils: Application to Murcia (Spain). In Third Interna-tional Conference on Unsaturated Soils; Jucá, J.F.T., de Campos, T.M.P., Marinho, F.A.M., Eds.; Balkema: Recife, Brazil, 2002; pp. 845–850.
  5. Bronswijk, J.J.B. Modeling of water balance, cracking and subsidence in clay soils. J. Hydrol. 1988, 97, 199–212.
  6. Kodikara, J.; Barbour, S.L.; Fredlund, D.G. Desiccation cracking of soil layers. In Unsaturated Soils for Asia; Rahardjo, H., Toll, D.G., Leong, E.C., Eds.; Balkema: Recife, Brazil, 2000; pp. 693–698.
  7. Ávila, G. Estudio de la Retracción Y el Agrietamiento de Arcillas. Aplicación a la Arcilla de Bogotá. Ph.D. Thesis, UPC-BarcelonaTech, Barcelona, Spain, 2004.
  8. Chertkov, V.Y. Modelling cracking stages of saturated soils as they dry and shrink. Eur. J. Soil Sci. 2002, 53, 105–118.
  9. Corte, A.; Higashi, A. Experimental Research on Desiccation Cracks in Soil; U.S. Army Materiel Command Cold Regions Research & Engineering Laboratory: Wilmette, IL, USA, 1960.
  10. Hu, L.B.; Hueckel, T.; Péron, H.; Laloui, L. Modeling Evaporation, Shrinkage and Cracking of Desiccating Soils. In Proceedings of the International Association for Computer Methods and Advances in Geomechanics (IACMAG), Goa, India, 1–6 October, 2008; pp. 1083–1090.
  11. Kodikara, J.K.; Nahlawi, H.; Bouazza, A. Modelling of curling in desiccation clay. Can. Geotech. J. 2004, 41, 560–566.
  12. Lakshmikantha, M.R. Experimental and Theoretical Analysis of Cracking in Drying Soils. Ph.D. Thesis, UPC-BarcelonaTech: Barcelona, Spain, 2009.
  13. Lakshmikantha, M.R.; Prat, P.C.; Ledesma, A. Characterization of crack networks in desiccating soils using image analysis techniques. In Numerical Models in Geomechanics X; Pande, G.N., Pietruszczak, S., Eds.; Balkema: Recife, Brazil, 2007; pp. 167–176.
  14. Lakshmikantha, M.R.; Prat, P.C.; Ledesma, A. Image analysis for the quantification of a developing crack network on a drying soil. Geotech. Test. J. 2009, 32, 505–515.
  15. Lakshmikantha, M.R.; Prat, P.C.; Ledesma, A. Experimental evidences of size-effect in soil cracking. Can. Geotech. J. 2012, 49, 264–284.
  16. Lakshmikantha, M.R.; Prat, P.C.; Ledesma, A. Evidences of hierarchy in cracking of drying soils. ASCE Geotech. Spec. Publ. 2013, 231, 782–789.
  17. Lau, J.T.K. Desiccation Cracking of Clay Soils. M.S. Thesis, University of Saskatchewan, Saskatoon, SK, Canada, 1987.
  18. Morris, P.H.; Graham, J.; Williams, D.J. Cracking in drying soils. Can. Geotech. J. 1992, 29, 263–277.
  19. Nahlawi, H.; Kodikara, J. Laboratory experiments on desiccation cracking of thin soil layers. Geotech. Geol. Eng. 2006, 24, 1641–1664.
  20. Péron, H.; Hueckel, T.; Laloui, L.; Hu, L.B. Fundamentals of desiccation cracking of fine-grained soils: Experimental charac-terisation and mechanisms identification. Can. Geotech. J. 2009, 46, 1177–1201.
  21. Tang, C.-S.; Shi, B.; Liu, C.; Gao, L.; Inyang, H. Experimental Investigation of the Desiccation Cracking Behavior of Soil Layers during Drying. J. Mater. Civ. Eng. 2011, 23, 873–878.
  22. Tang, C.-S.; Shi, B.; Liu, C.; Suo, W.-B.; Gao, L. Experimental characterization of shrinkage and desiccation cracking in thin clay layer. Appl. Clay Sci. 2011, 52, 69–77.
  23. Vogel, H.J.; Hoffmann, H.; Roth, K. Studies of crack dynamics in clay soil. I: Experimental methods, results, and morphological quantification. Geoderma 2005, 125, 203–211.
  24. Levatti, H.U. Estudio Experimental Y Análisis Numérico de la Desecación en Suelos Arcillosos. Ph.D. Thesis, UPC-BarcelonaTech: Barcelona, Spain, 2015. Available online: (accessed on 12 January 2022).
  25. Trabelsi, H.; Jamei, M.; Zenzri, H.; Olivella, S. Crack patterns in clayey soils: Experiments and modelling. Int. J. Numer. Anal. Methods Geomech. 2012, 36, 1410–1433.
  26. Sima, J.; Jiang, M.; Zhou, C. Modeling desiccation cracking in thin clay layer using three-dimensional discrete element method. Proceedings Am. Inst. Phys. 2013, 1542, 245–248.
  27. Amarasiri, A.L.; Kodikara, J.; Costa, S. Numerical modelling of desiccation cracking. Int. J. Numer. Anal. Methods Geomech. 2011, 35, 82–96.
  28. Sánchez, M.; Manzoli, O.L.; Guimarães, L.N. Modeling 3-D desiccation soil crack networks using a mesh fragmentation tech-nique. Comput. Geotech. 2014, 62, 27–39.
  29. Gui, Y.; Zhao, G.F. Modelling of laboratory soil desiccation cracking using DLSM with a two-phase bond model. Comput. Ge-otech. 2015, 69, 578–587.
  30. Yan, C.; Luo, Z.; Zheng, Y.; Ke, W. A 2D discrete moisture diffusion model for simulating desiccation fracturing of soil. Eng. Anal. Bound. Elem. 2022, 138, 42–64.
  31. Nikolic, M. Discrete element model for the failure analysis of partially saturated porous media with propagating cracks rep-resented with embedded strong discontinuities. Comput. Methods Appl. Mech. Eng. 2022, 390, 27.
  32. Lee, F.H.; Lo, K.W.; Lee, S.L. Tension crack development in soils. J. Geotech. Eng. ASCE 1988, 114, 915–929.
  33. ITASCA. FLAC3D: Fast Lagrangian Analysis of Continua in Three Dimensions, Version 2.0 Manual; Itasca Consulting Group: Minneapolis, MN, USA, 1997.
  34. ITASCA. UDEC: Distinct-Element Modeling of Jointed and Blocky Material in 2D; Itasca Consulting Group: Minneapolis, MN, USA, 2017.
  35. ITASCA. PFC3D: Particle Flow Code in Three Dimensions, Version 3.0 Manual; Itasca Consulting Group: Minneapolis, MN, USA, 2003.
  36. Musielak, G.; Sliwa, T. Modeling and Numerical Simulation of Clays Cracking During Drying. Dry. Technol. 2015, 33, 1758–1767.
  37. Zienkiewicz, O.C.; Chan, A.H.C.; Pastor, M.; Paul, D.K.; Shiomi, T. Static and Dynamic Behaviour of Soils: A Rational Ap-proach to Quantitative Solutions. I. Fully Saturated Problems. Proceedings of the Royal Society of London Series A. Math. Phys. Sci. 1990, 429, 285–309.
  38. Levatti, H.U.; Prat, P.C.; Ledesma, A. Numerical and experimental study of initiation and propagation of desiccation cracks in clayey soils. Comput. Geotech. 2019, 105, 155–167. Available online: (accessed on 12 January 2022).
This entry is offline, you can click here to edit this entry!