Arun Raina

A phase-field model for fracture in biological tissues

This work presents a recently developed phase-field model of fracture equipped with anisotropic crack driving force to model failure phenomena in soft biological tissues at finite deformations. The phase-field models present a promising and innovative approach to thermodynamically consistent modeling of fracture, applicable to both rate-dependent or rate-independent brittle and ductile failure modes. One key advantage of diffusive crack modeling lies in predicting the complex crack topologies where methods with sharp crack discontinuities are known to suffer. The starting point is the derivation of a regularized crack surface functional that

Phase field modeling of ductile fracture at finite strains: A variational gradient-extended plasticity-damage theory

{\varGamma }

A homogenization approach for nonwoven materials based on fiber undulations and reorientation

Γ-converges to a sharp crack topology for vanishing length-scale parameter. A constitutive balance equation of this functional governs the crack phase-field evolution in a modular format in terms of a crack driving state function. This allows flexibility to introduce alternative stress-based failure criteria, e.g., isotropic or anisotropic, whose maximum value from the deformation history drives the irreversible crack phase field. The resulting multi-field problem is solved by a robust operator split scheme that successively updates a history field, the crack phase field and finally the displacement field in a typical time step. For the representative numerical simulations, a hyperelastic anisotropic free energy, typical to incompressible soft biological tissues, is used which degrades with evolving phase field as a result of coupled constitutive setup. A quantitative comparison with experimental data is provided for verification of the proposed methodology.

A phase‐field model for chemo‐mechanical induced fracture in lithium‐ion battery electrode particles

The ever increasing demand of advanced engineered products also pushes the strengths of the materials used to their theoretical limits. It becomes crucially important to understand the behavior of such materials during failure for an efficient and safe design of the product. This thesis aims at the physical-based numerical modeling of complex failure phenomena in engineering materials, categorized into hard matter and soft matter. In Part I of this thesis, a modification of the well established strong discontinuity approach to model failure phenomena in hard matter by extending it to multiple levels is proposed. This is achieved by the resolution of the overall problem into a main boundary value problem and identified sub-domains based on the concepts of domain decomposition. Those sub- domains are subsequently adaptively discretized during run-time and comprise the so- called sub-boundary value problem to be solved simultaneously with the main boundary value problem. To model failure, only the sub-elements of those sub-boundary value problems are treated by the strong discontinuity approach which, depending on their state of stress, may develop cracks and shear bands. A single finite element of the main boundary value problem can therefore simulate the propagation of multiple propagating strong discontinuities specially arising for simulations of crack branching. The solutions of the different sub-boundary value problems are transferred to the main boundary value problem based on concepts of domain decomposition. The applied boundary conditions are also modified to account for the possible multiple jumps in the displacement fields. It is shown through the simulation of solids undergoing dynamic fracture that the modification allows to predict the onset of crack branching without the need for any artificial crack branching criterion. A close agreement with experiments of the simulation results in terms of micro- and macro branching in addition to studying certain key parameters like critical velocity, dynamic stress intensity factor, and the strain energy release rate at branching is found. In Part II of this thesis, failure phenomena in soft matter is modeled for which an advanced homogenization approach to model the highly anisotropic and non-linear stiffening response at finite strains is developed first. The constituent one-dimensional elements are modeled as linear elastic, by experimental justification, which are modified in the lower strain regime to account for the inherent fiber undulations and the associated fiber unfolding phenomena. Reorientation of these fibers is identified as one primary mechanism for the overall macroscopic stiffening which is achieved by a new bijective mapping asymptotically aligning these fibers with the maximum loading direction in the referential orientation space. A rate-independent evolution law for this map is sought by a physically motivated assumption to maintain the overall elastic framework of the proposed formulation. A closed form solution to the new evolution law is also presented which allows faster computation of updating orientations without resorting to numerical integration or storing history variables. The unit vectors upon reorientation in the referential orientation space are then mapped to the spatial orientation space by the macro deformation gradient to compute the macroscopic Kirchhoff stress and the associated spatial elasticity modulus. A direct comparison of the numerical results with the experimental results from the literature is made which demonstrates the predictive capabilities of the proposed formulation. Finally, the finite deformation extended strong discontinuity approach is utilized to simulate boundary value problems of failure in nonwoven felts. The simulation results of failure show a satisfactory agreement with the experimental data from literature. Die standing steigende Nachfrage nach hochentwickelten technischen Produkten treibt die Zuverlassigkeit der verwendeten Materialien an ihre Grenzen. Um ein effizientes und sicheres Design von Produkten zu gewahrleisten ist es daher von entscheidender Bedeutung, das Verhalten von Materie im Versagensfall zu verstehen. Die vorliegende Arbeit behandelt die physikalisch-basierte numerische Modellierung von komplexen Versagensszenarien von Ingenieurswerkstoffen. Diese lassen sich grob in harte und weiche Materie kategorisieren. Im ersten Teil der Arbeit wird eine Modifikation der wohlbegrundeten Methode der starken Diskontinuitaten vorgenommen, die die Modellierung des Versagens von harten Materie auf meheren Ebenen erlaubt. Die Methode beruht auf einem Gebietzerlegungsansatz, der das zu losende Problem aufteilt in ein ubergeordnetes (Haupt- )Randwertproblem (RWP) sowie weitere identifizierte Sub-Gebiete. Diese Sub-Gebiete werden wahrend der Laufzeit adaptiv diskretisiert und stellen somit ein Sub-Randwertproblem dar, das simultan mit dem ubergeordneten RWP gelost wird. Fur die Simulation des Versagensverhaltens werden ausschlieslich die Sub-Probleme herangezogen, die abhangig vom Spannungszustand Risse und Scherbander ausbilden konnen. Ein einziges finites Element des ubergeordneten RWPs kann somit den Fortschritt mehrerer starker Diskontinuitaten abbilden, was insbesondere bei der Simulation von Rissverzeigungen eine wichtige Rolle spielt. Basierend auf dem Konzept der Gebietszerlegung werden die Losungen der verschiedenen Sub-RWPs an das ubergeordnete RWP transferiert. Die aufgebrachten Randbedingungen werden ebenfalls modifiziert um moglichen Sprungen des Verschiebungsfeldes Rechnung zu tragen. Es werden Simulationen von Materialien mit dynamischem Versagen prasentiert, die auzeigen, dass diese Modifikation dazu in der Lage ist, Rissverzweigungen ohne die zusatzliche Auswertung eines Rissverzweigungskriteriums abzubilden. Eine starke Ubereinstimmung mit experimentellen Daten im Hinblick auf mikroskopische und makroskopische Rissverzweigungseffekte kann beobachtet werden. Gleiches gilt fur wichtige Schlusselparameter wie etwa die kritische Geschwindigkeit, den dynamischen Spannungsintensitatsfaktor sowie die Energiefreisetzungsrate. Im zweiten Teil der Arbeit wird die Modellierung von Versagensszenarien in weichen Materie behandelt. Dieser basiert auf einem Homogenisierungsansatz, der die ausgepragt anisotrope und nichtlinear versteifende Materialantwort bei finiten Deformationen berucksichtigt. Motiviert durch experimentelle Beobachtungen werden eingebettete eindimensionale Elemente als linear elastisches Medium modelliert. Fur den Anfangsbereich der Lastaufbringung wird eine Modifikation vorgenommen, die der inharenten Wellenform der Fasern und den dazu assoziierten Entfaltungsvorgangen Rechnung tragt. Es wird gezeigt, dass die Reorientierung der Fasern einen wesentlichen Mechanismus der effektiven Versteifung des Materials darstellt. Die Integration ebenjenes Mechanismus in das Modell erfolgt in Form einer neuen bijektiven Abbildungsvorschrift, die die Fasern im referentiellen Orientierungsraum asymptotisch in Belastungsrichtung ausrichtet. Es wird eine physikalisch motivierte ratenunabhangige Evolutionsgleichung dieser Abbildungsvorschrift angestrebt, die die zugrundeliegende elastische Struktur der Formulierung erhalt. In diesem Zusammenhang wird eine geschlossene Losung einer neuen Evolutionsgleichung bereitgestellt, die eine schnelle Aufdatierung der Orientierungen ohne Zuhilfenahme von numerischen Integrationsalgorithmen bzw. ohne Vorhalten von geschichtsabhangigen Variablen ermoglicht. Die referentiellen Einheitsvektoren der Orientierungen werden mittels des Deformationsgradienten auf die aktuelle Konfiguration abgebildet und dienen sodann der Berechnung der makroskopischen Kirchhoff-Spannungen sowie dem dazu assoziierten raumlichen Elastizitatsmodul. Der Vergleich mit experimentellen Daten aus der Literatur unterstreicht die pradiktiven Fahigkeiten der vorgeschlagenen Formulierung. Schlieslich wird die hier auf den Bereich finiter Deformationen ausgeweitete Methode der starken Diskontinuitaten auf die Simulation des Versagens von Filz angewandt. Auch hier zeigen die Simuationsergebnisse eine zufriedenstellende Ubereinstimmung mit Daten aus der Literatur.