Die Verwendung von lang kragenden Schaftfräsern im Bereich der Hochgeschwindigkeitsbearbeitung birgt besondere Herausforderungen bezüglich der Prozessdynamik. In diesem Zusammenhang werden verschiedene kontinuumsmechanische Berechnungsmodelle für Werkzeuge mit Hohlschaft vorgestellt. Dabei wird eine teilweise Füllung des Schaftes mit einer fließfähigen Ausgleichsmasse zum Zweck des automatischen Wuchtens berücksichtigt. Ausgehend von der Verformungskinematik wird die systembeschreibende Variationsformulierung mit Hilfe des Hamilton'schen Prinzips hergeleitet. Dabei wird auch auf den Einfluss von stochastisch verteilten Unwuchten, geometrischen Nichtlinearitäten und Schubdeformationen eingegangen. Zur Ortsdiskretisierung werden sowohl lokale als auch globale Methoden angewendet und miteinander verglichen. Die Auswertung stellt den Einfluss von verschiedenen geometrischen sowie prozessbedingten Parametern auf die Eigenfrequenzen, stationäre Deformation, Stabilität sowie Zeitlösung dar.:1. Einleitung 1.1. Problemstellung und Motivation der Arbeit 1.2. Stand der Technik 1.2.1. Hochgeschwindigkeitsfräsen 1.2.2. Verwendung lang kragender Schaftfräser 1.3. Thema und Aufbau der Arbeit 2. Theoretische Grundlagen 2.1. Kontinuumsmechanische Grundbegriffe 2.2. Spannungen und konstitutive Gleichungen 2.3. Prinzip von Hamilton 2.4. Lösungstheorie 2.4.1. Anfangswertprobleme 2.4.2. Randwertprobleme 2.5. Stochastische Grundbegriffe 3. Balkenmodelle 3.1. Verformungskinematik des Balkens 3.2. Variationsformulierung 3.3. Modellierung der Unwucht 3.4. Globale Diskretisierung 3.4.1. Stationäre Lage und Linearisierung 3.4.2. Ortsfunktionen 3.5. Lokale Diskretisierung 3.6. Anmerkungen zur schubweichen Formulierung 3.7. Berechnungsergebnisse 3.7.1. Ruhendes Werkzeug 3.7.2. Rotierendes Werkzeug 4. Schalenmodelle 4.1. Verformungskinematik der Schale 4.2. Variationsformulierung 4.3. Globale Diskretisierung 4.3.1. Stationäre Lage und Linearisierung 4.4. Lokale Diskretisierung mittels FEM 4.4.1. Konforme flache Schalenelemente 4.5. Anmerkungen zur schubweichen Formulierung 4.6. Berechnungsergebnisse 4.6.1. Ruhender Schaft 4.6.2. Rotierender Schaft 5. Zusammenfassung und Ausblick 6. Verzeichnisse 6.1. Quellenverzeichnis 6.2. Symbolverzeichnis 6.3. Abbildungsverzeichnis 6.4. Tabellenverzeichnis A. Feldgleichungen und Ableitungen der Ansätze für die Balkenmodelle B. Anmerkungen zum Timoshenko-Balken C. Feldgleichungen und Ableitungen der Ansätze für die Schalenmodelle D. Anmerkungen zur Mindlin-Reissner-Schale / The use of long slender end mills for high-speed-cutting (HSC) holds special requirements with respect to the system dynamics. In this context, several tool models in the area of continuum mechanics are presented. Especially hollow tool shafts, with a fluid medium inside, for the purpose of automatic balancing are considered. Starting with the kinematics of deformation, Hamilton's principle is used to evaluate the variational formulation. Therefore, also the influence of a stochastic distributed unbalance, geometrical nonlinearities and shear deformations are discussed. For space discretisation local as well as global approaches are used and compared with each other. Following up on this, results are presented, which show the influence of different geometrical and process-related parameters due to the eigenfrequencies, stationary deformation, stability and time solution.:1. Einleitung 1.1. Problemstellung und Motivation der Arbeit 1.2. Stand der Technik 1.2.1. Hochgeschwindigkeitsfräsen 1.2.2. Verwendung lang kragender Schaftfräser 1.3. Thema und Aufbau der Arbeit 2. Theoretische Grundlagen 2.1. Kontinuumsmechanische Grundbegriffe 2.2. Spannungen und konstitutive Gleichungen 2.3. Prinzip von Hamilton 2.4. Lösungstheorie 2.4.1. Anfangswertprobleme 2.4.2. Randwertprobleme 2.5. Stochastische Grundbegriffe 3. Balkenmodelle 3.1. Verformungskinematik des Balkens 3.2. Variationsformulierung 3.3. Modellierung der Unwucht 3.4. Globale Diskretisierung 3.4.1. Stationäre Lage und Linearisierung 3.4.2. Ortsfunktionen 3.5. Lokale Diskretisierung 3.6. Anmerkungen zur schubweichen Formulierung 3.7. Berechnungsergebnisse 3.7.1. Ruhendes Werkzeug 3.7.2. Rotierendes Werkzeug 4. Schalenmodelle 4.1. Verformungskinematik der Schale 4.2. Variationsformulierung 4.3. Globale Diskretisierung 4.3.1. Stationäre Lage und Linearisierung 4.4. Lokale Diskretisierung mittels FEM 4.4.1. Konforme flache Schalenelemente 4.5. Anmerkungen zur schubweichen Formulierung 4.6. Berechnungsergebnisse 4.6.1. Ruhender Schaft 4.6.2. Rotierender Schaft 5. Zusammenfassung und Ausblick 6. Verzeichnisse 6.1. Quellenverzeichnis 6.2. Symbolverzeichnis 6.3. Abbildungsverzeichnis 6.4. Tabellenverzeichnis A. Feldgleichungen und Ableitungen der Ansätze für die Balkenmodelle B. Anmerkungen zum Timoshenko-Balken C. Feldgleichungen und Ableitungen der Ansätze für die Schalenmodelle D. Anmerkungen zur Mindlin-Reissner-Schale

In der vorliegenden Dissertation wird auf der Grundlage einer kontinuumsmechanischen Herangehensweise die Formulierung eines mechanischen Modells in Verbindung mit der Umsetzung in ein Schalenelement im Rahmen der Finite-Element-Methode zur Simulation des Tragverhaltens geschichteter Flächentragwerksstrukturen unter Berücksichtigung der Schädigungsart Delamination vorgestellt. Grundlage des Mehrschichten-Modells ist die Entwicklung einer geometrisch nichtlinearen oberflächenorientierten Schalentheorie mit schub- und dickenelastischem Verhalten ausgehend von der vollständigen Kinematik einer Multidirektor-Theorie. Der Oberflächenbezug gewährleistet eine auf Kontaktprobleme angepasste mechanische Modellbildung. Innerhalb der Schichten wurde ein Materialgesetz für linear elastisches, orthotropes Material verwendet, dessen Dreidimensionalität durch die Schalenformulierung nicht eingeschränkt wird. Das Hauptaugenmerk der Arbeit liegt auf der Entwicklung eines auf verschiedene Materialien anpassbaren Schichten-Verbundmodells. Das Versagen des Schichtenverbundes – Delamination genannt – wurde durch ein einfaches Spannungskriterium beschrieben. Die Delamination wird durch Modifikation der kinematischen Bedingungen diskret berücksichtigt. Zur Beschreibung des Tragverhaltens nach Ausbildung der Delamination wurde ein als „innerer Kontakt“ bezeichnetes Kontakt-Modell entwickelt, durch das Adhäsion zwischen den Schichten berücksichtigt werden kann. Das vorliegende Schalenmodell kann bei Berücksichtigung von Delamination auf Probleme, in denen kleine Relativverschiebungen zu erwarten sind, für beliebige elastische Materialien angewendet werden. Der Rahmen, in dem diese Arbeit entstand, gab den hauptsächlichen Einsatzbereich, die Simulation von Flächentragwerksstrukturen mit einer Verstärkungsschicht aus textilbewehrtem Feinbeton, vor. / This publication introduces, in a continuum-mechanical approach, the formulation of a mechanical model in connection with the transformation into a shell element using the finite element method for the simulation of the load-bearing behaviour of laminated shell structures thereby considering delamination as a type of damage. This multi-layer model is based upon the development of a geometrically nonlinear surface-related shell theory with shear-elastic behaviour and variable thickness, beginning with the complete kinematics of a multi-director theory. The surface relationship ensures a mechanical modelling which is adaptable for contact problems. A linear-elastic orthotropic material law, whose three-dimensionality is not restricted by the shell formulation, applies within the layers. The main focus of the thesis is on the development of a layer-bond model that can be adjusted for different materials. The debonding of layers – called delamination – is described by a simple stress criterion. Delamination is discretely taken into account by modifying the kinematic conditions. A contact model, called „inner contact“, that can be used to account for adhesion between layers, has been developed to describe the load-bearing behaviour after delamination has occurred. The present shell model is restricted to elastic material behaviour and can preferably be applied to such problems where small relative displacements are expected. The environment, in which this research has been conducted, established the primary of application area, which is the simulation of shell structures within a strengthening layer comprised of textile-reinforced concrete.

Der Gegenstand der vorliegenden Arbeit ist die Simulation von Betonstrukturen beliebiger Geometrie unter überwiegender Zugbelastung. Die Modellierung erfolgt auf Makroebene als Kontinuum und zur Lösung des mechanischen Feldproblems wird die Finite-Elemente-Methode verwendet. Ein neues Materialmodell für Beton und eine Erweiterung der Bruchenergetischen Regularisierung werden vorgestellt. Die Arbeit ist in zwei Teile gegliedert. Im ersten Teil wird ein lokales, anisotropes Schädigungsmodell abgeleitet, wobei als Schädigungsvariable ein symmetrischer Tensor zweiter Stufe gewählt wird. Die Verwendung einer Normalenregel im Raum der dissipativen Kräfte zur Bestimmung der Schädigungsevolution und die Definition der Schädigungsgrenzflächen im Raum der dissipativen Kräfte gewährleisten die Gültigkeit der Hauptsätze der Thermodynamik und des Prinzips der maximalen Dissipationsrate. Vorteilhaft ist die Symmetrie der Materialtangente, die sich aus diesem Vorgehen ergibt. Eine Formulierung mit drei entkoppelten Schädigungsgrenzflächen wird vorgeschlagen. Eine wichtige Forderung bei der Ableitung des Materialmodells war die Verwendung einer möglichst geringen Anzahl von Materialparametern, welche darüber hinaus aus wenigen Standardversuchen bestimmbar sein sollten. Das Schädigungsmodell enthält als Materialparameter den Elastizitätsmodul, die Querdehnzahl, die Zugfestigkeit und die auf eine Einheitsfläche bezogene Bruchenergie. Im zweiten Teil der Arbeit stehen Lokalisierung und Regularisierung im Fokus der Betrachtungen. Aufgrund der lokalen Formulierung des Materialmodells tritt bei Finite-Elemente Simulationen eine Netzabhängigkeit der Simulationsergebnisse auf. Um dieser Problematik zu begegnen und netzunabhängige Simulationen zu erreichen, werden Regularisierungstechniken angewendet. In dieser Arbeit wird die Bruchenergetische Regularisierung eingesetzt, die durch die Einführung einer äquivalenten Breite in ein lokal formuliertes Stoffgesetz gekennzeichnet ist. Die spezielle Wahl eines Wertes für die äquivalente Breite beruht auf der Forderung, dass in der Simulation die korrekte Bruchenergie je Einheitsfläche für den Bruchprozess verbraucht wird, d.h. die Energiedissipation der Realität entspricht. In vorliegender Arbeit wird die neue These aufgestellt, dass die Energiedissipation nur für den Fall korrekt abgebildet wird, wenn die im Stoffgesetz enthaltene äquivalente Breite in jedem Belastungsinkrement der Breite des Bereiches entspricht, in dem in der Simulation Energie dissipiert wird. In einer Simulation wird in den Bereichen Energie dissipiert, in denen die Schädigung im aktuellen Belastungsinkrement zunimmt. In vorliegender Arbeit werden die energiedissipierenden Bereiche daher als Pfad der Schädigungsrate bezeichnet. Um Erkenntnisse über die Entwicklung des Pfades der Schädigungsrate über den Belastungsverlauf zu erhalten, wurden umfangreiche Untersuchungen anhand von Simulationen eines beidseitig gekerbten Betonprobekörpers unter kombinierter Zug-Schubbeanspruchung durchgeführt, wobei die gewählten Werte für die äquivalente Breite variiert wurden. Es wurde stets eine Diskretisierung mit linearen Verschiebungselementen verwendet, wobei die Bereiche mit zu erwartender Schädigung feiner und regelmäßig mit Elementen quadratischer Geometrie diskretisiert wurden. Die Ergebnisse der Untersuchungen zeigen, dass die Breite des Pfades der Schädigungsrate abhängig ist von der Schädigung am betrachteten Materialpunkt, dem von Schädigungsrichtung und Elementkante eingeschlossenen Winkel, der Elementgröße und den Materialparametern. Um die geforderte Übereinstimmung von äquivalenter Breite und der Breite des Pfades der Schädigungsrate zu erreichen, werden neue Ansätze für die äquivalente Breite vorgeschlagen, die die erwähnten Einflüsse berücksichtigen. Simulationen unter Verwendung der neuen Ansätze für die äquivalente Breite führen zu einer guten Übereinstimmung von äquivalenter Breite und der Breite des Pfades der Schädigungsrate in der Simulation. Die Ergebnisse der Simulationen, wie z.B. Last-Verformungsbeziehung und Rissverläufe, sind netzunabhängig und stimmen gut mit den experimentellen Beobachtungen überein. Basierend auf den gewonnenen Erkenntnissen wird eine Erweiterung der Bruchenergetischen Regularisierung vorgeschlagen: die Adaptive Bruchenergetische Regularisierung. Im abschließenden Kapitel der Arbeit werden mit der vorgeschlagenen Theorie, dem neuen Schädigungsmodell und der Adaptiven Bruchenergetischen Regularisierung, noch zwei in der Literatur gut dokumentierte Versuche simuliert. Die Simulationsergebnisse entsprechen den experimentellen Beobachtungen. / This doctoral thesis deals with the simulation of predominantly tensile loaded plain concrete structures. Concrete is modeled on the macro level and the Finite Element Method is applied to solve the resulting mechanical field problem. A new material model for concrete based on continuum damage mechanics and an extended regularization technique based on the fracture energy approach are presented. The thesis is subdivided into two parts. In the first part, a local, anisotropic damage model for concrete is derived. This model uses a symmetric second-order tensor as the damage variable, which enables the simulation of orthotropic degradation. The validity of the first and the second law of thermodynamics as well as the validity of the principle of maximum dissipation rate are required. Using a normal rule in the space of the dissipative forces, which are the thermodynamically conjugated variables to the damage variables, and the definition of the loading functions in the space of the dissipative forces guarantee their validity. The suggested formulation contains three decoupled loading functions. A further requirement in the derivation of the model was the minimization of the number of material parameters, which should be determined by a small number of standard experiments. The material parameters of the new damage model are the Young’s modulus, the Poisson’s ratio, the tensile strength and the fracture energy per unit area. The second part of the work focuses on localization and regularization. If a Finite Element simulation is performed using a local material model for concrete, the results of the Finite Element simulation are mesh-dependent. To attain mesh-independent simulations, a regularization technique must be applied. The fracture energy approach, which is characterized by introducing a characteristic length in a locally formulated material model, is used as regularization technique in this work. The choice of a value for the characteristic length is founded by the requirement, that the fracture energy per unit area, which is consumed for the fracture process in the simulation, must be the same as in experiment, i.e. the energy dissipation must be correct. In this dissertation, the new idea is suggested that the correct energy dissipation can be only attained if the characteristic length in the material model coincides in every loading increment with the width of the energy-dissipating zone in the simulation. The energy-dissipating zone in a simulation is formed by the integration points with increasing damage and obtains the name: damage rate path. Detailed investigations based on simulations of a double-edge notched specimen under mixed-mode loading are performed with varying characteristic lengths in order to obtain information concerning the evolution of the damage rate path during a simulation. All simulations were performed using displacement-based elements with four nodes. The range with expected damage was always finer and regularly discretized. The results of the simulations show that the width of the damage rate path depends on the damage at the specific material point, on the angle between damage direction and element edges, on the element size and on the material parameters. Based on these observations, new approaches for the characteristic length are suggested in order to attain the coincidence of the characteristic length with the width of the damage rate path. Simulations by using the new approaches yield a sufficient coincidence of the characteristic length with the width of the damage rate path. The simulations are mesh-independent and the results of the simulation, like load-displacement curves or crack paths, correspond to the experimental results. Based on all new information concerning the regularization technique, an extension of the fracture energy approach is suggested: the adaptive fracture energy approach. The validity and applicability of the suggested theory, the new anisotropic damage model and the adaptive fracture energy approach, are verified in the final chapter of the work with simulations of two additional experiments, which are well documented in the literature. The results of the simulations correspond to the observations in the experiments.

Zur Analyse der Beanspruchungen in textilverstärkten Elastomerbauteilen wie Luftfedern, Reifen, Riemen und Schläuchen sind Berechnungsmodelle mit einer feinen Balance zwischen Genauigkeit und Effizienz erforderlich. Die großen Deformationen, stark anisotropen Struktureigenschaften und kleinen Abmessungen der Festigkeitsträger gegenüber denen des Bauteils bedürfen einerseits einer detaillierten Modellierung, andererseits sind die kritischen Bereiche in diesen Bauteilen räumlich stark begrenzt, sodass eine Reduktion des Berechnungsaufwands erstrebenswert ist. Diese Modellreduktion führt zu Simulationen mit geringer Rechenzeit, die für eine praxistaugliche Optimierung von Hybridcord-Elastomer-Verbunden unerlässlich sind. Die beiden Hauptschwerpunkte der vorliegenden Arbeit bilden die kontinuumsmechanische Modellierung von Hybridcorden und die Erstellung repräsentativer Volumenelemente hochbeanspruchter Hybridcord-Elastomer-Verbunde. Aufbauend auf einem anisotropen Plastizitätsmodell zur Erfassung der Reibung in Multifilamentgarnen stellt ein Finite-Elemente-Modell zur Simulation der Verzwirnung von Hybridcorden das Fundament der Arbeit dar. Anhand experimenteller Ergebnisse aus Zug- und Torsionsversuchen sowie einem Vergleich mit Querschnittsaufnahmen wird gezeigt, dass das Modell die komplexen Eigenschaften eines Hybridcords abbilden kann. Die Grundlage der repräsentativen Volumenelemente stellt eine Erweiterung der klassischen periodischen Randbedingungen dar, die eine Berücksichtigung von Krümmungen und Drucklasten ermöglicht. Das Modell eignet sich daher, die Beanspruchungen in den hochbelasteten Bereichen textilverstärkter Elastomerbauteile wie der Rollfalte einer Luftfeder effizient zu analysieren. Mittels Parameterstudien werden abschließend Hybridcorde und Hybridcord-Elastomer-Verbunde untersucht und einige Hinweise für eine optimale Gestaltung hinsichtlich minimaler Beanspruchungen des Elastomers, des Hybridcords sowie der Grenzfläche gegeben. / The analysis of stresses and strains within textile-reinforced rubber components like air springs, tyres, driving belts, and tubes requires accurate as well as efficient computational models. On the one hand, the large deformations, the composite's strongly anisotropic properties, and the large ratio between the size of the cords and the composite necessitate a precise modeling. On the other hand, the highly loaded parts of the components are spatially confined and thus a reduction of the computational effort is desirable. These reduced models are efficient enough for performing engineering-oriented optimizations. The two main priorities of this work are the continuum mechanical modeling of hybrid cords and the development of representative volume elements of highly loaded hybrid cord-rubber composites. Based on an anisotropic plasticity model, which takes the frictional sliding between the filaments within multifilament yarns into account, a finite element model for the simulation of the twisting process of a hybrid cord is the fundament of this work. A comparison with experimental results from tensile and torsional tests as well as images of cross sections validate the proposed hybrid cord model. The basis of the representative volume element is the extension of the classical periodic boundary conditions, which now enable to take the curvature and pressure load into account. Thus, the model is suitable to analyze the highly loaded parts of hybrid cord-rubber composites like the rolling lobe of an air spring. Finally, the set-ups of hybrid cords and hybrid cord-rubber composites are analyzed by means of parameter studies to obtain a minimized loading of the rubber, yarns, and their interface.

Zur Analyse der Beanspruchungen in textilverstärkten Elastomerbauteilen wie Luftfedern, Reifen, Riemen und Schläuchen sind Berechnungsmodelle mit einer feinen Balance zwischen Genauigkeit und Effizienz erforderlich. Die großen Deformationen, stark anisotropen Struktureigenschaften und kleinen Abmessungen der Festigkeitsträger gegenüber denen des Bauteils bedürfen einerseits einer detaillierten Modellierung, andererseits sind die kritischen Bereiche in diesen Bauteilen räumlich stark begrenzt, sodass eine Reduktion des Berechnungsaufwands erstrebenswert ist. Diese Modellreduktion führt zu Simulationen mit geringer Rechenzeit, die für eine praxistaugliche Optimierung von Hybridcord-Elastomer-Verbunden unerlässlich sind. Die beiden Hauptschwerpunkte der vorliegenden Arbeit bilden die kontinuumsmechanische Modellierung von Hybridcorden und die Erstellung repräsentativer Volumenelemente hochbeanspruchter Hybridcord-Elastomer-Verbunde. Aufbauend auf einem anisotropen Plastizitätsmodell zur Erfassung der Reibung in Multifilamentgarnen stellt ein Finite-Elemente-Modell zur Simulation der Verzwirnung von Hybridcorden das Fundament der Arbeit dar. Anhand experimenteller Ergebnisse aus Zug- und Torsionsversuchen sowie einem Vergleich mit Querschnittsaufnahmen wird gezeigt, dass das Modell die komplexen Eigenschaften eines Hybridcords abbilden kann. Die Grundlage der repräsentativen Volumenelemente stellt eine Erweiterung der klassischen periodischen Randbedingungen dar, die eine Berücksichtigung von Krümmungen und Drucklasten ermöglicht. Das Modell eignet sich daher, die Beanspruchungen in den hochbelasteten Bereichen textilverstärkter Elastomerbauteile wie der Rollfalte einer Luftfeder effizient zu analysieren. Mittels Parameterstudien werden abschließend Hybridcorde und Hybridcord-Elastomer-Verbunde untersucht und einige Hinweise für eine optimale Gestaltung hinsichtlich minimaler Beanspruchungen des Elastomers, des Hybridcords sowie der Grenzfläche gegeben.:Inhaltsverzeichnis Abkürzungs- und Symbolverzeichnis VIII 1 Einleitung 1 2 Grundlagen der Mathematik und der Mechanik 6 2.1 Tensoralgebra und -analysis 6 2.2 Nichtlineare Kontinuumsmechanik 11 2.3 Nichtlineare Finite-Elemente-Methode 16 3 Einordnung in den Stand der Forschung 22 4 Experimentelle Untersuchungen 26 4.1 Charakterisierung der Standardcorde 26 4.2 Charakterisierung der Hybridcorde 33 5 Materialmodelle für Multi lamentgarne 38 5.1 Anisotropes Plastizitätsmodell der Filamentreibung 38 5.2 Numerische Lösung der Materialgleichungen 43 5.3 Analytische Lösung für reibungsfreies Gleiten 48 5.4 Modellierung des thermischen Schrumpfens 50 6 FEM-basierte Modellierung von Hybridcorden 53 6.1 Simulation der Verzwirnung eines Standardcords 53 6.2 Erweiterung des Berechnungsmodells auf Hybridcorde 60 6.3 Analytisches Modell der Geometrie eines Hybridcords 65 6.4 Qualitative Charakterisierung des Hybridcordmodells 74 6.5 Parameteridenti kation und Validierung 83 6.6 Optimierungsbeispiele 92 7 Schalenartige RVEs für Cord-Elastomer-Verbunde 96 7.1 Geometrie der Axial- und der Kreuzlage 96 7.2 Erweiterte periodische Randbedingungen 98 7.3 E ektive Schaleneigenschaften 111 7.4 Berücksichtigung der Drucklast 118 7.5 Diskretisierung der RVEs 122 7.6 Submodelltechnik 128 7.7 Parameterstudien an Hybridcord-Elastomer-Verbunden 135 8 Zusammenfassung und Ausblick 146 Literaturverzeichnis 151 / The analysis of stresses and strains within textile-reinforced rubber components like air springs, tyres, driving belts, and tubes requires accurate as well as efficient computational models. On the one hand, the large deformations, the composite's strongly anisotropic properties, and the large ratio between the size of the cords and the composite necessitate a precise modeling. On the other hand, the highly loaded parts of the components are spatially confined and thus a reduction of the computational effort is desirable. These reduced models are efficient enough for performing engineering-oriented optimizations. The two main priorities of this work are the continuum mechanical modeling of hybrid cords and the development of representative volume elements of highly loaded hybrid cord-rubber composites. Based on an anisotropic plasticity model, which takes the frictional sliding between the filaments within multifilament yarns into account, a finite element model for the simulation of the twisting process of a hybrid cord is the fundament of this work. A comparison with experimental results from tensile and torsional tests as well as images of cross sections validate the proposed hybrid cord model. The basis of the representative volume element is the extension of the classical periodic boundary conditions, which now enable to take the curvature and pressure load into account. Thus, the model is suitable to analyze the highly loaded parts of hybrid cord-rubber composites like the rolling lobe of an air spring. Finally, the set-ups of hybrid cords and hybrid cord-rubber composites are analyzed by means of parameter studies to obtain a minimized loading of the rubber, yarns, and their interface.:Inhaltsverzeichnis Abkürzungs- und Symbolverzeichnis VIII 1 Einleitung 1 2 Grundlagen der Mathematik und der Mechanik 6 2.1 Tensoralgebra und -analysis 6 2.2 Nichtlineare Kontinuumsmechanik 11 2.3 Nichtlineare Finite-Elemente-Methode 16 3 Einordnung in den Stand der Forschung 22 4 Experimentelle Untersuchungen 26 4.1 Charakterisierung der Standardcorde 26 4.2 Charakterisierung der Hybridcorde 33 5 Materialmodelle für Multi lamentgarne 38 5.1 Anisotropes Plastizitätsmodell der Filamentreibung 38 5.2 Numerische Lösung der Materialgleichungen 43 5.3 Analytische Lösung für reibungsfreies Gleiten 48 5.4 Modellierung des thermischen Schrumpfens 50 6 FEM-basierte Modellierung von Hybridcorden 53 6.1 Simulation der Verzwirnung eines Standardcords 53 6.2 Erweiterung des Berechnungsmodells auf Hybridcorde 60 6.3 Analytisches Modell der Geometrie eines Hybridcords 65 6.4 Qualitative Charakterisierung des Hybridcordmodells 74 6.5 Parameteridenti kation und Validierung 83 6.6 Optimierungsbeispiele 92 7 Schalenartige RVEs für Cord-Elastomer-Verbunde 96 7.1 Geometrie der Axial- und der Kreuzlage 96 7.2 Erweiterte periodische Randbedingungen 98 7.3 E ektive Schaleneigenschaften 111 7.4 Berücksichtigung der Drucklast 118 7.5 Diskretisierung der RVEs 122 7.6 Submodelltechnik 128 7.7 Parameterstudien an Hybridcord-Elastomer-Verbunden 135 8 Zusammenfassung und Ausblick 146 Literaturverzeichnis 151

The classical Cauchy-Boltzmann theory of continuum mechanics requires that the dimension, over which macroscopic gradients occur, are much larger than characteristic length scales of the microstructure. For this reason, the classical continuum theory comes to its limits for very small specimens or if material degradation leads to a localisation of deformations into bands, whose width is determined by the microstructure itself. Deviations from the predictions of the classical theory of continuum mechanics are referred to as size effects. It is well-known, that generalised continuum theories can describe size effects in principle. Especially micromorphic theories gain increasing popularity due its favorable numerical implementation. However, the formulation of the additionally necessary constitutive equations is a problem. For linear-elastic behavior, the number of material parameters increases considerably compared to the classical theory. The experimental determination of these parameters is thus very difficult. For nonlinear and history-dependent processes, even the qualitative structure of the constitutive equations can hardly be assessed solely on base of phenomenological considerations. Homogenisation methods are a promising approach to solve this problem. The present thesis starts with a critical review on the classical theory of homogenisation and the approaches on micromorphic homogenisation which are available in literature. On this basis, a theory is developed for the homogenisation of a classical Cauchy-Boltzmann continuum at the microscale towards a micromorphic continuum at the macroscale. In particular, the micro-macro-relations are specified for all macroscopic kinetic and kinematic field quantities. On the microscale, the corresponding boundary-value problem is formulated, whereby kinematic, static or periodic boundary conditions can be used. No restrictions are imposed on the material behavior, i. e. it can be linear or nonlinear. The special cases of the micropolar theory (Cosserat theory), microstrain theory and microdilatational theorie are considered. The proposed homogenisation method is demonstrated for several examples. The simplest example is the uniaxial case, for which the exact solution can be specified. Furthermore, the micromorphic elastic properties of a porous, foam-like material are estimated in closed form by means of Ritz' method with a cubic ansatz. A comparison with partly available exact solutions and FEM solutions indicates a qualitative and quantitative agreement of sufficient accuracy. For the special cases of micropolar and microdilatational theory, the material parameters are specified in the established nomenclature from literature. By means of these material parameters the size effect of an elastic foam structure is investigated and compared with corresponding results from literature. Furthermore, micromorphic damage models for quasi-brittle and ductile failure are presented. Quasi-brittle damage is modelled by propagation of microcracks. For the ductile mechanism, Gurson's limit-load approach on the microscale is extended by microdilatational terms. A finite-element implementation shows, that the damage model exhibits h-convergence even in the softening regime and that it thus can describe localisation.:1 Introduction 2 Literature review: Micromorphic theory and strain-gradient theory 2.1 Variational approach 2.1.1 Cauchy-Boltzmann continuum 2.1.2 Second gradient theory / Strain gradient theory 2.1.3 Micromorphic theory 2.1.4 Method of virtual power 2.2 Homogenisation approaches 2.2.1 Classical theory of homogenisation 2.2.2 Strain-gradient theory by Gologanu, Kouznetsova et al. 2.2.3 Micromorphic theory by Eringen 2.2.4 Average field theory by Forest et al. 2.3 Scope of the present thesis 3 Homogenisation towards a micromorphic continuum 3.1 Thermodynamic considerations and generalized Hill-Mandel lemma 3.2 Surface operator and kinetic micro-macro relations 3.3 Kinematic micro-macro relations 3.4 Porous material 3.5 Kinematic and periodic boundary conditions 3.6 Special cases 3.6.1 Strain-gradient theory / Second gradient theory 3.6.2 Micropolar theory 3.6.3 Microstrain theory 3.6.4 Microdilatational theory 4 Elastic Behaviour 4.1 Uniaxial case 4.2 Upper bound estimates by Ritz' Method 4.3 Isotropic porous material 4.4 Micropolar theory 4.5 Microdilatational theory 4.6 Size effect in simple shear 5 Damage Models 5.1 Quasi-brittle damage 5.2 Microdilatational extension of Gurson’s model of ductile damage 5.2.1 Limit load analysis for rigid ideal-plastic material 5.2.2 Phenomenological extensions 5.2.3 FEM implementation 5.2.4 Example 6 Discussion / Die klassische Cauchy-Boltzmann-Kontinuumstheorie setzt voraus, dass die Abmessungen, über denen makroskopische Gradienten auftreten, sehr viele größer sind als charakteristische Längenskalen der Mikrostruktur. Aus diesem Grund stößt die klassische Kontinuumstheorie bei sehr kleinen Proben ebenso an ihre Grenzen wie bei Schädigungsvorgängen, bei denen die Deformationen in Bändern lokalisieren, deren Breite selbst von der Längenskalen der Mikrostruktur bestimmt wird. Abweichungen von Vorhersagen der klassischen Kontinuumstheorie werden als Größeneffekte bezeichnet. Es ist bekannt, dass generalisierte Kontinuumstheorien Größeneffekte prinzipiell beschreiben können. Insbesondere mikromorphe Theorien erfreuen sich auf Grund ihrer vergleichsweise einfachen numerischen Implementierung wachsender Beliebtheit. Ein großes Problem stellt dabei die Formulierung der zusätzlich notwendigen konstitutiven Gleichungen dar. Für linear-elastisches Verhalten steigt die Zahl der Materialparameter im Vergleich zur klassischen Theorie stark an, was deren experimentelle Bestimmung sehr schwierig macht. Bei nichtlinearen und lastgeschichtsabhängigen Prozessen lässt sich selbst die qualitative Struktur der konstitutiven Gleichungen ausschließlich auf Basis phänomenologischer Überlegungen kaum erschließen. Homogenisierungsverfahren stellen einen vielversprechenden Ansatz dar, um dieses Problem zu lösen. Die vorliegende Arbeit gibt zunächst einen kritischen Überblick über die klassische Theorie der Homogenisierung sowie die im Schrifttum verfügbaren Ansätze zur mikromorphen Homogenisierung. Auf dieser Basis wird eine Theorie zur Homogenisierung eines klassischen Cauchy-Boltzmann-Kontinuums auf Mikroebene zu einem mikromorphen Kontinuum auf der Makroebene entwickelt. Insbesondere werden Mikro-Makro-Relationen für alle makroskopischen kinetischen und kinematischen Feldgrößen angegebenen. Auf der Mikroebene wird das entsprechende Randwertproblem formuliert, wobei kinematische, statische oder periodische Randbedingungen verwendet werden können. Das Materialverhalten unterliegt keinen Einschränkungen, d. h., dass es sowohl linear als auch nichtlinear sein kann. Die Sonderfälle der mikropolaren Theorie (Cosserat-Theorie), Mikrodehnungstheorie und mikrodilatationalen Theorie werden erarbeitet. Das vorgeschlagene Homogenisierungsverfahren wird für eine Reihe von Beispielen demonstriert. Als einfachstes Beispiel dient der einachsige Fall, für den die exakte Lösung angegebenen werden kann. Weiterhin werden die mikromorphen, elastischen Eigenschaften eines porösen, schaumartigen Materials mittels des Ritz-Verfahrens mit einem kubischen Ansatz in geschlossener Form abgeschätzt. Ein Vergleich mit teilweise verfügbaren exakten Lösungen sowie FEM-Lösungen weist eine qualitative und quantitative Übereinstimmung hinreichender Genauigkeit aus. Für die Sonderfälle mikropolaren und mikrodilatationalen Theorien werden die Materialparameter in der im Schrifttum üblichen Nomenklatur angegebenen. Mittels dieser Materialparameter wird der Größeneffekt in einer elastischen Schaumstruktur untersucht und mit entsprechenden Ergebnissen aus dem Schrifttum verglichen. Desweiteren werden mikromorphe Schädigungsmodelle für quasi-sprödes und duktiles Versagen vorgestellt. Quasi-spröde Schädigung wird durch das Wachstum von Mikrorissen modelliert. Für den duktilen Mechanismus wird der Ansatz von Gurson einer Grenzlastanalyse auf Mikroebene um mikrodilatationale Terme erweitert. Eine Finite-Elemente-Implementierung zeigt, dass das Schädigungsmodell auch im Entfestigungsbereich h-Konvergenz aufweist und die Lokalisierung beschreiben kann.:1 Introduction 2 Literature review: Micromorphic theory and strain-gradient theory 2.1 Variational approach 2.1.1 Cauchy-Boltzmann continuum 2.1.2 Second gradient theory / Strain gradient theory 2.1.3 Micromorphic theory 2.1.4 Method of virtual power 2.2 Homogenisation approaches 2.2.1 Classical theory of homogenisation 2.2.2 Strain-gradient theory by Gologanu, Kouznetsova et al. 2.2.3 Micromorphic theory by Eringen 2.2.4 Average field theory by Forest et al. 2.3 Scope of the present thesis 3 Homogenisation towards a micromorphic continuum 3.1 Thermodynamic considerations and generalized Hill-Mandel lemma 3.2 Surface operator and kinetic micro-macro relations 3.3 Kinematic micro-macro relations 3.4 Porous material 3.5 Kinematic and periodic boundary conditions 3.6 Special cases 3.6.1 Strain-gradient theory / Second gradient theory 3.6.2 Micropolar theory 3.6.3 Microstrain theory 3.6.4 Microdilatational theory 4 Elastic Behaviour 4.1 Uniaxial case 4.2 Upper bound estimates by Ritz' Method 4.3 Isotropic porous material 4.4 Micropolar theory 4.5 Microdilatational theory 4.6 Size effect in simple shear 5 Damage Models 5.1 Quasi-brittle damage 5.2 Microdilatational extension of Gurson’s model of ductile damage 5.2.1 Limit load analysis for rigid ideal-plastic material 5.2.2 Phenomenological extensions 5.2.3 FEM implementation 5.2.4 Example 6 Discussion

Der Gegenstand der vorliegenden Arbeit ist die Simulation von Betonstrukturen beliebiger Geometrie unter überwiegender Zugbelastung. Die Modellierung erfolgt auf Makroebene als Kontinuum und zur Lösung des mechanischen Feldproblems wird die Finite-Elemente-Methode verwendet. Ein neues Materialmodell für Beton und eine Erweiterung der Bruchenergetischen Regularisierung werden vorgestellt. Die Arbeit ist in zwei Teile gegliedert. Im ersten Teil wird ein lokales, anisotropes Schädigungsmodell abgeleitet, wobei als Schädigungsvariable ein symmetrischer Tensor zweiter Stufe gewählt wird. Die Verwendung einer Normalenregel im Raum der dissipativen Kräfte zur Bestimmung der Schädigungsevolution und die Definition der Schädigungsgrenzflächen im Raum der dissipativen Kräfte gewährleisten die Gültigkeit der Hauptsätze der Thermodynamik und des Prinzips der maximalen Dissipationsrate. Vorteilhaft ist die Symmetrie der Materialtangente, die sich aus diesem Vorgehen ergibt. Eine Formulierung mit drei entkoppelten Schädigungsgrenzflächen wird vorgeschlagen. Eine wichtige Forderung bei der Ableitung des Materialmodells war die Verwendung einer möglichst geringen Anzahl von Materialparametern, welche darüber hinaus aus wenigen Standardversuchen bestimmbar sein sollten. Das Schädigungsmodell enthält als Materialparameter den Elastizitätsmodul, die Querdehnzahl, die Zugfestigkeit und die auf eine Einheitsfläche bezogene Bruchenergie. Im zweiten Teil der Arbeit stehen Lokalisierung und Regularisierung im Fokus der Betrachtungen. Aufgrund der lokalen Formulierung des Materialmodells tritt bei Finite-Elemente Simulationen eine Netzabhängigkeit der Simulationsergebnisse auf. Um dieser Problematik zu begegnen und netzunabhängige Simulationen zu erreichen, werden Regularisierungstechniken angewendet. In dieser Arbeit wird die Bruchenergetische Regularisierung eingesetzt, die durch die Einführung einer äquivalenten Breite in ein lokal formuliertes Stoffgesetz gekennzeichnet ist. Die spezielle Wahl eines Wertes für die äquivalente Breite beruht auf der Forderung, dass in der Simulation die korrekte Bruchenergie je Einheitsfläche für den Bruchprozess verbraucht wird, d.h. die Energiedissipation der Realität entspricht. In vorliegender Arbeit wird die neue These aufgestellt, dass die Energiedissipation nur für den Fall korrekt abgebildet wird, wenn die im Stoffgesetz enthaltene äquivalente Breite in jedem Belastungsinkrement der Breite des Bereiches entspricht, in dem in der Simulation Energie dissipiert wird. In einer Simulation wird in den Bereichen Energie dissipiert, in denen die Schädigung im aktuellen Belastungsinkrement zunimmt. In vorliegender Arbeit werden die energiedissipierenden Bereiche daher als Pfad der Schädigungsrate bezeichnet. Um Erkenntnisse über die Entwicklung des Pfades der Schädigungsrate über den Belastungsverlauf zu erhalten, wurden umfangreiche Untersuchungen anhand von Simulationen eines beidseitig gekerbten Betonprobekörpers unter kombinierter Zug-Schubbeanspruchung durchgeführt, wobei die gewählten Werte für die äquivalente Breite variiert wurden. Es wurde stets eine Diskretisierung mit linearen Verschiebungselementen verwendet, wobei die Bereiche mit zu erwartender Schädigung feiner und regelmäßig mit Elementen quadratischer Geometrie diskretisiert wurden. Die Ergebnisse der Untersuchungen zeigen, dass die Breite des Pfades der Schädigungsrate abhängig ist von der Schädigung am betrachteten Materialpunkt, dem von Schädigungsrichtung und Elementkante eingeschlossenen Winkel, der Elementgröße und den Materialparametern. Um die geforderte Übereinstimmung von äquivalenter Breite und der Breite des Pfades der Schädigungsrate zu erreichen, werden neue Ansätze für die äquivalente Breite vorgeschlagen, die die erwähnten Einflüsse berücksichtigen. Simulationen unter Verwendung der neuen Ansätze für die äquivalente Breite führen zu einer guten Übereinstimmung von äquivalenter Breite und der Breite des Pfades der Schädigungsrate in der Simulation. Die Ergebnisse der Simulationen, wie z.B. Last-Verformungsbeziehung und Rissverläufe, sind netzunabhängig und stimmen gut mit den experimentellen Beobachtungen überein. Basierend auf den gewonnenen Erkenntnissen wird eine Erweiterung der Bruchenergetischen Regularisierung vorgeschlagen: die Adaptive Bruchenergetische Regularisierung. Im abschließenden Kapitel der Arbeit werden mit der vorgeschlagenen Theorie, dem neuen Schädigungsmodell und der Adaptiven Bruchenergetischen Regularisierung, noch zwei in der Literatur gut dokumentierte Versuche simuliert. Die Simulationsergebnisse entsprechen den experimentellen Beobachtungen. / This doctoral thesis deals with the simulation of predominantly tensile loaded plain concrete structures. Concrete is modeled on the macro level and the Finite Element Method is applied to solve the resulting mechanical field problem. A new material model for concrete based on continuum damage mechanics and an extended regularization technique based on the fracture energy approach are presented. The thesis is subdivided into two parts. In the first part, a local, anisotropic damage model for concrete is derived. This model uses a symmetric second-order tensor as the damage variable, which enables the simulation of orthotropic degradation. The validity of the first and the second law of thermodynamics as well as the validity of the principle of maximum dissipation rate are required. Using a normal rule in the space of the dissipative forces, which are the thermodynamically conjugated variables to the damage variables, and the definition of the loading functions in the space of the dissipative forces guarantee their validity. The suggested formulation contains three decoupled loading functions. A further requirement in the derivation of the model was the minimization of the number of material parameters, which should be determined by a small number of standard experiments. The material parameters of the new damage model are the Young’s modulus, the Poisson’s ratio, the tensile strength and the fracture energy per unit area. The second part of the work focuses on localization and regularization. If a Finite Element simulation is performed using a local material model for concrete, the results of the Finite Element simulation are mesh-dependent. To attain mesh-independent simulations, a regularization technique must be applied. The fracture energy approach, which is characterized by introducing a characteristic length in a locally formulated material model, is used as regularization technique in this work. The choice of a value for the characteristic length is founded by the requirement, that the fracture energy per unit area, which is consumed for the fracture process in the simulation, must be the same as in experiment, i.e. the energy dissipation must be correct. In this dissertation, the new idea is suggested that the correct energy dissipation can be only attained if the characteristic length in the material model coincides in every loading increment with the width of the energy-dissipating zone in the simulation. The energy-dissipating zone in a simulation is formed by the integration points with increasing damage and obtains the name: damage rate path. Detailed investigations based on simulations of a double-edge notched specimen under mixed-mode loading are performed with varying characteristic lengths in order to obtain information concerning the evolution of the damage rate path during a simulation. All simulations were performed using displacement-based elements with four nodes. The range with expected damage was always finer and regularly discretized. The results of the simulations show that the width of the damage rate path depends on the damage at the specific material point, on the angle between damage direction and element edges, on the element size and on the material parameters. Based on these observations, new approaches for the characteristic length are suggested in order to attain the coincidence of the characteristic length with the width of the damage rate path. Simulations by using the new approaches yield a sufficient coincidence of the characteristic length with the width of the damage rate path. The simulations are mesh-independent and the results of the simulation, like load-displacement curves or crack paths, correspond to the experimental results. Based on all new information concerning the regularization technique, an extension of the fracture energy approach is suggested: the adaptive fracture energy approach. The validity and applicability of the suggested theory, the new anisotropic damage model and the adaptive fracture energy approach, are verified in the final chapter of the work with simulations of two additional experiments, which are well documented in the literature. The results of the simulations correspond to the observations in the experiments.

Die vorliegende Arbeit befasst sich mit der kontinuumsmechanischen Formulierung des Aushärteverhaltens polymerer Werkstoffe sowie der Implementierung und Simulation von Aushärtestoffgesetzen im Rahmen der Finite-Elemente-Methode. Auf Basis eines allgemeinen Modellierungsrahmens wird ein konkretisiertes Stoffgesetz für die Nachbildung von Aushärteprozessen eines acrylischen Knochenzements entwickelt. Darüber hinaus werden verschiedene Finite-Elemente-Simulationen zum klinischen Verfahren der Vertebroplastie präsentiert. / This work deals with the continuum mechanical formulation of curing phenomena in polymers as well as the implementation and simulation of curing models within the finite element method. Based on a general modelling framework, a specified material model for the simulation of curing processes in an acrylic bone cement is developed. Moreover, different finite element simulations regarding the clinical procedure of vertebroplasty are presented.

Die vorliegende Arbeit widmet sich der Aufgabe makroskopische Berechnungsmodelle zur Beschreibung des Drahtsägens zu erarbeiten. Ziel ist es, die wesentlichen Effekte abzubilden und den Einfluss von Prozessparametern auf die Dynamik des Systems zu bestimmen. Ein zentraler Punkt ist die Modellierung des bewegten Sägedrahtes. Durch die dem Kontinuum an den Auflagern aufgeprägte Führungsbewegung sind einerseits die Randbedingungen und andererseits ortsfest auf den Draht wirkende Lasten nichtmateriell. Die korrekte kinematische Beschreibung dieses Sachverhaltes ist essentielle Grundlage für die spätere Anwendung des Prinzips von HAMILTON. Durch die Führungsbewegung, die Formulierung der Kontaktkräfte als Folgelasten und durch explizit zeitabhängige Systemparameter ergibt sich ein kompliziertes Systemverhalten. Die dargestellten Berechnungsergebnisse umfassen Studien zu stationären Lagen, die Berechnung von Eigenfrequenzen, Stabilitätsnachweise des dynamischen Grundzustandes, die Bestimmung von Zeitlösungen und die Simulation des Materialabtrages beim Einschnitt. / The aim of the present thesis is to generate macroscopic models to describe the wire sawing process. The principal purpose is to illustrate basic effects and to investigate the influence of important process parameters relating to the dynamics of the system. A fundamental point is the modeling of the moving wire. Because of the axially movement of the continuum the boundary conditions and spatial acting loads are non-material. The precise kinematical description of this issue is the pre-condition for the correct evaluation of HAMILTON’s principle to characterize the dynamics of the system. The resultant complex system behavior is a consequence of the movement of the wire, of the formulation of the contact forces as follower loads and of explicitly time-dependent model parameters. The results of research contain studies of steady state equilibrium solutions and the proof of their LJAPUNOW stability, the calculation of eigenfrequencies, steady state time solutions under harmonically oscillating contact forces and the simulation of the material removal during the cutting process.

Die vorliegende Arbeit beschäftigt sich mit Festkörperkontinuumsmechanik und Metall- bzw. Kristallplastizität auf verschiedenen Längenskalen. Diesbezüglich besteht die Arbeit aus drei größeren Teilen. Im ersten Teil werden Verformungsvorgänge mit expliziter FEM (Finite-Elemente-Methode) und einem makroskopischen phänomenologischen Modell der Viskoplastizität simuliert. Hierbei wird sich auf das Gleichkanalwinkelpressen (ECAP) eines Metallbarrens und die Stauchung einer sogenannten Crashbox konzentriert. In beiden Fällen gelingt es, die im Experiment bereits beobachtete Lokalisierung der Verformung korrekt wiederzugeben. Da bei den Simulationen die konkrete Mikrostruktur des Materials vernachlässigt wird, werden diese Lokalisierungsphänomene als verformungsinduziert angesehen. Der zweite Teil beschäftigt sich mit der Erweiterung des viskoplastischen Modells, sodass mikroskopische Vorgänge der Gitterdefektstruktur des Materials berücksichtigt werden können. Dazu wird ein Modell des dynamischen Verhaltens von Versetzungspopulationen entwickelt und an das makroskopische viskoplastische Modell gekoppelt. Auf diese Weise können Aspekte der sogenannten Kornfeinung – einem komplexen Strukturbildungsprozess von Versetzungen und anderen Gitterdefekten – erfasst werden. Allerdings kann die für die makroskopischen Eigenschaften entscheidende Bildung von Subkorngrenzen auf diese Weise nicht abgebildet werden. Um dies zu erreichen, wird im dritten Teil der Arbeit eine mesoskopische Theorie der Kristallplastizität mit kontinuierlich verteilten Versetzungen verwendet und weiterentwickelt. Hierbei werden die für eine Subkornbildung wesentlichen Freiheitsgrade hinzugenommen, die Anzahl phänomenologischer Ansätze und zugehöriger Materialparameter aber so klein wie möglich gehalten. Mit dieser Kontinuumsversetzungstheorie (KVT) gelingt es, die Bildung von Subkorngrenzen bei großen plastischen Verformungen eines Kristallits zu verfolgen. Bei den impliziten FEM-Simulationen wird ebenfalls eine Lokalisierung beobachtet, allerdings in Bezug auf die Aktivität der Versetzungen in verschiedenen Gleitebenen. Dementsprechend wird dieses Lokalisierungsphänomen als versetzungsinduziert angesehen. Der Beitrag der vorliegenden Arbeit liegt zum einen in der Aufarbeitung und Gegenüberstellung unterschiedlicher methodischer Herangehensweisen zur Modellierung verformungs- und versetzungsinduzierter Lokalisierungsphänomene. Zum anderen wird eine Analyse und Vereinheitlichung der geometrisch linearen KVT nach Berdichevsky & Le vorgenommen. Wie sich dabei zeigt, verhindern inhärente kinematische Einschränkungen der Theorie die Simulation einer Subkornbildung. Aus diesem Grund wird die konsistente geometrisch nichtlineare KVT von Gurtin aufgegriffen und erweitert. Mit einem daraus abgeleiteten elastisch und plastisch anisotropen Modell der Einkristallviskoplastizität wird der Nachweis erbracht, dass die Subkornbildung damit simuliert werden kann. Darüber hinaus wird eine Aufbereitung und Synthese von Algorithmen zur numerischen Lösung der zugehörigen Feldgleichungen mittels der Methode der finiten Differenzen und der finiten Elemente geliefert. Zudem werden beide Näherungsverfahren in Bezug auf Vor- und Nachteile sowie thermodynamische Konsistenz bei der Anwendung auf Mehrfeldprobleme miteinander verglichen. / The present thesis deals with solid continuum mechanics applied to metal and crystal plasticity on different length scales. In this respect, the work consists of three larger parts. In the first part, deformation processes are simulated with explicit FEM (Finite Element Method) and a macroscopic phenomenological model of viscoplasticity. Here the focus is on the Equal-Channel Angular Pressing (ECAP) of a metal billet and the compression of a so-called crash box. In both cases it is possible to correctly reproduce the localization of the deformation as already observed in the experiment. Since the concrete microstructure of the material is neglected in the simulations, these localization phenomena are regarded as deformation-induced. The second part deals with the extension of the viscoplastic model so that microscopic processes of the lattice defect structure of the material can be considered. A model of the dynamic behavior of dislocation populations is developed and coupled to the macroscopic viscoplastic model. In this way, aspects of the so-called grain refinement – a complex structure formation process of dislocations and other lattice defects – can be captured. However, the formation of subgrain boundaries, which is decisive for the macroscopic properties, cannot be predicted in this way. To achieve this, a mesoscopic theory of crystal plasticity with continuously distributed dislocations is used and further developed in the third part of the thesis. Here, the degrees of freedom essential for subgrain formation are added, while the number of phenomenological approaches and associated material parameters are kept as small as possible. With this continuum dislocation theory it is possible to follow the formation of subgrain boundaries during large plastic deformations of a crystallite. In the implicit FEM simulations, localization is also observed, but with respect to the dislocation activity in different slip planes. Accordingly, this localization phenomenon is considered dislocation-induced. The contribution of the present work lies on the one hand in the review and comparison of different methodical approaches to the modeling of deformation- and dislocation-induced localization phenomena. On the other hand, an analysis and unification of the geometrically linear continuum dislocation theory according to Berdichevsky & Le is carried out. As it turns out, inherent kinematic limitations of the theory prevent the simulation of subgrain formation. For this reason the consistent geometrically non-linear continuum dislocation theory from Gurtin is adopted and extended. With the derived model of elastically and plastically anisotropic single crystal viscoplasticity it is proven that subgrain formation can be simulated. Moreover, a preparation and synthesis of algorithms for the numerical solution of the associated field equations using the method of finite differences and finite elements is provided. In addition, both approximation methods are compared in terms of advantages and disadvantages as well as thermodynamic consistency when applied to multi-field problems.