Loss of uniqueness of numerical solutions of the borehole problem modelled with enhanced media

It is well known that, initial boundary value problems involving constitutive equations modeling the degradation of the strength of materials are not well posed, which renders computations questionable. To overcome this issue it is necessary to enhance the models by incorporating some internal length. It has been shown that such an enhancement restores the objectivity of the computation as spurious mesh dependency is avoided. However, at least for simple problems (e.g. one dimensional ones), it has been proven that uniqueness of the underlying mathematical problem itself is not restored. Moreover numerical modeling of element tests yields several solutions. This paper demonstrates that several numerical solutions can be obtained also for less simple problems, namely the borehole problems. Even when a defect is introduced in the computed problems, different numerical solutions are found. Contrary to the one dimensional problem there is no proof that this loss of uniqueness comes from the underlying mathematical problem. It is our opinion that this is an inherent property of initial boundary value problems where, broadly speaking, strong degradation of the mechanical properties is modeled. In any case, it is necessary to be aware of this issue. For problems involving constitutive equation modeling strength degradation, it is important to try to find other solutions than the one obtained by using routinely a numerical code. The failure patterns of the different solutions found are however similar to experimental observations. This possible loss of uniqueness can then be seen as a counterpart of the difficulties encountered when attempting to reproduce experiments. This is crucial when dealing with geomaterials.

[1]  René Chambon,et al.  Uniqueness studies in boundary value problems involving some second gradient models , 2004 .

[2]  Frédéric Collin,et al.  SWITCHING DEFORMATION MODES IN POST-LOCALIZATION SOLUTIONS WITH A QUASIBRITTLE MATERIAL , 2006 .

[3]  E. Aifantis On the Microstructural Origin of Certain Inelastic Models , 1984 .

[4]  R. Chambon,et al.  Etude de la localisation unidimensionnelle à l'aide d'un modèle de second gradient , 1996 .

[5]  S. Xia,et al.  A nonlocal damage theory , 1987 .

[6]  R. D. Mindlin Second gradient of strain and surface-tension in linear elasticity , 1965 .

[7]  I. Vardoulakis,et al.  Bifurcation analysis of deep boreholes: II. Scale effect , 1989 .

[8]  C. Truesdell,et al.  The Non-Linear Field Theories Of Mechanics , 1992 .

[9]  S. Radi Etude numérique de la localisation à l'aide de modèles de second gradient : Perte d'unicité et évolution de la taille de la zone localisée , 2005 .

[10]  P. Germain,et al.  The Method of Virtual Power in Continuum Mechanics. Part 2: Microstructure , 1973 .

[11]  R. D. Mindlin Micro-structure in linear elasticity , 1964 .

[12]  René Chambon,et al.  A strain space gradient plasticity theory for finite strain , 2004 .

[13]  Robert Charlier,et al.  An algorithm and a method to search bifurcation points in non‐linear problems , 2001 .

[14]  R. Chambon,et al.  Void ratio evolution inside shear bands in triaxial sand specimens studied by computed tomography , 1996 .

[15]  I. Vardoulakis,et al.  Bifurcation Analysis in Geomechanics , 1995 .

[16]  René Chambon,et al.  Plastic continuum with microstructure, Local Second Gradient Theories for Geomaterials : Localization Studies , 2001 .

[17]  Gioacchino Viggiani,et al.  Volumetric Digital Image Correlation Applied to X‐ray Microtomography Images from Triaxial Compression Tests on Argillaceous Rock , 2007 .

[18]  D. Caillerie,et al.  One-dimensional localisation studied with a second grade model , 1998 .

[19]  Jian Guo Yu,et al.  Computational modelling of the localized deformation associated with borehole breakout in quasi-brittle materials , 2003 .

[20]  René Chambon,et al.  Large strain finite element analysis of a local second gradient model: application to localization , 2002 .