Experimental and theoretical analysis of cracking in drying soils

The thesis focuses on the experimental and theoretical aspects of the process of cracking in drying soils. The results and conclusions were drawn from an exhaustive experimental campaign characterised by innovative multidisciplinary aspects incorporating Fracture Mechanics and classical Soil mechanics, aided with image analysis techniques. A detailed study of the previous works on the topic showed the absence of large scale fully monitored laboratory tests, while the existing studies were performed on already cracked soil, rather than focusing on the process of cracking. Also, the absence of a multidisciplinary approach was strongly felt. The soil used in the present work is a Barcelona silty clay, extensively studied previously for its Thermo-Hydro-Mechanical (THM) behaviour. Tensile strength and fracture toughness are two important fracture parameters that were not characterized previously and which were determined for the soil used in the investigation. The effect of moisture content on tensile strength and fracture toughness was established. The variation of tensile strength with the degree of saturation was theoretically explored considering the concept of unsaturated cohesion, while the concepts of activation energy and rate process theory were used to explain the relation between the fracture toughness and moisture content. Preliminary desiccation tests were conducted to gain basic knowledge of the cracking behaviour of Barcelona soil. Further experiments were conducted in a controlled laboratory environment with rectangular holding trays of similar geometry with five different surface areas. Crack initiation, temporal evolution of cracks resulting in the final crack pattern, meeting and bifurcation of cracks resulting in intersections, etc. were studied in detail. One of the objectives of the desiccation tests with similar geometry was to check the applicability of fracture mechanics. The results of the experiments showed existence of size-effect in soil cracking, thus justifying the use of fracture mechanics as a framework to model crack propagation. A novel equipment capable of simulating the different combinations of climatic conditions with monitoring by sensors and a digital camera was designed and constructed. This enivironmental chamber can hold large specimens up to 80 cm in diameter and 20 cm thick. The results from experiments carried out with the environmental chamber have provided new knowledge about the process and mechanisms of desiccation cracking. Some observations made during the experiments prompted more detailed experiments on some aspects of crack formation. The morphology of cracks, fissures, and the presence of spiral and ripple-like cracks on the bottom surface demanded to explore a different theoretical point of view into the mechanism of cracking in soils. The possibility of using a theory based on either classical or unsaturated soil mechanics is critically examined, with emphasis on the importance of surrounding ambient conditions and seasonal variations. Finally, detailed macro-morphology analysis of the crack patterns was performed. A theory was developed based on the morphology of crack patterns to explain the crack formation as a result of successive domain divisions or hierarchical nature of crack formation. To test the proposed theory, temporal evolution of cracks and crack pattern were studied in detail, showing the existence of hierarchy very clearly. An attempt has been made to explore a new way of analysing the cracking in soils as a process in transition from a disordered (at crack initiation) to a deterministic behaviour (at the end of drying resulting in a complex pattern). Finally, the possible mechanisms for the formation of cracks at the bottom of the specimens were examined, with emphasis on curling and synaeresis processes.La tesis se centra en los aspectos teóricos y experimentales del proceso de agrietamiento de suelos por desecación. Los resultados y las conclusiones fueron extraídos de una exhaustiva campaña de experimentación caracterizada por su aspecto multidisciplinario, incorporando las bases de la mecánica de fractura y la mecánica del suelo con la ayuda de técnicas de análisis de imagen. Un estudio detallado de los trabajos anteriores sobre el tema puso de manifiesto la ausencia de experimentos a gran escala monitorizados en el laboratorio. Los estudios existentes se han realizado sobre suelos agrietados en lugar de centrarse en suelos en proceso de agrietamiento. Es remarcable, además, la ausencia de un enfoque multidisciplinario. El suelo utilizado en el presente trabajo es una arcilla limosa de Barcelona, ampliamente estudiada por su comportamiento Termo-Hidro-Mecánico (THM). Se recopilaron las características geológicas y geotécnicas reportadas en anteriores tesis doctorales y se determinaron dos nuevos importantes parámetros de fractura tales como la resistencia a tracción y la tenacidad de fractura del suelo estudiado. Fueron utilizados equipos específicos para cada parámetro. Se estableció el efecto del contenido de humedad en el suelo sobre la resistencia a tracción y tenacidad de fractura. Se estudió de forma teórica la variación de la resistencia a tracción con el grado de saturación, teniendo en cuenta el concepto de cohesión no saturada. Los conceptos de energía de activación y teoría de tasa de proceso se utilizaron para explicar la relación entre la tenacidad de fractura y el contenido de humedad. Para conocer el comportamiento del suelo estudiado en cuanto a desecación y agrietamiento se realizaron unos ensayos preliminares. Adicionalmente, se realizó otra serie de experimentos con bandejas de geometría similar pero con diferente superficie, estudiando con precisión detalles como la iniciación y la evolución temporal de las grietas, sus intersecciones y el patrón final de agrietamiento. En general los resultados de los ensayos de desecación confirmaron los datos conocidos, excepto algunas observaciones que no correspondían a ningún resultado de estudios previos publicados. Uno de los objetivos de los ensayos de desecación con geometrías similares fue el de estudiar la aplicabilidad de la teoría de Mecánica de Fractura. Los experimentos demostraron la existencia del efecto del tamaño en el agrietamiento del suelo, justificando el uso de la mecánica de fractura como marco teórico de referencia para modelar el proceso de propagación de las fisuras en suelos. Se diseñó y construyó un nuevo equipo especial para el estudio del fenómeno del agrietamiento por desecación, la cámara ambiental, capaz de simular las diferentes combinaciones de condiciones climáticas con control por sensores y una cámara fotográfica digital integrada. En ella se pudieron estudiar muestras de gran tamaño, hasta 80 cm de diámetro y 20 cm de espesor. Algunas observaciones llevadas a cabo durante los experimentos mostraron la necesidad de realizar estudios ms detallados y de cambiar el punto de vista teórico del mecanismo de agrietamiento de suelos. Especialmente, teniendo en cuenta la presencia de grietas en forma de espiral en la superficie inferior de la muestra. La posibilidad de utilizar teorías basadas tanto en la mecánica de suelos clásica como en la no saturada fue examinada desde un punto de vista crítico, con énfasis en la importancia de las condiciones ambientales y las variaciones estacionales. Por último, se llevó a cabo un estudio sobre la macro-morfología detallada del patrón de agrietamiento. Una teoría basada en la morfología de las grietas fue desarrollada para explicar la formación de las grietas como resultado de las sucesivas divisiones de áreas integrales o de la naturaleza jerárquica de la formación de grietas. Para probar la nueva teoría, se estudió la evolución temporal de las grietas y el patrón de fisuras en detalle y los resultados mostraron la existencia de una jerarquía clara. Se realizó un intento de explorar una nueva forma de enfocar el agrietamiento de suelos como una conducta desordenada o determinista. Finalmente, los posibles mecanismos para la formación de grietas en la parte inferior de la muestra fueron examinados, dando importancia a los fenómenos de deformación por flexión y sinéresis.La tesis es centra en els aspectes teòrics i experimentals del procés d’esquerdament en sòls per dessecació. Els resultats i les conclusions varen ser extretes d’una exhaustiva campanya d’experimentació caracteritzada pel seu aspecte multidisciplinari, incorporant les bases de la mecànica de fractura i la mecànica del sòl amb l’ajuda de tècniques d’anàlisis d’imatge. Un estudi detallat dels treballs anteriors sobre el tema va posar de manifest l’absència d’experiments a gran escala monitoritzats al laboratori. Els estudis existents s’han realitzat sobre sols esquerdats en comptes de centrar-se en sòls en procés d’esquerdament. És remarcable, a més, l’absència d’un enfocament multidisciplinar. El sòl utilitzat en el present treball és una argila llimosa de Barcelona, àmpliament estudiada pel seu comportament Termo-Hidro-Mecànic (THM). Es varen recopilar les característiques geològiques i geotècniques descrites en anteriors treballs i es varen determinar dos nous importants paràmetres de fractura tals com la resistència a tracció i la tenacitat de fractura del sòl estudiat. Van ser necessaris equips específics per a cada un dels paràmetres. Es va establir l’efecte del contingut d’humitat del sòl sobre la resistència a tracció i la tenacitat de fractura. Es va explorar de forma teòrica la variació de la resistència a tracció amb el grau de saturació, tenint en compte el concepte de cohesió no saturada. Els conceptes d’energia d’activació i la teoria de velocitat dels processos es varen utilitzar per a explicar la relació entre la tenacitat de fractura i el contingut d’humitat. Per a conèixer el comportament del sòl estudiat en quan a dessecació i esquerdament es varen realitzar uns assaigs preliminars. Addicionalment, es varen realitzar una altra sèrie d’experiments amb safates de geometria similar però amb diferents superfícies, estudiant amb precisió detalls com l’inici i l’evolució temporal de les esquerdes, llurs interseccions i el patró final d’esquerdament. En general els resultats dels assaigs de dessecació varen confirmar les dades conegudes, exceptuant algunes observacions que no corresponien a cap resultat d’estudis previs publicats. Un dels objectius dels assaigs de dessecació amb geometries similars va ser el d’estudiar l’aplicabilitat de la teoria de mecànica de fractura. Els experiments varen demostrar l’existència del efecte del tamany en l’esquerdament de sòls, justificant l’ús de la teoria de la mecànica de fractura com a marc de referència per a un model teòric de propagació d’esquerdes en sòls. Es va dissenyar i construir un nou equip especial per a l’estudi del fenomen d’esquerdament per dessecació, la cambra ambiental, capaç de simular diferents combinacions de condicions climàtiques amb control per sensors i una càmera digital integrada. Amb ella es varen poder estudiar mostres de gran tamany, fins a 80 cm de diàmetre i 20 cm d’espessor. Algunes observacions dutes a terme durant els experiments varen mostrar la necessitat de realitzar estudis més detallats i de canviar el punt de vista teòric del mecanisme de d’esquerdament en sòls. Especialment, tenint en compte la presència d’esquerdes en forma d’espiral a la superfície inferior de la mostra. La possibilitat d’utilitzar teories basades tant en la mecànica de sòls clàssica com en la no saturada va ser examinada des d’un punt de vista crític, emfatitzant la importància de les condicions ambientals i les variacions estacionals. Per últim, es va estudiar la macro-morfologia detallada del patró d’esquerdament. Una teoria basada en la morfologia de les esquerdes va ser desenvolupada per explicar la formació de les esquerdes com a resultat de les successives divisions d’àrees integrals o de la naturalesa jeràrquica de la formació d’esquerdes. Per a provar la nova teoria, es va estudiar l’evolució temporal de les esquerdes i el patró d’esquerdament en detall i les resultats varen mostrar l’existència d’una jerarquia clara. Es va realitzar un intent d’explorar una nova manera d’observar l’esquerdament de sòls com una conducta desordenada o determinista. Finalment, els possibles mecanismes per a la formació d’esquerdes a la part inferior de la mostra varen ser examinats, donant importància als fenòmens de deformació per flexió i la sinèresi

Coupled h-m fracture interaction using fem with zero-thickness interface elements

Intensive hydraulic fracturing is a procedure employed for low permeability reservoir stimulation. This technique consists of generating a sequence of regularly spaced parallel fractures (multi-stage fracturing). The generation of a fracture involves the modification of the local stress state, and therefore, in the case of multi-stage fracturing, the propagation of a certain fracture can be affected by the injection sequence, as it has been observed with microseismicity monitoring [1]. This paper describes a study of this technique by means of the Finite Element Method with zero-thickness interface elements for the geo-mechanical modelling of discontinuities [2]. The technique consists in inserting interface elements in between standard elements to allow jumps in the displacement solution fields. For the mechanical problem, their kinematic constitutive variables are relative displacements, and the corresponding static variables are stress tractions. The relationship between variables is controlled via a fracture-based constitutive law with elasto-plastic structure [3]. Concerning the hydraulic problem, the interface formulation includes both the longitudinal flow (with a longitudinal conductivity parameter strongly dependent on the fracture aperture), as well as and the transversal flow across the element [4]. Previous work by the authors focused on the validation of the method, the analysis a single fracture plane problem [5, 6]. In this case the method is extended to allow free propagation of fractures in any direction, by means of inserting interface elements between all continuum elements. The results presented in this paper analyse the effect of material properties, in particular fracture characterization, in the propagation and the effect of different major to minor principal horizontal stress ratio, on the trajectory and interaction of the fractures

Avoiding fracture instability in wedge splitting tests by means of numerical simulations

In this paper, unstable fracture propagation obtained in a in-house performed experimental Wedge Splitting Test (WST) is simulated by means of the FEM and fracturebased zero-thickness interface elements. In order to obtain a specimen geometry suitable for a stable WST without modifying the remaining significant parameters of the test (machine stiffness and control parameter), additional simulations were performed varying the length of the specimen notch, until a load-COD (Crack Opening Displacement) curve without snap-back was obtained. Finally, a new experimental WST with the modified geometry was carried out leading to a stable load-COD curve. In the simulations, elastic continuum elements were used to represent the rock, the steel loading plates and the testing machine compliance via an “equivalent spring”, whereas interface elements were used for the notch and along the potential crack path. The interface elements representing the notch were equipped with linear elastic constitutive law, with very low elastic stiffness Kn and Kt so that they do not oppose any significant resistance to opening. For the interface elements along the fracture path, an elastoplastic constitutive model with fracture energy-based evolution laws was used

Fracture analog of the sub-Andean Devonian of southern Bolivia: Lidar applied to Abra Del Condor, in G. Zamora, K. R. McClay, and V. A. Ramos, eds., Petroleum basins and hydrocarbon potential of the Andes of Peru and Bolivia

Tight fractured sandstones of the Devonian Huamampampa Formation are associated with large gas discoveries in the sub-Andean fold-and-thrust belt of southern Bolivia. A LIDAR-based fracture characterization of the Abra del Condor backlimb anticline, a structuralstratigraphic analog, is used as the basis for a fracture stratigraphy determination. Fracture characterization using LIDAR is integrated with outcrop scanlines and is framed by stratigraphy and structural positions within this thrust-related anticline. SEFL software was used to process LIDAR data, dividing the outcrop by orientations. A workflow to extract modeled fracture planes and their associated orientations, lengths, and heights results in five fracture sets, partially validated by fracture outcrop scanlines. Multiple virtual scanlines are used to measure fracture intensity, identify fracture stratigraphic units, and define fracture-associated parameters of abundance and size distribution