Cette thèse porte sur l'étude expérimentale de la photosensibilité de verres de chalcogénures (ChG) sous la forme de couches minces. Plus particulièrement, elle est dédiée à l’étude des modifications photoinduites de leurs propriétés optiques ainsi qu’aux changements structuraux qui y sont liés au niveau atomique. Une étude systématique des propriétés des ChG sélectionnés dans le système vitreux Ge-As-S a été réalisée en fonction de la concentration relative des éléments Ge, As et S, de l’épaisseur des couches minces déposées ainsi que des différentes conditions expérimentales d’exposition au faisceau laser. Tout d’abord nous nous sommes intéressés au band gap optique du matériau, au décalage du bord d'absorption et au changement de sa pente qui sont les résultats d’arrangements atomiques complexes dans le réseau désordonné du ChG. Ensuite, les résultats expérimentaux ont démontré que la composition vitreuse Ge25As30S45 possède la plus forte photosensibilité et notamment la valeur la plus élevée de biréfringence photo-induite (PIB) parmi les verres des systèmes Ge-As-S et As-S. La conversion de liaisons homopolaires (Ge-Ge, As-As) à hétéropolaires (Ge-S, As-S) a de plus été mise en évidence pour expliquer ce phénomène. En outre, la modélisation théorique simple que nous avons proposée avec une certaine approximation, montre que la valeur locale du PIB peut être d’un ordre de grandeur plus élevée que sa valeur moyenne. Les changements dynamiques d’absorption photo-induite étudiés pour différentes conditions expérimentales sont caractérisés par de forts changements asymétriques et non-monotones durant l'excitation et la relaxation. Ces changements ont été décrits par un modèle phénoménologique unipolaire que nous avons proposé et qui est basé sur certaines conversions séquentielles de liaisons se produisant après le franchissement d’une barrière énergétique donnée (estimée sur la base de nos mesures). Puis cette photosensibilité élevée des couches minces Ge-As-S a été utilisée pour l'enregistrement d’un réseau polarisé et pour la fabrication d’une lentille à gradient d’indice (GRIN) sur la surface, obtenus par irradiation laser à une longueur d’onde correspondante à la valeur de son band gap optique. La variation des efficacités de diffraction maximale obtenues pour les hologrammes scalaires et vectoriels a été discutée en considérant les différentes unités structurales identifiées et le rôle des transitions électroniques directes et indirectes dans ces deux types de réseaux. La stabilité thermique des hologrammes vectoriels a été montrée expérimentalement grâce à l’ajout de l’élément germanium Ge dans la composition de la couche mince. Enfin, les forces optiques des lentilles obtenues ainsi que les distorsions de front d'onde et l’effet de vieillissement ont été caractérisés à l’aide de capteurs Shack Hartmann. / This PhD thesis refers to the experimental study of photosensitivity of chalcogenide glassy (ChG) thin films and their induced structural changes at the atomic level. A systematic study of the ChG properties is presented as a function of the elemental composition in the selected Ge-As-S system and the film thickness. More particularly, the goals of this work were to evidence and characterize the photoinduced birefringence and dichroism effects, to investigate the mechanisms involved and to correlate experimental observations with theoretical modeling. The first part of the work was dedicated to the study of the optical properties, specifically the optical band gap of the prepared composition within the Ge-As-S vitreous system to reveal the most appropriate composition for further photoinduced effects examination. The shift and slope change observed for the absorption edge (associated with the optical band gap) according to the film thickness resulted from complex atomic (re)arrangements in the ChG network. The experiments carried out for the photoinduced effects have permitted to determine the best composition to be Ge25As30S45 among the Ge-As-S and As-S glasses in terms of higher photosensitivity and higher value of photoinduced birefringence (PIB) produced by the conversion from homopolar (Ge-Ge, As-As) to heteropolar (Ge-S, As-S) bonds. Moreover, the simple theoretical model proposed herein showed, with some approximation, that the local value of the PIB in these ChG thin films may be one order of magnitude higher than its average value. Then, the dynamic study of the photoinduced absorption revealed a strong asymmetric and non-monotonic behavior as a function of the irradiation laser power. To account for this specific behavior, a new unipolar phenomenological model is proposed based on sequential bond conversions occurring beyond an estimated energetic barrier. The photoinduced anisotropy of these ChG Ge-As-S thin films was then used to record polarization gratings and gradient index lenses (GRIN). The maximum diffraction efficiencies achieved between scalar and vector holograms was discussed by means of involved structural units and the role played by indirect and direct electronic transitions. In addition, an improved thermal stability of the recorded vector holograms was experimentally shown after incorporation of germanium Ge into the material composition. The optical performance of the obtained lenses as well as the wave front distortions, aging effect and so on were studied by means of Shack Hartmann sensor.

Cette thèse dresse un cadre d'analyse qui permet de répondre aux questions suivantes: quelle est l'évolution de la projection longitudinale des faisceaux laser vectoriels dans les milieux anisotropes, et quelle est son influence dans les processus non linéaires d'interaction photonique? La projection vectorielle longitudinale est généralement négligée lors de la propagation dans les milieux anisotropes et en optique paraxiale. Or elle devient importante hors de ce régime, habituellement en focalisant à grand angle les faisceaux laser. On décrit le champ électromagnétique exact d'un faisceau laser non paraxial Laguerre-Gauss élégant TM₀₁ à pincement critique. Les conventions anisotropes adoptées généralisent le classement des milieux biaxiaux selon leur état normalisé. L'orientation énergétiquement orthogonale du repère propagatif propre permet la description de l'évolution des champs vectoriels selon les modes lents et rapides transverses et longitudinaux. L'imperméabilité tensorielle dans ce repère propre procure tous les paramètres propagatifs essentiels dont le gradient diamétral. Le modèle non linéaire développé est une extension tensorielle non paraxiale exacte des notions scalaires et paraxiales disponibles dans la littérature. On développe l'outillage nécessaire pour les représentations complexes par phaseurs temporels des projections des champs lors du mélange ondulatoire vectoriel. On effectue la représentation par phaseurs confluents dans un repère propagatif propre d'ondes planes, ce qui permet de calculer l'évolution de chacune des fréquences angulaires à l'intérieur du milieu non linéaire anisotrope. Cette approche opérationnelle est nécessaire pour prévoir correctement les échanges photoniques lors des processus non linéaires non paraxiaux qui impliquent la projection longitudinale. On applique ce modèle à la solution de l'équation d'onde vectorielle anisotrope non linéaire dans le cas de la génération de seconde harmonique sous approximation de pompe non dépeuplée. De nouveaux processus non linéaires sont décrits, faisant intervenir les modes propagatifs longitudinaux. Chaque processus est caractérisé par une fonction de transfert modale qui reflète l'efficacité de pompage et la cohérence de l'interaction. Le faisceau TM₀₁ pourra avantageusement être utilisé comme pompe dans les processus non linéaires modélisés. Le modèle peut également s'appliquer à toute autre onde de pompage.

In der vorliegenden Arbeit wird ein gemischtes Element gegen Locking-Effekte untersucht. Dazu wird ein Fünf-Feld-Hu-Washizu-Funktional (CoFEM-Element) für lineare und quadratische Hexaeder-Elemente unter einer hyperelastischen, isotropen, polykonvexen sowie einer transversal-isotropen Materialformulierung implementiert. Die resultierenden nichtlinearen Gleichungen werden mithilfe eines Mehrebenen-NEWTON-RAPHSON-Verfahren unter Beachtung einer konsistenten Linearisierung gelöst. Als repräsentatives Beispiel der numerischen Untersuchungen dient der einseitig eingespannte Cook-Balken mit einer quadratischen Druckverteilung am Rand. Zur Beurteilung des CoFEM-Elements wird das räumliche Konvergenzverhalten für unterschiedliche Polynomgrade und für verschiedene Netze unter Beachtung der algorithmischen Effizienz untersucht. / This paper presents a mixed finite element formulation of Hu-Washizu type (CoFEM) designed to reduce locking effects with respect to a linear and quadratic approximation in space. We consider a hyperelastic, isotropic, polyconvex material formulation as well as transverse isotropy. The resulting nonlinear algebraic equations are solved with a multilevel NEWTON-RAPHSON method. As a numerical example serves a cook-like cantilever beam with a quadratic distribution of in-plane load on the Neumann boundary. We analyze the spatial convergence with respect to the polynomial degree of the underlying Lagrange polynomials and with respect to the level of mesh refinement in terms of algorithmic efficiency.

The variation of rock strength and stiffness, known as mechanical anisotropy, is expected at different scales: large (rock mass) - or small (intact rock) - scales. It is always mandatory for engineering applications built either on or in anisotropic rock masses to investigate the strength and deformation behavior of those masses. To achieve this goal, continuum-based constitutive models are presented to analyze the mechanical anisotropy. One of both implemented models is named ‘Transubi model’ which considers the transverse isotropic elasticity into bi-linear Mohr-Coulomb strain hardening/softening plastic framework. Experimental investigations and numerical simulations focused mainly on the influence of the mechanical anisotropy on the plastic zoning around excavated openings in laminated rocks. Later, the Transubi model was applied to a tunnel excavated in a shaly facies formation of bedded argillaceous Opalinus clay in an URL (FE-tunnel) to obtain the short-term stability insights. Overall, the research outcomes may have a prospective impact regarding the understanding of anisotropy of laminated, bedded and foliated rocks which improves the deformation behaviour predictability using continuum-based numerical modeling tools.

The present work describes investigations on the anisotropic strength behavior of rocks in the splitting tensile test (Brazilian test). Three transversely isotropic rocks (gneiss, slate and sandstone) were studied in the Lab. A total of more than 550 indirect tensile strength tests were conducted, with emphasis was placed on the investigation of the influence of the spatial position of anisotropic weakness plane to the direction of the load on the fracture strength and fracture or fracture mode. In parallel, analytical solutions were evaluated for stress distribution and developed 3D numerical models to study the stress distribution and the fracture mode at the transversely isotropic disc. There were new findings on the fracture mode of crack propagation, the influence of the disc thickness, the influence of the applying loading angle and angle of the loading-foliation for transversely isotropic material. / Inhalt der Arbeit sind Untersuchungen zum anisotropen Festigkeitsverhalten von Gesteinen beim Spaltzugversuch (Brazilian Test). Laborativ wurden drei transversalisotrope Gesteine (Granit, Schiefer und Sandstein) untersucht. Insgesamt wurden mehr als 550 Spaltzugversuche durchgeführt, wobei der Schwerpunkt auf die Untersuchung des Einflusses der räumlichen Lage der Anisotropieebene zur Richtung des Lasteintrages auf die Bruchfestigkeit und das Bruchbild bzw. den Bruchmodus gelegt wurde. Parallel dazu wurden analytische Lösungen zur Spannungsverteilung ausgewertet sowie numerische 3D-Modelle entwickelt, um die Spannungsverteilung sowie den Bruchmodus bei einer transversalisotropen Scheibe zu untersuchen. Es wurden neue Erkenntnisse zum Bruchmodus, der Rissausbreitung, des Einflusses der Scheibendicke, dem Einfluss des Lasteinleitungswinkel sowie des Winkels Lasteintrag - Anisotropieebene für transversalisotropes Material gewonnen.

This thesis covers different problems that arise due to crystal and pseudospin anisotropy present in 3d and 5d transition metal oxides. We demonstrate that the methods of computational quantum chemistry can be fruitfully used for quantitative studies of such problems. In Chapter 2, Chapter 3, and Chapter 7 we show that it is possible to reliably calculate local multiplet splittings fully ab initio, and therefore help to assign peaks in experimental spectra to corresponding electronic states. In a situation of large number of peaks due to low local symmetry such assignment using semi-empirical methods can be very tedious and non-unique. Moreover, in Chapter 4 we present a computational scheme for calculating intensities as observed in the resonant inelastic X-ray scattering and X-ray absorption experiments. In our scheme highly-excited core-hole states are calculated explicitly taking into account corresponding orbital relaxation and electron polarization. Computed Cu L-edge spectra for the Li2CuO2 compound reproduce all features present in experiment. Unbiased ab initio calculations allow us to unravel a delicate interplay between the distortion of the local ligand cage around the transition metal ions and the anisotropic electrostatic interactions due to second and farther coordination shells. As shown in Chapter 5 and Chapter 6 this interplay can lead to the counter intuitive multiplet structure, single-ion anisotropy, and magnetic g factors. The effect is quite general and may occur in compounds with large difference between charges of metal ions that form anisotropic environment around the transition metal, like Ir 4+ in plane versus Sr 2+ out of plane in the case of Sr2IrO4. An important aspect of the presented study is the mapping of the quantum chemistry results onto simpler physical models, namely extended Heisenberg model, providing an ab initio parametrization. In Chapter 5 we employ the effective Hamiltonian technique for extracting parameters of the anisotropic Heisenberg model with single-ion anisotropy in the case of quenched orbital moment and second-order spin-orbit coupling. Calculated strong easy-axis anisotropy of the same order of magnitude as the symmetric exchange is consistent with experimentally-observer all-in/all-out magnetic order. In Chapter 6 we introduce new flavour of the mapping procedure applicable to systems with first-order spin-orbit coupling, such as 5d 5 iridates based on analysis of the wavefunction and interaction with magnetic field. In Chapter 6 and Chapter 7 we use this new procedure to obtain parameters of the pseudospin anisotropic Heisenberg model. We find large antisymmetric exchange leading to the canted antiferromagnetic state in Sr2IrO4 and nearly ideal one-dimensional Heisenberg behaviour of the CaIrO3, both agree very well with experimental findings.

The acoustics of small cavities raises interest of the scientific community since it involves particular damping mechanisms. In fluid dynamics, when a small perturbation is propagating within a Newtonian and heat-conducting fluid bounded by a rigid and isothermal surface, viscous and thermal dissipative mechanisms are generated near the walls. Such effects can have significant impact on the acoustic behaviour of the system.<p>Several types of practical applications can be cited, among which: hearing aids, micro-electro-mechanical systems (transducers, microphones and loud-speakers), absorbing materials made of thin capillary nets or small pores, dissipative silencers, thermo-acoustic heat exchangers, or any kind of device bringing into play small resonant cavities filled with a dissipative fluid (micro-acoustics).<p><p>This study focuses on appropriated reductions of the physical equations, in order to enhance the efficiency of the numerical resolution without adversely affecting the accuracy. Moreover, the proposed strategies lead to numerically stable systems as they involve only one scalar partial order differential equation (or equivalent fluid equation). The emphasis is put on the physical aspect of those reductions, their range of applicability, benefits and drawbacks.<p>Two new reduced models are proposed to estimate the visco-thermal acoustic wave propagation. A first extension deals with waveguide geometries and relax the hypothesis of the fluid at rest. The second original formulation addresses visco-thermal acoustics in 3D arbitrary geometries. This model is based on different considerations coming from existing techniques as well as the estimation of a wall-distance field.<p><p>A second part aims at studying the acoustic behaviour of biphasic materials and more specifically poro-elastic materials. This type of acoustic component is widely used in industry because of their good absorbing properties in the medium- and high-frequency <p>ranges.<p>A preliminary bibliographic research deals with the derivation of the set of partial order differential equations that account for both fluid/structure interactions and the anisotropy of a given poro-elastic material. It has been shown that transversely orientated capillary materials (for instance catalyst substrates) can be simulated using the proposed reduction technique.<p>At last, the modelling of the acoustic transmission between two domains separated by perforated or micro-perforated plates or thin plates of poro-elastic materials is discussed. The analogy between the rigid perforated plate models with an equivalent fluid formulation has been presented. As a result, this model has been extended in order to account for flexural effects of the solid part.<p><p><p>Ce travail porte sur l'étude de certains phénomènes d'amortissements intervenant dans l'acoustique des petites cavités. En méchanique des fluides, lorsqu'une petite perturbation se propage au sein d'un fluide newtonien et caloporteur borné par un mur rigide et isotherme, ces mécanismes dissipatifs particuliers se localisent aux abords des parois et jouent un rôle significatif dans certaines situations.<p>Parmi les exemples d'applications pratiques, il est possible de citer les appareils d'aide auditive, les systèmes microélectromécaniques (transducteurs, microphones et haut-parleurs), les matériaux absorbants constitués de fins réseaux capillaires ou de pores aux dimensions réduites, les systèmes de silencieux, d'échangeurs de chaleur thermo-acoustiques ou tout autre appareil mettant en jeu des cavités résonantes aux dimensions réduites (micro-acoustique).<p><p>L'étude proposée ici se focalise sur des stratégies de réduction appropriées des équations physiques, ceci afin d'améliorer l'efficacité du modèle tout en conservant une précision acceptable. Les techniques présentées aboutissent à des systèmes numériquement stables mettant en jeu une seule équation scalaire (ou équation fluide équivalent). Ainsi, l'accent est porté sur l'aspect physique des réductions, leurs domaines d'application, avantages et inconvénients.<p>Deux modèles originaux sont proposés afin de prédire la propagation acoustique visco-thermique. Une première extension permet d'évaluer la pression acoustique au sein de géométries particulières de type guides d'onde en présence d'un écoulement hydrodynamique. La seconde formulation présentée s'intéresse à l'acoustique dans des domaines 3D arbitraires. Cette méthode se base sur des considérations conjointes de modèles réduits existants ainsi que sur l'estimation d'un champ de distance à la plus proche paroi.<p><p>Dans une seconde partie, nous nous proposons d'étudier le comportement acoustique de matériaux biphasique et plus précisément les matériaux poro-élastiques (très utilisés dans l'industrie en raison de leurs caractéristiques absorbantes dans les domaines des moyennes et hautes fréquences).<p>Une étude bibliographique préliminaire nous a permis d'exprimer l'ensemble des équations aux dérivées <p>partielles modélisant à la fois les interactions fluide/structure et l'anisotropie générale des matériaux. <p>Cette réflexion nous a permis d'aboutir à un modèle de matériau isotrope transverse intéressant, combinant le modèle fluide proposé et la formulation acousto-élastique équivalente. Ainsi la modélisation de structures capillaires orientées (comme les matériaux utilisés dans les catalyseurs automobiles) s'en trouve grandement simplifiée.<p>Enfin, la transmission acoustique intervenant entre deux domaines fluides séparés par une plaque perforée ou micro-perforée ou bien une couche de matériau poreux a été étudiée. L'analogie entre les modèlisations existantes et un modèle générique de fluide équivalent a été mise en évidence. Pour finir, cette formulation a été étendue afin de prendre en compte les effets de flexion de la partie solide. / Doctorat en Sciences de l'ingénieur / info:eu-repo/semantics/nonPublished

This thesis covers different problems that arise due to crystal and pseudospin anisotropy present in 3d and 5d transition metal oxides. We demonstrate that the methods of computational quantum chemistry can be fruitfully used for quantitative studies of such problems. In Chapter 2, Chapter 3, and Chapter 7 we show that it is possible to reliably calculate local multiplet splittings fully ab initio, and therefore help to assign peaks in experimental spectra to corresponding electronic states. In a situation of large number of peaks due to low local symmetry such assignment using semi-empirical methods can be very tedious and non-unique. Moreover, in Chapter 4 we present a computational scheme for calculating intensities as observed in the resonant inelastic X-ray scattering and X-ray absorption experiments. In our scheme highly-excited core-hole states are calculated explicitly taking into account corresponding orbital relaxation and electron polarization. Computed Cu L-edge spectra for the Li2CuO2 compound reproduce all features present in experiment. Unbiased ab initio calculations allow us to unravel a delicate interplay between the distortion of the local ligand cage around the transition metal ions and the anisotropic electrostatic interactions due to second and farther coordination shells. As shown in Chapter 5 and Chapter 6 this interplay can lead to the counter intuitive multiplet structure, single-ion anisotropy, and magnetic g factors. The effect is quite general and may occur in compounds with large difference between charges of metal ions that form anisotropic environment around the transition metal, like Ir 4+ in plane versus Sr 2+ out of plane in the case of Sr2IrO4. An important aspect of the presented study is the mapping of the quantum chemistry results onto simpler physical models, namely extended Heisenberg model, providing an ab initio parametrization. In Chapter 5 we employ the effective Hamiltonian technique for extracting parameters of the anisotropic Heisenberg model with single-ion anisotropy in the case of quenched orbital moment and second-order spin-orbit coupling. Calculated strong easy-axis anisotropy of the same order of magnitude as the symmetric exchange is consistent with experimentally-observer all-in/all-out magnetic order. In Chapter 6 we introduce new flavour of the mapping procedure applicable to systems with first-order spin-orbit coupling, such as 5d 5 iridates based on analysis of the wavefunction and interaction with magnetic field. In Chapter 6 and Chapter 7 we use this new procedure to obtain parameters of the pseudospin anisotropic Heisenberg model. We find large antisymmetric exchange leading to the canted antiferromagnetic state in Sr2IrO4 and nearly ideal one-dimensional Heisenberg behaviour of the CaIrO3, both agree very well with experimental findings.

In der vorliegenden Arbeit wird ein gemischtes Element gegen Locking-Effekte untersucht. Dazu wird ein Fünf-Feld-Hu-Washizu-Funktional (CoFEM-Element) für lineare und quadratische Hexaeder-Elemente unter einer hyperelastischen, isotropen, polykonvexen sowie einer transversal-isotropen Materialformulierung implementiert. Die resultierenden nichtlinearen Gleichungen werden mithilfe eines Mehrebenen-NEWTON-RAPHSON-Verfahren unter Beachtung einer konsistenten Linearisierung gelöst. Als repräsentatives Beispiel der numerischen Untersuchungen dient der einseitig eingespannte Cook-Balken mit einer quadratischen Druckverteilung am Rand. Zur Beurteilung des CoFEM-Elements wird das räumliche Konvergenzverhalten für unterschiedliche Polynomgrade und für verschiedene Netze unter Beachtung der algorithmischen Effizienz untersucht. / This paper presents a mixed finite element formulation of Hu-Washizu type (CoFEM) designed to reduce locking effects with respect to a linear and quadratic approximation in space. We consider a hyperelastic, isotropic, polyconvex material formulation as well as transverse isotropy. The resulting nonlinear algebraic equations are solved with a multilevel NEWTON-RAPHSON method. As a numerical example serves a cook-like cantilever beam with a quadratic distribution of in-plane load on the Neumann boundary. We analyze the spatial convergence with respect to the polynomial degree of the underlying Lagrange polynomials and with respect to the level of mesh refinement in terms of algorithmic efficiency.

The present work describes investigations on the anisotropic strength behavior of rocks in the splitting tensile test (Brazilian test). Three transversely isotropic rocks (gneiss, slate and sandstone) were studied in the Lab. A total of more than 550 indirect tensile strength tests were conducted, with emphasis was placed on the investigation of the influence of the spatial position of anisotropic weakness plane to the direction of the load on the fracture strength and fracture or fracture mode. In parallel, analytical solutions were evaluated for stress distribution and developed 3D numerical models to study the stress distribution and the fracture mode at the transversely isotropic disc. There were new findings on the fracture mode of crack propagation, the influence of the disc thickness, the influence of the applying loading angle and angle of the loading-foliation for transversely isotropic material.:ACKNOWLEDGMENTS 5 ABSTRACT 7 TABLE OF CONTENTS 9 LIST OF FIGURES 13 LIST OF TABLES 19 I. INTRODUCTION 21 Objective of this work 22 Scope of work 23 Research procedure 23 Significance of the work 24 Layout 24 1 STATE OF THE ART 27 1.1 Review of the Brazilian tensile strength test 27 1.1.1 General overview 27 1.1.2 Development of the Brazilian tensile strength test 29 1.1.3 The Brazilian tensile strength test on anisotropic rocks 31 1.1.4 Summary 32 1.2 Analytical aspects 33 1.2.1 Hypotheses for the conventional Brazilian test 34 1.2.2 Failure criteria 36 1.2.3 Crack initiation and propagation 39 1.2.4 Summary 41 1.3 Numerical considerations 41 1.3.1 Numerical methods 42 1.3.2 Summary 42 1.4 Conclusion 43 2 DIAMETRAL COMPRESSION IN A SOLID DISC – COMPILATION OF ANALYTICAL AND SEMI-ANALYTICAL SOLUTIONS 45 2.1 Introduction 45 2.2 Diametral compressive stress distribution in an isotropic elastic disc 45 2.2.1 Elastic theory of line load 46 2.2.2 2D analytical solutions 47 2.2.3 3D disc under line and diametral compressive distributed loads 55 2.2.4 3D solution under diametral compressive distributed load 56 2.3 Stress and strain in an isotropic solid disc 59 2.4 Stress and strain in anisotropic rocks 61 2.5 Conclusion 65 3 LABORATORY TESTS 69 3.1 Introduction 69 3.2 Laboratory test program 70 3.3 Sample preparation 71 3.4 Ultrasonic measurements 72 3.5 Uniaxial and triaxial compression tests 73 3.5.1 Uniaxial compression test 73 3.5.2 Triaxial compression tests 74 3.6 Brazilian tensile strength tests 76 3.6.1 Test apparatus 76 3.6.2 Laboratory test results 77 3.6.3 Interpretation of the test results 89 3.7 Conclusion 96 4 NUMERICAL SIMULATION OF ISOTROPIC MATERIALS - COMPARISON WITH ANALYTICAL SOLUTIONS 97 4.1 Introduction 97 4.2 Numerical simulation of isotropic materials 97 4.2.1 FLAC3D simulation program 97 4.2.2 Simulation procedure 98 4.2.3 Numerical model setup 98 4.2.4 Influence of mesh type 99 4.2.5 Influence of specimen thickness 100 4.2.6 Influence of Poisson’s ratio 102 4.2.7 Influence of loading angle (2) 106 4.2.8 Comparison of 3D analytical and numerical results 110 4.2.9 Influence of stress concentration at the loading jaws 112 4.3 Comparison with experimental results of Postaer Sandstone (FG.Ss) 112 4.4 Conclusion 114 5 NUMERICAL SIMULATION OF ANISOTROPIC MATERIALS - COMPARISON WITH LABORATORY TESTS 117 5.1 Introduction 117 5.2 General procedure for simulating the Brazilian test using FLAC3D 117 5.2.1 Conceptual model 119 5.2.2 Boundary Conditions 119 5.2.3 Numerical model set-up 120 5.3 Constitutive model 121 5.3.1 Choice of constitutive model 121 5.3.2 Bilinear Strain-Hardening/Softening Ubiquitous-Joint Model [98] 121 5.4 Parameter calibration 124 5.4.1 Material parameters used 124 5.4.2 Contact between disc and loading jaws 126 5.4.3 Post-failure deformation properties 128 5.4.4 Tension cut-off 129 5.5 Numerical simulation results 131 5.5.1 Introduction 131 5.5.2 Stress distribution and failure state 133 5.5.3 Stress state in an isotropic elastic medium with arbitrary orientation planes 136 5.5.4 Plasticity states 139 5.5.5 Damage and fracture process 141 5.5.6 Fracture patterns – Comparison of lab results and numerical simulations 148 5.6 Tensile strength – Comparison of lab results and numerical simulations 149 5.6.1 Tensile strength of Le.Gs Gneiss 150 5.6.2 Tensile strength of My.Sc Slate 155 5.7 Summary and Review 159 5.7.1 Potential failure state deduced from pure elastic considerations 159 5.7.2 Tensile strength distribution 160 5.7.3 Tensile strength – determining the anisotropy factor 161 5.7.4 Tensile strength – different procedures - different results 163 6 CONCLUSION AND RECOMMENDATIONS 165 APPENDICES 171 Appendix 3.1 - Fracture patterns in FG.Ss samples 171 Appendix 3.2 - Fracture patterns in FG.Gs samples 177 Appendix 3.3 - Fracture patterns in Le.Gs samples 183 Appendix 3.4 - Fracture patterns in My.Sc samples 190 Appendix 4.1 - Influence of loading angle 197 Appendix 4.2 - Influence of material properties 203 Appendix 5.1 - Failure zone state in Le.Gs Gneiss 209 Appendix 5.2: Failure zone state in My.Sc Slate 216 REFERENCES 223 / Inhalt der Arbeit sind Untersuchungen zum anisotropen Festigkeitsverhalten von Gesteinen beim Spaltzugversuch (Brazilian Test). Laborativ wurden drei transversalisotrope Gesteine (Granit, Schiefer und Sandstein) untersucht. Insgesamt wurden mehr als 550 Spaltzugversuche durchgeführt, wobei der Schwerpunkt auf die Untersuchung des Einflusses der räumlichen Lage der Anisotropieebene zur Richtung des Lasteintrages auf die Bruchfestigkeit und das Bruchbild bzw. den Bruchmodus gelegt wurde. Parallel dazu wurden analytische Lösungen zur Spannungsverteilung ausgewertet sowie numerische 3D-Modelle entwickelt, um die Spannungsverteilung sowie den Bruchmodus bei einer transversalisotropen Scheibe zu untersuchen. Es wurden neue Erkenntnisse zum Bruchmodus, der Rissausbreitung, des Einflusses der Scheibendicke, dem Einfluss des Lasteinleitungswinkel sowie des Winkels Lasteintrag - Anisotropieebene für transversalisotropes Material gewonnen.:ACKNOWLEDGMENTS 5 ABSTRACT 7 TABLE OF CONTENTS 9 LIST OF FIGURES 13 LIST OF TABLES 19 I. INTRODUCTION 21 Objective of this work 22 Scope of work 23 Research procedure 23 Significance of the work 24 Layout 24 1 STATE OF THE ART 27 1.1 Review of the Brazilian tensile strength test 27 1.1.1 General overview 27 1.1.2 Development of the Brazilian tensile strength test 29 1.1.3 The Brazilian tensile strength test on anisotropic rocks 31 1.1.4 Summary 32 1.2 Analytical aspects 33 1.2.1 Hypotheses for the conventional Brazilian test 34 1.2.2 Failure criteria 36 1.2.3 Crack initiation and propagation 39 1.2.4 Summary 41 1.3 Numerical considerations 41 1.3.1 Numerical methods 42 1.3.2 Summary 42 1.4 Conclusion 43 2 DIAMETRAL COMPRESSION IN A SOLID DISC – COMPILATION OF ANALYTICAL AND SEMI-ANALYTICAL SOLUTIONS 45 2.1 Introduction 45 2.2 Diametral compressive stress distribution in an isotropic elastic disc 45 2.2.1 Elastic theory of line load 46 2.2.2 2D analytical solutions 47 2.2.3 3D disc under line and diametral compressive distributed loads 55 2.2.4 3D solution under diametral compressive distributed load 56 2.3 Stress and strain in an isotropic solid disc 59 2.4 Stress and strain in anisotropic rocks 61 2.5 Conclusion 65 3 LABORATORY TESTS 69 3.1 Introduction 69 3.2 Laboratory test program 70 3.3 Sample preparation 71 3.4 Ultrasonic measurements 72 3.5 Uniaxial and triaxial compression tests 73 3.5.1 Uniaxial compression test 73 3.5.2 Triaxial compression tests 74 3.6 Brazilian tensile strength tests 76 3.6.1 Test apparatus 76 3.6.2 Laboratory test results 77 3.6.3 Interpretation of the test results 89 3.7 Conclusion 96 4 NUMERICAL SIMULATION OF ISOTROPIC MATERIALS - COMPARISON WITH ANALYTICAL SOLUTIONS 97 4.1 Introduction 97 4.2 Numerical simulation of isotropic materials 97 4.2.1 FLAC3D simulation program 97 4.2.2 Simulation procedure 98 4.2.3 Numerical model setup 98 4.2.4 Influence of mesh type 99 4.2.5 Influence of specimen thickness 100 4.2.6 Influence of Poisson’s ratio 102 4.2.7 Influence of loading angle (2) 106 4.2.8 Comparison of 3D analytical and numerical results 110 4.2.9 Influence of stress concentration at the loading jaws 112 4.3 Comparison with experimental results of Postaer Sandstone (FG.Ss) 112 4.4 Conclusion 114 5 NUMERICAL SIMULATION OF ANISOTROPIC MATERIALS - COMPARISON WITH LABORATORY TESTS 117 5.1 Introduction 117 5.2 General procedure for simulating the Brazilian test using FLAC3D 117 5.2.1 Conceptual model 119 5.2.2 Boundary Conditions 119 5.2.3 Numerical model set-up 120 5.3 Constitutive model 121 5.3.1 Choice of constitutive model 121 5.3.2 Bilinear Strain-Hardening/Softening Ubiquitous-Joint Model [98] 121 5.4 Parameter calibration 124 5.4.1 Material parameters used 124 5.4.2 Contact between disc and loading jaws 126 5.4.3 Post-failure deformation properties 128 5.4.4 Tension cut-off 129 5.5 Numerical simulation results 131 5.5.1 Introduction 131 5.5.2 Stress distribution and failure state 133 5.5.3 Stress state in an isotropic elastic medium with arbitrary orientation planes 136 5.5.4 Plasticity states 139 5.5.5 Damage and fracture process 141 5.5.6 Fracture patterns – Comparison of lab results and numerical simulations 148 5.6 Tensile strength – Comparison of lab results and numerical simulations 149 5.6.1 Tensile strength of Le.Gs Gneiss 150 5.6.2 Tensile strength of My.Sc Slate 155 5.7 Summary and Review 159 5.7.1 Potential failure state deduced from pure elastic considerations 159 5.7.2 Tensile strength distribution 160 5.7.3 Tensile strength – determining the anisotropy factor 161 5.7.4 Tensile strength – different procedures - different results 163 6 CONCLUSION AND RECOMMENDATIONS 165 APPENDICES 171 Appendix 3.1 - Fracture patterns in FG.Ss samples 171 Appendix 3.2 - Fracture patterns in FG.Gs samples 177 Appendix 3.3 - Fracture patterns in Le.Gs samples 183 Appendix 3.4 - Fracture patterns in My.Sc samples 190 Appendix 4.1 - Influence of loading angle 197 Appendix 4.2 - Influence of material properties 203 Appendix 5.1 - Failure zone state in Le.Gs Gneiss 209 Appendix 5.2: Failure zone state in My.Sc Slate 216 REFERENCES 223