The problem of a thin elastomeric layer confined between two stiff adherends arises in numerous applications such as microelectronics, bio-inspired adhesion and the manufacture of soft biomedical products. A common requirement is that the debonding of the elastomeric layer from the adherends be controlled to avoid undesirable failure modes. This level of control may necessitate understanding the collective role of the interfacial adhesion, material properties, part geometries, and loading conditions on the debonding. Analytical and numerical approaches using the finite element method and a cohesive zone model (CZM) for the interfacial debonding are used in this dissertation to delineate the role of the afore-mentioned parameters on the initiation and propagation of debonding for both rigid and non-rigid adherends. Extensively studied in the dissertation is the debonding of a semi-infinite relatively stiffer adherend from an elastomer layer with its other surface firmly bonded to a rigid base. The adherend is pulled upwards by applying normal displacements either on its entire unbonded surface or on the edge of its part overhanging from the elastomer layer. The adherend and the elastomeric layer materials are assumed to be linear elastic, homogeneous and isotropic and the elastomer is assumed to be incompressible. Viscoelasticity of the elastomer is considered in the first part of the work. Plane strain deformations of the system with a bilinear traction-separation (TS) relation in the CZM are analyzed. Two non-dimensional numbers, one related to the layer confinement and the other to the interfacial TS parameters, are found to determine if debonding initiates at interior points in addition to at corner points on the adherend/elastomer interface, and if adhesion-induced instability is exhibited. This work is extended to axisymmetric problems in which debonding can take place at both interfaces. Motivated by an industrial demolding problem, numerical experiments are conducted to derive insights into preferential debonding at one of the two interfaces, including for curved adherends. Results reported herein should help engineers design an elastomer layer sandwiched between two adherends for achieving desired failure characteristics. / Ph. D.

On s'intéresse à la résolution numérique de problèmes de contact et de fissuration en dynamique. Le problème de contact envisagé est le problème de Signorini avec ou sans frottement de Coulomb. Quant au problème de fissuration, il s'agit d'un modèle de zone cohésive avec trajet de fissuration pré-défini. Ces problèmes se caractérisent par la présence d'une condition aux limites non-régulière et se formulent comme des inéquations variationnelles d'évolution ou des inclusions différentielles. Pour les résoudre numériquement, nous combinons, comme il est courant en dynamique des solides, une discrétisation en espace par éléments finis et des schémas d'intégration en temps (de types différences finies). Pour le problème de contact, nous commençons par comparer les principales méthodes proposées dans la littérature. Nous étudions ensuite plus particulièrement la méthode dite de masse modifiée récemment introduite par H. Khenous, P. Laborde et Y. Renard. Nous en proposons une variante semi-explicite. Par ailleurs, nous prouvons un résultat de convergence des solutions semi-discrètes en espace vers une solution continue dans le cas d'un problème de Signorini sans frottement et d'un matériau viscoélastique. Nous analysons également les methodes semi-discrètes en espace et totalement discrètes dans le cas d'un problème de Signorini avec frottement de Coulomb. Pour le problème de fissuration dynamique, la non-régularité de la condition aux limites rend impossible ou peu robuste l'utilisation de schémas totalement explicites. Nous proposons donc des schémas où cette condition aux limites est traitée de façon implicite. Enfin, nous présentons et analysons des méthodes de lagrangien augmenté pour la résolution numérique du problème de fissuration en statique / The present work deals with the numerical solution of dynamic contact and fracture problems. The contact problem is a Signorini problem with or without Coulomb friction. The fracture problem uses a cohesive zone model with a prescribed crack path. These problems are characterized by a non-regular boundary condition and can be formulated with evolutionary variational inequations or differential inclusions. For the numerical solution, we combine, as usual in solid dynamics, a finite element discretization in space and time-integration schemes. For the contact problem, we begin by comparing the main methods proposed in the literature. We then focus on the so-called modified mass method recently introduced by H. Khenous, P. Laborde et Y. Renard, for which we propose a semi-explicit variant. In addition, we prove a convergence result of the space semi-discrete solutions to a continuous solution in the frictionless viscoelastic case. We also analyze the space semi-discrete and fully discrete problems in the friction Coulomb case. For the dynamic fracture problem, using a fully explicit scheme is impossible or not robust enough. Therefore, we propose time-integration schemes where the boundary condition is treated in an implicit way. Finally, we present and analyze augmented Lagrangian methods for static fracture problems

Zyklische Kohäsivzonenmodelle beschreiben irreversibles Separationsverhalten und Schädigungsakkumulation unter zyklischer Belastung. In der vorliegenden Arbeit wird die Formulierung zyklischer Kohäsivzonenmodelle systematisiert und ihr Potenzial zur Simulation von Ermüdungsvorgängen analysiert. Die Kohäsivspannungs-Separations-Beziehungen werden auf Basis etablierter thermodynamischer Konzepte der Schädigungsmechanik aufgestellt. Zyklische Schädigungsakkumulation wird über die Entwicklungsgleichung der Schädigungsvariablen unter Berücksichtigung einer zustandsabhängigen Dauerfestigkeit beschrieben. Das Kohäsivzonenmodell wird erfolgreich für die Simulation von Werkstoffermüdung angewandt. Numerisch mithilfe der Methode der finiten Elemente erzeugte Rissfortschrittskurven bilden das experimentell beobachtete Ermüdungsrisswachstumsverhalten in allen Bereichen ab. Über Parameterstudien wird der Einfluss der einzelnen Modellparameter ermittelt. Darüber hinaus wird die Anwendung des zyklischen Kohäsivzonenmodells auf die Simulation von Wöhler-Versuchen vorgestellt und der Probengrößeneffekt auf das Ermüdungsverhalten untersucht. Der Zusammenhang zwischen den lokalen Beanspruchungszuständen in der Kohäsivzone und dem vorhergesagten globalen Versagensverhalten wird aufgeklärt. Die gewonnenen Erkenntnisse bilden die Grundlage für ein Konzept zur Identifikation der Kohäsivparameter, das auf der Auswertung von Wöhler- und Rissfortschrittskurven beruht. / Cyclic cohesive zone models describe irreversible separation behaviour and damage accumulation under cyclic loading. In the present thesis, the formulation of cyclic cohesive zone models is systemised and their potential to simulate fatigue processes is analysed. The relation between traction and separation is described based on established thermodynamical concepts of damage mechanics. Cyclic damage accumulation is controlled by a damage evolution equation taking into account a state-dependent endurance limit. The cohesive zone model is applied successfully to the simulation of material fatigue. Fatigue crack growth rate curves, which were obtained numerically by means of the finite element method, reproduce the experimentally observed behaviour in all stages. The influences of the particular parameters of the model are determined by parametric studies. In addition, simulations of uniaxial fatigue tests using the cyclic cohesive zone model are presented. Furthermore, the size effect on the fatigue behaviour is investigated. The relation between the local states within the cohesive zone and the predicted global failure modes is explained. These findings form the foundation for a concept of parameter identification which bases on the evaluation of Wöhler-curves and fatigue crack growth rate curves.

Ce travail concerne des expériences mécaniques, des analyses numériques et des modélisations analytiques de la rupture cohésives (Mode I), vis-à-vis de l’étude du comportement mis en évidence par la courbe de Résistance (Courbe-R) et l’effet d’échelle de structures entaillées en bois massif. Des expériences de fissuration sont combinées à des analyses numériques pour déterminer les propriétés de rupture au moyen d’une procédure appelée Théorie de la Mécanique de la Rupture Linéaire Élastique équivalente (TMRLE), basée sur la complaisance de la structure. La courbe-R, obtenue à partir des expériences, selon une méthode de correction du poids propre, montre l’existence d’un domaine endommagé (Zone de Processus de Rupture) de taille non négligeable se développant en fond de fissure. Dans des conditions de fissuration stationnaire, ce domaine atteint une taille critique, et l’énergie nécessaire pour faire propager la fissure avec ce domaine endommagé (par unité de surface de rupture), reste constante. Le taux de libération de l’énergie de fissuration ainsi attendu, joue un rôle important en Mécanique de la Rupture, car il est possible simuler le comportement quasi-fragile du matériau en combinaison avec les autres propriétés de cohésion. La loi d’effet d’échelle de Bažant, utilisée pour prévoir l’influence de la taille sur la contrainte nominale, est estimée à partir de la réunion de deux comportements asymptotiques réalisés sur de petites tailles (Analyse limite ou RdM) et des grandes tailles. Une procédure analytique est présentée pour évaluer le comportement asymptotique additionnel exhibé par la contrainte nominale dans le régime intermédiaire, de façon plus exacte. Une validation numérique est présentée, et l’information expérimentale vient confirmer ce comportement asymptotique. / This work concerns the mechanical testing, numerical analysis and modelling of cohesive fracture (Mode I) on the purpose to study the Resistance-curve behaviour and the size effect in wooden notched structures. The mechanical testing is combined with the numerical analysis to evaluate fracture properties by means of an equivalent LEFM approach based on the structure compliance. The Resistance-curve being revealed from the experiments, by means of a self-weight compensation method, correction puts into evidence that a non-negligible damaged domain (Fracture Process Zone) is under development in the crack front during the loading process. This being the case, among other fracture parameters issued from the Resistance-curve, the critical (asymptotic) energy release rate is determined, turning possible to use it in combination with other cohesive crack properties in the crack modelling (in Mode I). Thus, for a given geometry it is possible to monitor the critical dimension being revealed by the Fracture Process Zone (FPZ) during the crack propagation. The well known Bažant’s size effect law provides the scaling of the nominal strength through the asymptotic matching performed both on the small (Strength Theory) and on the large (LEFM) structure sizes. An analytical procedure is proposed to determine an additional asymptotic regime in the intermediate size range through a more accurate manner. Numerical validations of the proposed procedure are made and experimental data is presented revealing the scaling of the nominal strength through an envelop of values.

La prédiction de l’état de fissuration est un enjeu crucial pour l’analyse des structures en béton armé, qui nécessite le recours à la modélisation et à la simulation numérique. Le calcul par éléments finis des structures en béton armé pose au moins deux problèmes majeurs :d’une part il existe peu de modèles permettant de traiter à la fois l’initiation, la propagation et l’ouverture des fissures, d’autre part le diamètre généralement faible des armatures métalliques par rapport aux dimensions des structures étudiées nécessite des maillages particulièrement fins. On propose donc des solutions à ces deux problématiques basées sur l’utilisation de fonctions de niveau (level set). L’endommagement et la fissuration du béton sont modélisés à l’aide de l’approche TLS (Thick Level Set). Cette méthode,développée en tant que méthode de régularisation des modèles d’endommagement locaux, utilise une level set afin d’introduire une longueur caractéristique. Cela permet de rendre aisée la localisation de la position des fissures, et donc d’enrichir le champ de déplacement parla méthode des éléments finis étendus (X-FEM) afin de modéliser l’ouverture des macro-fissures. Concernant la modélisation des armatures, une nouvelle approche multidimensionnelle est proposée. Une représentation volumique des armatures par la méthode X-FEM est utilisée dans les zones d’intérêt afin d’obtenir des résultats précis tout en simplifiant la procédure de maillage, tandis qu’une représentation linéique est utilisée dans le reste de la structure afin de réduire le nombre de degrés de liberté du calcul. La méthode de transition développée ici permet d’assurer la cohérence des résultats obtenus / Prediction of cracking is a key point for the analysis ofreinforced concrete structures, which requires the use of Modeling and numerical simulation. The analysis of reinforced concrete structures using the finite element method raises two issues: on one hand, few models areable to deal with the initiation, the propagation and the opening of cracks, on the other hand the diameter of thereinforcements which is usually small compared to the dimensions of the structures necessitates very fine meshes. Some solutions to these two problematics areproposed, based on the use of level set functions.Damage and cracking of concrete are modeled using theThick Level Set (TLS) approach. This method,developped as a mean to regularize local damagemodels, uses a level set to introduce a characteristic length. It makes the location of the cracks easy, whichallows to enrich the displacement field with the eXtendedFinite Element Method (X-FEM) in order to model the macro-cracks opening. Concerning the modeling of thereinforcements, a new multidimensionnal approach isproposed. A volumic representation of the reinforcements with the X-FEM method is used in the zones of interest to get accurate results while simplifying the meshing process, whereas a lineic representation isused elsewhere to decrease the number of degrees of freedom. The developed transition method insures the consistency of the results.

Le présent travail de thèse fait suite à une première étude où une stratégie d’identification des paramètres et formes des lois de zones cohésives a été élaborée pour des matériaux homogènes. L’extension au cas de matériaux présentant des hétérogénéités nécessite d’accéder localement aux champs de contraintes.Ainsi, l’objectif principal de cette étude est de mettre au point une méthode de caractérisation locale des propriétés mécaniques et des contraintes. Cette méthode est basée sur l’erreur en relation de comportement combinée à l’exploitation de la richesse des mesures de champs cinématiques planes et plus particulièrement des champs de déformations, obtenus par dérivation numérique des champs de déplacements. Cette mesure cinématique est réalisée par une technique de corrélation d’images numériques enrichie.La méthode d’identification est basée sur la minimisation itérative d’une norme énergétique faisant intervenir le tenseur élastoplastique sécant. Différentes simulations numériques ont illustré la capacité de la procédure à identifier localement des champs de propriétés hétérogènes et sa robustesse et sa stabilité vis-à-vis du bruit de mesure, du choix du jeu de paramètres d’initialisation de l’algorithme et de la finesse du maillage.Pour finir, des essais plans avec différentes géométries d’éprouvettes ont été effectués et un essai a été mis au point pour obtenir de manière maîtrisée un état initial très hétérogène. Les résultats d’identification élastoplastique multilinéaire ont montré la capacité de la méthode à identifier les lois de comportements locales sur ce matériau hétérogène. / The present work follows a first approach where a strategy for identifying the shape and the parameters of cohesive-zone laws has been developed for homogeneous materials. The extension of this method to heterogeneous material requires the knowledge of the local stress state.The study aims at developing a local characterization method for mechanical properties and stresses. This method is based on the constitutive equation gap principles and relies on the knowledge of mechanical kinematic fields and particularly of the strain fields. These fields are obtained by the numerical differentiation of displacement fields measured by digital image correlation.This identification method is based on the iterative minimization of an energy norm involving the secant elastoplastic tensor. Various numerical simulations were used to illustrate the performance of the procedure for locally identifying heterogeneous property fields, and to characterize its robustness and its stability with respect to noise to the values of the algorithm initialization parameter and to the mesh refinement.Finally, various experimental tests with different specimen geometries were performed and a test has been developed to obtain a controlled heterogeneous initial state. The multilinear elastoplastic identification results showed the ability of the method to identify the local behavior properties on heterogeneous materials.

In the present thesis the crack propagation in the ductile-brittle transition region is studied on two scales with deterministic models. In the macroscopic model the ductile failure is described by a non-local Gurson-model whereas the discrete void microstructure is resolved around the crack tip in the microscopic model. The failure by cleavage is not evaluated by means of a post-processing criterion but is modeled equivalently using a cohesive zone model on both scales. Thus, cleavage is not a priori identified with unstable crack propagation but the transition between stable and unstable mode of propagation is a result of the simulation. The problem of handling completely failed material within the framework of non-local damage models is pointed out. A method to overcome this problem is proposed and successfully applied. The case of contained plastic yielding at the crack tip is addressed with a modified-boundary layer model. The macroscopic simulations reproduce many features which are known from experiments like the formation of stretch zones, cleavage after initial ductile tearing, pop-ins with crack arrest, among others. The microscopic simulations substantiate the understanding of the macroscopically observed behavior. Systematic parameter studies are performed. Starting with considerations on the limit cases like pure ductile failure or the lower-ductile brittle transition region allows to separate the effects of the different constitutive parameters. Based on these results, a methodology is proposed to extract the macroscopic material parameters from experiments. This scheme is successfully applied to experimental data from literature. The results show that the behavior of a low-constraint specimen can be reliably predicted with the parameters extracted from a high-constraint specimen. / In der vorliegenden Arbeit wird die Rissausbreitung im spröd-duktilen Übergangsbereich auf zwei Skalen mittels deterministischer Modelle untersucht. Das duktile Versagen wird im makroskopischen Modell durch ein nichtlokales Gurson-Modell beschrieben, während im mikroskopischen Modell die Porenmikrostruktur im Bereich um die Rissspitze diskret aufgelöst wird. Das mögliche Versagen durch Spaltbruch wird nicht, wie üblich, nachträglich durch ein spannungsbasiertes Kriterium bewertet. Stattdessen wird der Spaltbruch auf beiden Skalen durch ein Kohäsivzonenmodell abgebildet. Somit wird die Spaltbruchinitiierung nicht a priori mit instabiler Rissausbreitung gleichgesetzt. Vielmehr ist die Stabilität der Rissausbreitung ein Ergebnis der Simulationen. Außerdem wird das Problem der der Handhabung vollständig ausgefallenen Materials im Rahmen nichtlokaler Schädigungsmodelle herausgestellt. Es wird eine Methode vorgestellt, dieses Problem zu behandeln und erfolgreich angewendet. In den Simulationen wird der Fall vollständig eingebetteten, plastischen Fließens untersucht. Die Simulationen mit dem makroskopischen Modell geben viele Effekte wieder, die aus Experimenten bekannt sind. Dazu zählen die Ausbildung von Stretchzonen, die Spaltbruchinitiierung nach anfänglichem, duktilem Reißen oder lokale Instabilitäten mit Rissarrest. Die mikroskopischen Simulationen tragen zum Verständnis des makroskopisch beobachteten Verhaltens bei. In der vorliegenden Arbeit werden systematische Parameterstudien durchgeführt. Zunächst werden Grenzfälle wie das rein duktile Versagens oder der Spaltbruch in Abwesenheit der Mikroporen untersucht, um die Einflüsse der einzelnen Materialparameter abzugrenzen. Ausgehend von diesen Ergebnissen wird eine Prozedur vorgeschlagen, die Materialparameter des makroskopischen Modells Schritt für Schritt aus Experimenten zu bestimmen. Diese Prozedur wird erfolgreich auf experimentelle Daten aus der Literatur angewendet. Die Ergebnisse zeigen, dass es das entwickelte Modell erlaubt, das Verhalten einer Bruchmechanikprobe mit geringer Dehnungsbehinderung an der Rissspitze mit denjenigen Materialparametern vorherzusagen, die an Proben mit einer hohen Dehnungsbehinderung ermittelt wurden.

Die Arbeit beschreibt eine Modellierung des Materialverhaltens von Faser-Kunststoff-Verbunden unter Berücksichtigung der lokalen Materialstruktur, den konstitutiven Eigenschaften der Verbundbestandteile sowie charakteristischer Schädigungsphönomene. Die Diskretisierung eines repräsentativen Ausschnitts der Materialstruktur erfolgt unter Verwendung der erweiterten Finiten-Elemente-Methode (XFEM). Sie ermöglicht die effiziente Modellierung des Steifigkeitssprunges an den inneren Materialgrenzen und deren Versagen. Der Verlauf der Elementgrenzen muss dabei nicht an die Materialstruktur angepasst werden. Für die Beschreibung der Dehnratenabhängigkeit der polymeren Matrix wird ein Modell der nichtlinearen fraktionalen Viskoelastizität angewendet. Die Kombination mit einem nichtlokalen Kontinuumsschädigungsmodell ermöglicht weiterhin die Modellierung einer verzerrungsgesteuerten Schädigung des Matrixwerkstoffs. Die Parametrisierung, Validierung des Gesamtmodells erfolgt anhand ausgewählter experimenteller Untersuchungen an einem unidirektional verstärkten Glasfaser-Polypropylen-Verbund. / The thesis addresses the modelling of the material behavior of fibre reinforced polymers. It systematically includes the influence of the local material structure, the mechanical behaviour of the consituents and characteristic damage phenomena. The diskretisation of a representative volume element of the material structure is based on the extended finite element method (XFEM). It allows for an efficient modelling of the stiffness jump at internal material boundaries as well as their damage. With the XFEM, the element boundaries are no longer required to coincide with the material structure. The approximation of the strain rate dependence of the polymeric matrix is based on a nonlinear, fractional viscoelasticity approach. Its combination with a nonlocal strain driven continuum damage modell allows for the modelling of damage effects. The parametrisation and validation of the overall approach is based on a comparison with experimental results for a unidirectional reinforced glass-fibre-polypropylene composite.

Double cantilever beam (DCB) specimens with distinct intermetallic microstructures and different geometries were fractured under different mode ratios of loading, ψ, to obtain critical strain energy release rate, Jc. The strain energy release rate at crack initiation, Jci, increased with phase angle, ψ, but remained unaffected by the joint geometry. However, the steady-state energy release rate, Jcs, increased with the solder layer thickness. Also, both the Jci and Jcs decreased with the thickness of the intermetallic compound layer. Next, mode I and mixed-mode fracture tests were performed on discrete (l=2 mm and l=5 mm) solder joints arranged in a linear array between two copper bars to evaluate the J = Jci (ψ) failure criteria using finite element analysis. Failure loads of both the discrete joints and the joints in commercial electronic assemblies were predicted reasonably well using the Jci from the continuous DCBs. In addition, the mode-I fracture of the discrete joints was simulated with a cohesive zone model which predicted reasonably well not only the fracture loads but also the overall load-displacement behavior of the specimen. Additionally, the Jci calculated from FEA were verified estimated from measured crack opening displacements in both the continuous and discrete joints. Finally, the pad-crater fracture mode of solder joints was characterized in terms of the Jci measured at various mode ratios, ψ. Specimens were prepared from lead-free chip scale package-PCB assemblies and fractured at low and high loading rates in various bending configurations to generate a range of mode ratios. The specimens tested at low loading rates all failed by pad cratering, while the ones tested at higher loading rates fractured in the brittle intermetallic layer of the solder. The Jci of pad cratering increased with the phase angle, ψ, but was independent of surface finish and reflow profile. The generality of the J =Jci(ψ) failure criterion to predict pad cratering fracture was then demonstrated by predicting the fracture loads of single lap-shear specimens made from the same assemblies.

Double cantilever beam (DCB) specimens with distinct intermetallic microstructures and different geometries were fractured under different mode ratios of loading, ψ, to obtain critical strain energy release rate, Jc. The strain energy release rate at crack initiation, Jci, increased with phase angle, ψ, but remained unaffected by the joint geometry. However, the steady-state energy release rate, Jcs, increased with the solder layer thickness. Also, both the Jci and Jcs decreased with the thickness of the intermetallic compound layer. Next, mode I and mixed-mode fracture tests were performed on discrete (l=2 mm and l=5 mm) solder joints arranged in a linear array between two copper bars to evaluate the J = Jci (ψ) failure criteria using finite element analysis. Failure loads of both the discrete joints and the joints in commercial electronic assemblies were predicted reasonably well using the Jci from the continuous DCBs. In addition, the mode-I fracture of the discrete joints was simulated with a cohesive zone model which predicted reasonably well not only the fracture loads but also the overall load-displacement behavior of the specimen. Additionally, the Jci calculated from FEA were verified estimated from measured crack opening displacements in both the continuous and discrete joints. Finally, the pad-crater fracture mode of solder joints was characterized in terms of the Jci measured at various mode ratios, ψ. Specimens were prepared from lead-free chip scale package-PCB assemblies and fractured at low and high loading rates in various bending configurations to generate a range of mode ratios. The specimens tested at low loading rates all failed by pad cratering, while the ones tested at higher loading rates fractured in the brittle intermetallic layer of the solder. The Jci of pad cratering increased with the phase angle, ψ, but was independent of surface finish and reflow profile. The generality of the J =Jci(ψ) failure criterion to predict pad cratering fracture was then demonstrated by predicting the fracture loads of single lap-shear specimens made from the same assemblies.