An orthotropic interface damage model for simulating drying processes in soils

The study of drying process in soils has received an increased attention in the last few years. This is very complex phenomenon that generally leads to the formation and propagation of desiccation cracks in the soil mass. In recent engineering applications, high aspect ratio elements have proved to be well suited to tackle this type of problem using finite elements. However, the modeling of interfaces between materials with orthotropic properties that generally exist in this type of problem using standard (isotropic) constitutive model is very complex and challenging in terms of the mesh generation, leading to very fine meshes that are intensive CPU demanding. A novel orthotropic interface mechanical model based on damage mechanics and capable of dealing with interfaces between materials in which the strength depends on the direction of analysis is proposed in this paper. The complete mathematical formulation is presented together with the algorithm suggested for its numerical implementation. Some simple yet challenging synthetic benchmarks are analyzed to explore the model capabilities. Laboratory tests using different textures at the contact surface between materials were conducted to evaluate the strengths of the interface in different directions. These experiments were then used to validate the proposed model. Finally, the approach is applied to simulate an actual desiccation test involving an orthotropic contact surface. In all the application cases the performance of the model was very satisfactory.

[1]  A. Misra Micromechanical model for anisotropic rock joints , 1999 .

[2]  Alberto Ledesma,et al.  Experimental and numerical analysis of desiccation of a mining waste , 2007 .

[3]  M. Zieliński,et al.  Precise observation of soil surface curling , 2014 .

[4]  J. C. Simo,et al.  An analysis of strong discontinuities induced by strain-softening in rate-independent inelastic solids , 1993 .

[5]  J. Oliver,et al.  Strong discontinuities and continuum plasticity models: the strong discontinuity approach , 1999 .

[6]  Esteban Samaniego,et al.  AES for multiscale localization modeling in granular media , 2011 .

[7]  Osvaldo L. Manzoli,et al.  Modeling of interfaces in two-dimensional problems using solid finite elements with high aspect ratio , 2012 .

[8]  Osvaldo L. Manzoli,et al.  Modeling 3-D desiccation soil crack networks using a mesh fragmentation technique , 2014 .

[9]  Michael Ortiz,et al.  An Efficient Adaptive Procedure for Three-Dimensional Fragmentation Simulations , 2001, Engineering with Computers.

[10]  Osvaldo L. Manzoli,et al.  On the use of finite elements with a high aspect ratio for modeling cracks in quasi-brittle materials , 2016 .

[11]  Yujing Jiang,et al.  Anisotropic shear behavior of closely jointed rock masses , 2014 .

[12]  Ignacio Carol,et al.  A meso‐level approach to the 3D numerical analysis of cracking and fracture of concrete materials , 2006 .

[13]  Kerstin Weinberg,et al.  A numerical approach to the analysis of failure modes in anisotropic plates , 2011 .

[14]  J. Kodikara,et al.  Laboratory experiments on desiccation cracking of thin soil layers , 2006 .

[15]  Pere C. Prat,et al.  Image Analysis for the Quantification of a Developing Crack Network on a Drying Soil , 2009 .

[16]  Alejandro Mota,et al.  Fracture and fragmentation of simplicial finite element meshes using graphs , 2006 .

[17]  Ignacio Carol,et al.  NORMAL/SHEAR CRACKING MODEL: APPLICATION TO DISCRETE CRACK ANALYSIS , 1997 .

[18]  Antonio Aguado,et al.  Meso-structural study of concrete fracture using interface elements. I: numerical model and tensile behavior , 2008 .

[19]  Rashid K. Abu Al-Rub,et al.  Prediction of micro and nanoindentation size effect from conical or pyramidal indentation , 2007 .

[20]  H. Moslemi,et al.  Modeling of cohesive crack growth using an adaptive mesh refinement via the modified-SPR technique , 2009 .

[21]  Antonio Aguado,et al.  Meso-structural study of concrete fracture using interface elements. II: compression, biaxial and Brazilian test , 2008 .

[22]  Lyesse Laloui,et al.  Fundamentals of desiccation cracking of fine-grained soils: experimental characterisation and mechanisms identification , 2009 .

[23]  Jonny Rutqvist,et al.  Hydro-mechanical model for wetting/drying and fracture development in geomaterials , 2014, Comput. Geosci..

[24]  Michael Kaliske,et al.  Discrete crack path prediction by an adaptive cohesive crack model , 2010 .

[25]  George A. Gazonas,et al.  The cohesive element approach to dynamic fragmentation: the question of energy convergence , 2007 .

[26]  Luís A.G. Bitencourt,et al.  An adaptive concurrent multiscale model for concrete based on coupling finite elements , 2018 .

[27]  A. Huespe,et al.  Stability and robustness issues in numerical modeling of material failure with the strong discontinuity approach , 2006 .

[28]  Chun Liu,et al.  Automatic quantification of crack patterns by image processing , 2013, Comput. Geosci..

[29]  Xiaopeng Xu,et al.  Numerical simulations of fast crack growth in brittle solids , 1994 .

[30]  A. Pandolfi,et al.  Numerical modelling of fracture in human arteries , 2008, Computer methods in biomechanics and biomedical engineering.

[31]  Enzo Martinelli,et al.  Zero-thickness interface model formulation for failure behavior of fiber-reinforced cementitious composites , 2012 .

[32]  Richard A. Regueiro,et al.  A finite element model of localized deformation in frictional materials taking a strong discontinuity approach , 1999 .

[33]  A. Caggiano,et al.  Multiscale failure analysis of fiber reinforced concrete based on a discrete crack model , 2012, International Journal of Fracture.

[34]  Jayantha Kodikara,et al.  Numerical modelling of desiccation cracking , 2011 .

[35]  Rong Tian,et al.  Improved XFEM: Accurate and robust dynamic crack growth simulation , 2016 .

[36]  D. Seyedi,et al.  FE modelling with strong discontinuities for 3D tensile and shear fractures: Application to underground excavation , 2016 .

[37]  G. Miller,et al.  Shear strength of unsaturated soil interfaces , 2009 .

[38]  Osvaldo L. Manzoli,et al.  Modeling of multiple cracks in reinforced concrete members using solid finite elements with high aspect ratio , 2014 .

[39]  Jonny Rutqvist,et al.  TOUGH-RDCA modeling of multiple fracture interactions in caprock during CO2 injection into a deep brine aquifer , 2014, Comput. Geosci..

[40]  Sewon Kim,et al.  Exploring desiccation cracks in soils using a 2D profile laser device , 2013 .

[41]  Sutthisak Phongthanapanich,et al.  Adaptive Delaunay triangulation with object-oriented programming for crack propagation analysis , 2004 .

[42]  A. Gens,et al.  AN INTERFACE ELEMENT FORMULATION FOR THE ANALYSIS OF SOIL-REINFORCEMENT INTERACTION , 1989 .