Vielfache Bruchstrukturen bei einfachen Beanspruchungen — Rechnerische Untersuchungen mit Hilfe der Finite-Element-Methode

Zusammenfassung Die Finite-Element-Methode ist ein neues, auf die Benutzung von Großrechenanlagen ausgerichtetes Rechenverfahren, das sich mit Erfolg bei der Lösung von Spannungs-Verformungs-Problemen anwenden läßt. Diese Methode erwies sich bei technischen Problemen allen anderen Rechenverfahren gegenüber als überlegen und hat in den letzten Jahren in nahezu allen Sparten des Ingenieurwesens, die Felsmechanik eingeschlossen, Eingang gefunden. Seitdem einige Finite-Element-Versionen (z. B.Malina, 1969) elastoplastisches Stoffverhalten und Spannungsumlagerungen beim lokalen Bruchvorgang erfassen können, eignet sich diese Methode auch zur Simulation bestimmter tektonischer Vorgänge.Rechnerisch untersucht wird eine Scheibe der Dicke 1 (ebener Verformungszustand), in der ein kurzer Anriß vorgegeben ist. Ab einer gewissen Beanspruchung pflanzt sich der Anriß durch die Anlage von Neubrüchen fort. Es entsteht eine Scherzone, die durch ein vielfältiges Bruchflächengefüge gekennzeichnet ist. Die entstehenden Bruchflächenmuster hängen auf empfindliche Weise von den Neigungsunterschieden der im Mohrschen Diagramm definierten Grenzkurven für das unzerbrochene Material bzw. für die Diskontinuitäten ab. Die Kohäsion bleibt ohne Einfluß auf das Bruchflächengefüge.Die Rechnung entwickelt zwei Generationen von Scherflächen, ohne daß die von außen auf den Probekörper einwirkende Belastung in Größe und Richtung geändert worden wäre. Nach der Anlage von primären Scherflächenpaaren kommt es zur Entwicklung sekundärer Scherflächenpaare. Primäre und sekundäre Scherflächenpaare haben voneinander abweichende Stellungen. Sie bilden zusammen komplexe Scherkörper, die je nach Materialeigenschaften als Scherlinsen oder als Pseudo-Mohrsche Scherflächenkörper ausgebildet sind. Nach dem Entlasten bleiben im untersuchten Körper beträchtliche Restspannungen zurück, die zur Umkehr der Relativbewegungen entlang bestehender Gleitflächen und zur Anlage von neuen Zugrissen führen.Die in der Rechnung gewonnenen Ergebnisse ermöglichen ein vertieftes Verständnis vieler tektonischer und geotechnischer Erscheinungen (z. B. syn- und antithetische Brüche, Störungen zweiter Ordnung, Scherlinsen, Fiederklüfte, Restspannungen, oberflächenparallele Klüfte, verschiedene Striemungsrichtungen auf einer Harnischfläche, ac-Klüfte, Bruchstrukturen im Dreiaxialversuch). Die Rechnung zeigt, daß vielfache Bruchstrukturen nicht immer durch mehrphasige, unterschiedlich gerichtete Beanspruchungen verursacht sein müssen. Es ist denkbar, daß sie auch bei einfachen Beanspruchungen entstehen können.

---

The finite-element-method is a new calculation method, which using big computers can successfully be applied to stress-strain problems. In solving technical problems this method has proved to be better than all other present calculation methods and has been successfully applied in nearly all branches of engineering science, including rock mechanics. Since some finite-element-versions (e.g.Malina, 1969) are able to regard elastoplastic material properties and stress rearrangements due to local fracture this method is also suitable for simulating certain tectonic processes.A plate of unit thickness (plane strain) is considered, in which a small crack already exists. At a certain state of stress, the crack propagates by creating new rupture surfaces. A shear zone develops, which is characterized by a complex fracture fabric. The developing patterns of shear and sliding surfaces change sensitively with the inclination of theMohr envelopes for the unfractured material and the discontinuities resp. Cohesion has no influence on the fracture fabric.The computation gives two generations of shear surfaces, without any change in the magnitude and direction of the external load. After the formation of primary pairs of shear surfaces secondary pairs develop. Primary and secondary shear surfaces differ in their orientation. Together they form complex shear bodies, which, depending on material properties, are shear lenses or quasi-Mohr-shear bodies. With the release of the external load, residual stresses are stored within the plate. These stresses cause a reversal in the sense of sliding movement and development of new fractures.The results of computation allow a better understanding of many tectonic and geotechnical phenomena (e. g. synthetic and antithetic fractures, second order faults, shear lenses, feather joints, residual stresses, surface-parallel joints, different striae orientations on slickensided surfaces, ac-joints, fracture patterns in triaxial compression tests). The computations show that multiple fracture patterns are not necessarily caused by poly-phase, multiaxial loading. It seems possible that these patterns can also develop under simple load conditions.

---

Résumé La méthode des éléments finis permet de résoudre avec beaucoup de succès les problèmes de contraintes et de déformations à l'aide d'ordinateurs de grande capacité. Cette méthode s'est révélée supérieure à tous les autres procédés de calcul pour résoudre les problèmes techniques et elle est appliquée depuis quelques années à presque toutes les branches des sciences de l'Ingénieur, y compris la mécanique des roches. Depuis que quelques versions de calcul des éléments finis (par exempleMalina, 1969) peuvent prendre en considération le comportement élasto-plastique des matériaux et les redistributions des contraintes dues à des ruptures locales, cette méthode convient aussi à la simulation des phénomènes tectoniques.Un disque d'épaisseur unitaire (déformations planes), dans lequel existe une courte fissure est soumis au calcul. A partir d'une certaine sollicitation, la fissure se prolonge par de nouvelles ruptures. Ainsi se forme une zone de cisaillement qui est caractérisée par une structure multiple des surfaces de rupture. Le dessin ainsi formé par ces surfaces de rupture dépend sensiblement des différences d'inclinaison des courbes limites telles qu'elles sont définies dans le diagramme deMohr, soit pour le matériau intact, soit pour les discontinuités. La »cohésion« n'a aucune influence sur la structure de ces surfaces de rupture.Le calcul développe deux générations de surfaces de rupture, sans que l'intensité et la direction de la charge extérieure affectant l'échantillon ne soient modifiées. Après le développement des couples de surfaces de rupture primaires, se forment des couples de surfaces de rupture secondaires dont les directions sont différentes des premières. Elles forment ensemble des corps de cisaillement complexe, qui, selon les propriétés du matériau, s'expriment sous la forme soit de lentilles de cisaillement, soit de corps » pseudo-Mohr « de surfaces de cisaillement. Après suppression de la charge, il persiste dans le modèle des contraintes résiduelles importantes qui entraînent des mouvements de retrait le long des surfaces de glissement existantes et de nouvelles fissures dues aux efforts de traction.Les résultats obtenus par le calcul permettent de mieux comprendre de nombreux phénomènes tectoniques et géotechniques (tels que les fractures synthétiques et antithétiques, les failles de deuxième ordre, les lentilles de cisaillement, les fissures réticulées, les contraintes résiduelles, les fissures parallèles à la surface, les diverses directions des stries de glissement sur les miroirs de faille, les fissures de type ac, les structures de fracture dan l'essai triaxial). Le calcul montre que les multiples structures de fracture ne sont pas toujours dues à des efforts multiphasés de directions différentes. On peut penser qu'elles se forment également sous des états de contraintes simples.

---

, Finite-Element, , - . . , (Malina, 1969), . 1 ( ), ; , . , . , Mohr'a , . .

---

---

Diese Arbeit ist im SFB 77 - Felsmechanik« -, Karlsruhe, entstanden. Der Autor dankt für vielfache Hilfe und Anregungen, insbesondere den Herren Dr. H.Malina und cand. geol. B.Hoffers.