L’analyse isogéométrique pour la résolution de problèmes de la mécanique du solide suscite de vifs intérêts depuis une dizaine d’année. En effet, cette méthode de discrétisation autorise la description exacte des géométries étudiées permettant ainsi de supprimer les erreurs dues à une mauvaise description du domaine spatial étudié. Cependant elle pose un problème théorique de propagation de raffinement lors de la localisation de maillage. Des méthodes pour contourner ce problème ont été proposée dans la littérature mais complexifient grandement la mise en œuvre de cette stratégie de résolution. Cette thèse propose une stratégie de raffinement localisé adaptatif en espace pour les problèmes de statique et en espace temps pour les problèmes de dynamique transitoire dans le cadre de l’analyse isogéométrique. Pour cela une méthode de localisation pour l’analyse isogéométrique en statique basée sur une résolution multigrille est tout d’abord développée pour des problèmes en deux dimensions. Elle présente l’avantage de contourner la problématique de propagation de raffinement de maillage due à l’analyse isogéométrique tout en étant plus simple à mettre en œuvre que les méthodes déjà existantes. De plus, l’utilisation de l’analyse isogéométrique permet de simplifier les procédures de raffinement lors de l’adaptation de maillage qui peuvent être complexes lors de l’utilisationd’éléments finis classiques. Une méthode de raffinement adaptatif espace temps basée sur une résolution multigrille est ensuite développée pour des problèmes en une dimension. Une étude sur la structure des opérateurs est proposée afin de choisir un intégrateur temporel adapté. Les performances de cette stratégies sont mises en évidence, puis une modification de la méthode de résolution est proposée afin de diminuer significativement les coûts de calculs associées à cette résolution. La méthode de raffinement adaptatif espace temps est appliquée à quelques exemples académiques afin de valider son bon comportement lors de la localisation. / Isogeometric analysis applied to structural mechanics problems is a topic of intense concerns for a decade. Indeed, an exact description of geometries studied is allowed by this discretization method suppressing errors due to a bad description of the spatial domain considered. However, a theoretical problem of refinement propagation appears during mesh localization. Local refinement methods for isogeometric analysis has been developed and implied a complexification of the implementation of such a resolution strategy. This PhD thesis expose a space adaptative refinement strategy for linear elastic problems and a space-time one for transient dynamic using isogeometric analysis. For this purpose, a localization method for isogeometric analysis based on a multigrid resolution is developed for 2D linear elastic problems. This method allow to circumvent mesh refinement propagation inherent to isogeometric analysis, and is easier to implement than existing methods. Moreover, the use of isogeometric analysis simplifies refinement procedures occuring during mesh adaptation and which can be really complex using classical finite element analysis. Then, a space-time adaptative refinement based on a multigrid resolution is developed for one dimensional in space problems. A study on operators structure is exposed in order to choose a well suited time integrator. This strategy's performances are highlighted, then an evolution of this method is set up in order to lower computational costs. The space-time adaptaptive refinement is applied to some academical examples to show it good behavior during localization.

Thesis deals with solving the problems of continuum mechanics by method of Isogeometric analysis. This relatively young method combines the advantages of precise NURBS geometry and robustness of the classical finite element method. The method is described on procedure of solving a plane Poissons boundary value problem. Solver is implemented in MatLab and algorithms are attached to the text.

Dans le contexte industriel actuel, où la simulation numérique joue un rôle majeur, de nombreux outils sont développés afin de rendre les calculs les plus performants et exacts possibles en utilisant les ressources numériques de façon optimale. Parmi ces outils, ceux non-intrusifs, c’est-à-dire ne modifiant pas les codes commerciaux disponibles mais permettant d’utiliser des méthodes de résolution avancées telles que l’analyse isogéométrique ou les couplages multi-échelles, apparaissent parmi les plus attirants pour les industriels. L’objectif de cette thèse est ainsi de coupler l’Analyse IsoGéométrique (AIG) et la Méthode des Éléments Finis (MEF) standard pour l’analyse de détails structuraux par une approche non-intrusive et certifiée. Dans un premier temps, on développe un lien global approché entre les fonctions de Lagrange, classiquement utilisées en éléments finis et les fonctions NURBS bases de l’AIG, ce qui permet d’implémenter des analyses isogéométriques dans un code industriel EF vu comme une boîte noire. Au travers d’exemples linéaires et non-linéaires implémentés dans le code industriel Code_Aster de EDF, nous démontrons l’efficacité de ce pont AIG\MEF et les possibilités d’applications industrielles. Il est aussi démontré que ce lien permet de simplifier l’implémentation du couplage non-intrusif entre un problème global isogéométrique et un problème local éléments finis. Ensuite, le concept de couplage non-intrusif entre les méthodes étant ainsi possible, une stratégie d’adaptation est mise en place afin de certifier ce couplage vis-à-vis d’une quantité d’intérêt. Cette stratégie d’adaptation est basée sur des méthodes d’estimation d’erreur a posteriori. Un estimateur global et des indicateurs d’erreur d’itération, de modèle et de discrétisation permettent de piloter la définition du problème couplé. La méthode des résidus est utilisée pour évaluer ces erreurs dans des cas linéaires, et une extension aux problèmes non-linéaires via le concept d’Erreur en Relation de Comportement (ERC) est proposée. / In the current industrial context where the numerical simulation plays a major role, a large amount of tools are developed in order to perform accurate and effective simulations using as less numerical resources as possible. Among all these tools, the non-intrusive ones which do not modify the existing structure of commercial softwares but allowing the use of advanced solving methods, such as isogeometric analysis or multi-scale coupling, are the more attractive to the industry. The goal of these thesis works is thus the coupling of the Isogeometric Analysis (IGA) with the Finite Element Method (FEM) to analyse structural details with a non-intrusive and certified approach. First, we develop an approximate global link between the Lagrange functions, commonly used in the FEM, and the NURBS functions on which the IGA is based. It’s allowed the implementation of isogeometric analysis in an existing finite element industrial software considering as a black-box. Through linear and nonlinear examples implemented in the industrial software Code_Aster of EDF, we show the efficiency of the IGA\FEM bridge and all the industrial applications that can be made. This link is also a key to simplify the non-intrusive coupling between a global isogeometric problem and a local finite element problem. Then, as the non-intrusive coupling between both methods is possible, an adaptive process is introduced in order to certify this coupling regarding a quantity of interest. This adaptive strategy is based on a posteriori error estimation. A global estimator and indicators of iteration, model and discretization error sources are computed to control the definition of the coupled problem. Residual base methods are performed to estimated errors for linear cases, an extension to the concept of constitutive relation errors is also initiated for non-linear problems.

Industrial development processes as well as research in physics, materials and engineering science rely on computer modelling and simulation techniques today. With increasing computer power, computations are carried out on multiple scales and involve the analysis of coupled problems. In this work, continuum modelling is therefore applied at different scales in order to facilitate a prediction of the effective material or structural behaviour based on the local morphology and the properties of the individual constituents. This provides valueable insight into the structure-property relations which are of interest for any design process. In order to obtain reasonable predictions for the effective behaviour, numerical models which capture the essential fine scale features are required. In this context, the efficient representation of discontinuities as they arise at, e.g. material interfaces or cracks, becomes more important than in purely phenomenological macroscopic approaches. In this work, two different approaches to the modelling of discontinuities are discussed: (i) a sharp interface representation which requires the localisation of interfaces by the mesh topology. Since many interesting macroscopic phenomena are related to the temporal evolution of certain microscopic features, (ii) diffuse interface models which regularise the interface in terms of an additional field variable and therefore avoid topological mesh updates are considered as an alternative. With the two combinations (i) Extended Finite Elemente Method (XFEM) + sharp interface model, and (ii) Isogeometric Analysis (IGA) + diffuse interface model, two fundamentally different approaches to the modelling of discontinuities are investigated in this work. XFEM reduces the continuity of the approximation by introducing suitable enrichment functions according to the discontinuity to be modelled. Instead, diffuse models regularise the interface which in many cases requires even an increased continuity that is provided by the spline-based approximation. To further increase the efficiency of isogeometric discretisations of diffuse interfaces, adaptive mesh refinement and coarsening techniques based on hierarchical splines are presented. The adaptive meshes are found to reduce the number of degrees of freedom required for a certain accuracy of the approximation significantly. Selected discretisation techniques are applied to solve a coupled magneto-mechanical problem for particulate microstructures of Magnetorheological Elastomers (MRE). In combination with a computational homogenisation approach, these microscopic models allow for the prediction of the effective coupled magneto-mechanical response of MRE. Moreover, finite element models of generic MRE microstructures are coupled with a BEM domain that represents the surrounding free space in order to take into account finite sample geometries. The macroscopic behaviour is analysed in terms of actuation stresses, magnetostrictive deformations, and magnetorheological effects. The results obtained for different microstructures and various loadings have been found to be in qualitative agreement with experiments on MRE as well as analytical results. / Industrielle Entwicklungsprozesse und die Forschung in Physik, Material- und Ingenieurwissenschaft greifen in einem immer stärkeren Umfang auf rechnergestützte Modellierungs- und Simulationsverfahren zurück. Die ständig steigende Rechenleistung ermöglicht dabei auch die Analyse mehrskaliger und gekoppelter Probleme. In dieser Arbeit kommt daher ein kontinuumsmechanischer Modellierungsansatz auf verschiedenen Skalen zum Einsatz. Das Ziel der Berechnungen ist dabei die Vorhersage des effektiven Material- bzw. Strukturverhaltens auf der Grundlage der lokalen Werkstoffstruktur und der Eigenschafen der konstitutiven Bestandteile. Derartige Simulationen liefern interessante Aussagen zu den Struktur-Eigenschaftsbeziehungen, deren Verständnis entscheidend für das Material- und Strukturdesign ist. Um aussagekräftige Vorhersagen des effektiven Verhaltens zu erhalten, sind numerische Modelle erforderlich, die wesentliche Eigenschaften der lokalen Materialstruktur abbilden. Dabei kommt der effizienten Modellierung von Diskontinuitäten, beispielsweise Materialgrenzen oder Rissen, eine deutlich größere Bedeutung zu als bei einer makroskopischen Betrachtung. In der vorliegenden Arbeit werden zwei unterschiedliche Modellierungsansätze für Unstetigkeiten diskutiert: (i) eine scharfe Abbildung, die üblicherweise konforme Berechnungsnetze erfordert. Da eine Evolution der Mikrostruktur bei einer derartigen Modellierung eine Topologieänderung bzw. eine aufwendige Neuvernetzung nach sich zieht, werden alternativ (ii) diffuse Modelle, die eine zusätzliche Feldvariable zur Regularisierung der Grenzfläche verwenden, betrachtet. Mit der Kombination von (i) Erweiterter Finite-Elemente-Methode (XFEM) + scharfem Grenzflächenmodell sowie (ii) Isogeometrischer Analyse (IGA) + diffuser Grenzflächenmodellierung werden in der vorliegenden Arbeit zwei fundamental verschiedene Zugänge zur Modellierung von Unstetigkeiten betrachtet. Bei der Diskretisierung mit XFEM wird die Kontinuität der Approximation durch eine Anreicherung der Ansatzfunktionen gemäß der abzubildenden Unstetigkeit reduziert. Demgegenüber erfolgt bei einer diffusen Grenzflächenmodellierung eine Regularisierung. Die dazu erforderliche zusätzliche Feldvariable führt oft zu Feldgleichungen mit partiellen Ableitungen höherer Ordnung und weist in ihrem Verlauf starke Gradienten auf. Die daraus resultierenden Anforderungen an den Ansatz werden durch eine Spline-basierte Approximation erfüllt. Um die Effizienz dieser isogeometrischen Diskretisierung weiter zu erhöhen, werden auf der Grundlage hierarchischer Splines adaptive Verfeinerungs- und Vergröberungstechniken entwickelt. Ausgewählte Diskretisierungsverfahren werden zur mehrskaligen Modellierung des gekoppelten magnetomechanischen Verhaltens von Magnetorheologischen Elastomeren (MRE) angewendet. In Kombination mit numerischen Homogenisierungsverfahren, ermöglichen die Mikrostrukturmodelle eine Vorhersage des effektiven magnetomechanischen Verhaltens von MRE. Außerderm wurden Verfahren zur Kopplung von FE-Modellen der MRE-Mikrostruktur mit einem Randelement-Modell der Umgebung vorgestellt. Mit Hilfe der entwickelten Verfahren kann das Verhalten von MRE in Form von Aktuatorspannungen, magnetostriktiven Deformationen und magnetischen Steifigkeitsänderungen vorhergesagt werden. Im Gegensatz zu zahlreichen anderen Modellierungsansätzen, stimmen die mit den hier vorgestellten Methoden für unterschiedliche Mikrostrukturen erzielten Vorhersagen sowohl mit analytischen als auch experimentellen Ergebnissen überein.

Interface balance laws are derived to describe transport across a phase interface. This is used to derive generalized conditions for phase nucleation and growth, valid even for solids with inhomogeneous surface stress.<div><br></div><div>An explicit interface tracking approach called Enriched Isogeometric Analysis (EIGA) is used to simulate phase evolution. Algebraic level sets are used as a measure of distance and for point projection, both necessary operations in EIGA. Algebraic level sets are observed to often fail for surfaces. Rectification measures are developed to make algebraic level sets more robust and applicable for general surfaces. The proposed methods are demonstrated on electromigration problems. The simulations are validated by modeling electromigration experiments conducted on Cu-TiN line structures.</div><div><br></div><div>To model topological changes, common in phase evolution problems, Boolean operations are performed on the algebraic level sets using R-functions. This is demonstrated on electromigration simulations on solids with multiple voids, and on a bubble coalescence problem. </div>

<p>Moving boundary problems such as solidification, crack propagation, multi-body contact or shape optimal design represent an important class of engineering problems. Common to these problems are one or more moving interfaces or boundaries. One of the main challenges associated with boundary evolution is the difficulty that arises when the topology of the geometry changes. Other geometric issues such as distance to the boundary, projected point on the boundary and intersection between surfaces are also important and need to be efficiently solved. In general, the present thesis is concerned with the geometric arrangement and behavioral analysis of evolving parametric boundaries immersed in a domain. </p> <p>The first problem addressed in this thesis is the packing of ellipsoidal fillers in a regular domain and to estimate their effective physical behavior. Particle packing problem arises when one generates simulated microstructures of particulate composites. Such particulate composites used as thermal interface materials (TIMs) motivates this work. The collision detection and distance calculation between ellipsoids is much more difficult than other regular shapes such as spheres or polyhedra.  While many existing methods address the spherical packing problems, few appear to achieve volume loading exceeding 60%. The packing of ellipsoidal particles is even more difficult than that of spherical particles due to the need to detect contact between the particles. In this thesis, an efficient and robust ultra-packing algorithm termed Modified Drop-Fall-Shake is developed. The algorithm is used to simulate the real mixing process when manufacturing TIMs with hundreds of thousands ellipsoidal particles. The effective thermal conductivity of the particulate system is evaluated using an algorithm based on Random Network Model. </p> <p><br></p> <p>In problems where general free-form parametric surfaces (as opposed to the ellipsoidal fillers) need to be evolved inside a regular domain, the geometric distance from a point in the domain to the boundary is necessary to determine the influence of the moving boundary on the underlying domain approximation. Furthermore, during analysis, since the driving force behind interface evolution depends on locally computed curvatures and normals, it is ideal if the parametric entity is not approximated as piecewise-linear. To address this challenge,  an algebraic procedure is presented here to find the level sets of rational parametric surfaces commonly utilized by commercial CAD systems. The developed technique utilizes the resultant theory to construct implicit forms of parametric Bezier patches, level sets of which are termed algebraic level sets (ALS). Boolean compositions of the algebraic level sets are carried out using the theory of R-functions. The algebraic level sets and their gradients at a given point on the domain can also be used to project the point onto the immersed boundary. Beginning with a first-order algorithm, sequentially refined procedures culminating in a second-order projection algorithm are described for NURBS curves and surfaces. Examples are presented to illustrate the efficiency and robustness of the developed method. More importantly, the method is shown to be robust and able to generate valid solutions even for curves and surfaces with high local curvature or G₀ continuity---problems where the Newton--Raphson method fails due to discontinuity in the projected points or because the numerical iterations fail to converge to a solution, respectively. </p> <p><br></p> <p>Next, ALS is also extended for boundary representation (B-rep) models that are popularly used in CAD systems for modeling solids. B-rep model generally contains multiple NURBS patches due to the trimming feature used to construct such models, and as a result are not ``watertight" or mathematically compatible at patch edges. A time consuming geometry clean-up procedure is needed to preprocess geometry prior to finite element mesh generation using a B-rep model, which can take up to 70% of total analysis time according to literature. To avoid the need to clean up geometry and directly provide link between CAD and CAE integration,  signed algebraic level sets using novel inner/outer bounding box strategy is proposed for point classification of B-rep model. Several geometric examples are demonstrated, showing that this technique naturally models single patch NURBS geometry as well, and can deal with multiple patches involving planar trimming feature and Boolean operation. During the investigation of algebraic level sets, a complex self-intersection problem is also reported, especially for three-dimensional surface. The self-intersection may occur within an interval of interest during implicitization of a curve or surface since the implicitized curve or surface is not trimmed and extends to infinity. Although there is no robust and universal solution the problem, two potential solutions are provided and discussed in this thesis.</p> <p><br></p> <p>In order to improve the computational efficiency of analysis in immersed boundary problems, an efficient local refinement technique for both mesh and quadrature  using the kd-tree data structure is further proposed. The kd-tree sub-division is theoretically proved to be more efficient against traditional quad-/oct-tree subdivision methods. In addition, an efficient local refinement strategy based on signed algebraic level sets is proposed to divide the cells. The efficiency of kd-tree based mesh refinement and adaptive quadrature is later shown through numerical examples comparing with oct-tree subdivision, revealing significant reduction of degrees of freedom and quadrature points.</p> <p><br></p> <p>Towards analysis of moving boundaries problems, an explicit interface tracking method termed enriched isogeometric analysis (EIGA) is adopted in this thesis. EIGA utilizes NURBS shape function for both geometry representation and field approximation. The behavior field is modeled by a weighted blending of the underlying domain approximation and enriching field, allowing high order continuity naturally. Since interface is explicitly represented, EIGA provides direct geometric information such as normals and curvatures. In addition, the blending procedure ensures strong enforced boundary conditions. An important moving boundary problem, namely, reaction-diffusion problem, is investigated using EIGA. In reaction-diffusion problems, the phase interfaces evolve due to chemical reaction and diffusion under multi-physics driven forces, such as mechanical, electrical, thermal, etc. Typical failure phenomenon due to reaction-diffusion problems include void formation and intermetallic compound (IMC) growth. EIGA is applied to study factors and behavior patterns in these failure phenomenon, including void size, current direction, current density, etc. A full joint simulation is also conducted to study the degradation of solder joint under thermal aging and electromigration. </p>

[pt] A Análise Isogeométrica (IGA) é um método de análise numérica de estruturas que surge com a proposta de unificação entre projeto e simulação, permitindo a criação de modelos computacionais que preservam a geometria exata do problema. Essa abordagem é possível por meio de uma classe de funções matemáticas denominadas NURBS (Non-Uniform Rational B-Splines), amplamente utilizadas em sistemas CAD (Computer-Aided Design) para modelagem de curvas e superfícies. Na análise isogeométrica, as mesmas funções que representam a geometria aproximam as variáveis de campo. Neste contexto, foi desenvolvido este trabalho que tem como objetivo fornecer uma ferramenta no âmbito da mecânica computacional para análise isogeométrica bidimensional de problemas de elasticidade linear, incluindo as etapas de modelagem, análise e visualização de resultados. O sistema é composto por dois softwares: FEMEP (Finite Element Method Educational Computer Program), desenvolvido em Python e responsável pela etapa de modelagem geométrica, e FEMOOLab (Finite Element Method Object-Oriented Laboratory), software MATLAB para análise e exibição de resultados. A ferramenta proposta apresenta uma interface gráfica de usuário (GUI) que permite a visualização e manipulação intuitiva de curvas NURBS com recursos avançados de modelagem, como interseção de curvas e recursos de reconhecimento de região que agilizam e simplificam o processo. Uma contribuição significativa deste trabalho reside na capacidade de gerar malhas isogeométricas não estruturadas, utilizando T-Splines baseadas em um algoritmo de decomposição de domínio. O sistema de código aberto permite a colaboração e o desenvolvimento contínuo pela comunidade de usuários e desenvolvedores. / [en] Isogeometric Analysis (IGA) is a numerical analysis method for structures that arises with the proposal of unification between design and simulation, allowing the creation of computational models that preserve the exact geometry of the problem. This approach is possible by a class of mathematical functions called NURBS (Non-Uniform Rational B-Splines), widely used in CAD (Computer-Aided Design) systems for modeling curves and surfaces. In isogeometric analysis, the same functions representing the geometry approximate the field variables. In this context, this work was developed to provide a tool within the scope of computational mechanics for two-dimensional isogeometric analysis of linear elasticity problems, including the steps of modeling, analysis, and visualization of results. The system consists of two software programs: FEMEP (Finite Element Method Educational Computer Program), developed in Python and responsible for the geometric modeling stage, and FEMOOLab (Finite Element Method Object-Oriented Laboratory), a MATLAB software for analysis and display of results. The proposed tool features a graphical user interface (GUI) that allows intuitive visualization and manipulation of NURBS curves with advanced modeling features such as curve intersection and region recognition features that streamline and simplify the process. A significant contribution of this work lies in the ability to generate non-structured isogeometric meshes, using T-Splines based on a domain decomposition algorithm. The open-source system allows collaboration and continuous development by the community of users and developers.

Au cours de la dernière décennie, la méthode d’analyse isogéométrique (AIG) a attiré l’attention des chercheurs grâce à ses capacités supérieures à la méthode standard des éléments finis (MEF). Le concept AIG utilise les mêmes fonctions de base que celles utilisées dans la conception assistée par ordinateur (CAO) pour l’approximation des champs inconnus tels que les déplacements, pression interstitielle ou la température dans la solution des éléments finis d’un problème thermo–hydro–mécanique (éventuellement couplé). Parmi les caractéristiques les plus importantes d’AIG, la régularité, le taux de convergence et surtout sa continuité intrinsèque d’ordre supérieur représentent une nette amélioration par rapport à la méthode standard des éléments finis, permettant d’obtenir des avantages computationnels significatifs en termes de précision de la solution et de efficacité.Ce travail tente d’exploiter les caractéristiques d’AIG pour la résolution numérique des problèmes hydromécaniques (HM) couplés dans les géomatériaux de second gradient de type poro–élastoplastiques partiellement saturés. D’une part, le modèle second gradient appartenant à la théorie des milieux continus avec microstructure assure l’objectivité des résultats en présence de phénomènes de localisation de la déformation en termes d’indépendance de maillage de la solution numérique, ce qui ne peut être réalisé avec des modèles constitutifs classiques qui n’implique pas l’intervention d’une longueur interne. D’autre part, la continuité C1 réalisable au moyen de fonctions de base AIG permet une implémentation directe de tels modèles constitutifs d’ordre supérieur, dans une formulation HM dérivée de l’approche de mélange classique. De plus, la régularité des fonctions de base AIG s’est révélée très efficace dans la modélisation de processus couplés caractérisés par de forts gradients hydrauliques – comme la simulation de la propagation d’un front de saturation dans une pente partiellement saturée. Dernier point, mais non des moindres, il convient de noter que, par rapport aux approches existantes basées sur les multiplicateurs de Lagrange, la méthode AIG pour résoudre les problèmes hydromécaniques (HM) couplés dans les matériaux du second gradient saturé et partiellement saturé permet une réduction considérable du nombre de degrés de libertés requis pour atteindre le même niveau de précision. Cela entraîne non seulement une augmentation significative de l’efficacité de calcul, mais permet également d’étendre la formulation du second gradient à l’analyse de problèmes réalistes en 3D, dont la solution a été présentée pour la première fois dans ce travail.La formulation poro–élastoplastique du second gradient développée dans ce travail est mise en œuvre dans le code orienté vers la recherche GeoPDEs, un code IAG–MEF open source écrit en Matlab et développé à l’Université de Pavia. Sur la base des résultats obtenus dans une large série de problèmes aux limites en 2D et 3D analysées dans ce travail, on peut conclure que la combinaison de AIG et d’élastoplasticité du second gradient représente un outil puissant pour la simulation numérique de problèmes géotechniques caractérisés par de forts couplages multiphysiques, un comportement fortement non linéaire du sol, et des gradients de déplacement et de pression interstitielle fortement localisés. / During the last decade, Isogeometric Analysis (IGA) has drawn the attention of the Finite Element community to its superior capabilities over the standard Finite Element Method (FEM). The IGA concept uses the same basis functions used in Computed Aided Design (CAD) for the approximation of the unknown fields such as displacements, pore pressure or temperature in the Finite Element solution of a (possibly coupled) thermo– hydro–mechanical problem. Among the most relevant features of IGA, its smoothness, its convergence rate and particularly its intrinsic higher–order continuity between elements represent a definite improvement over the standard FEM, which allow to obtain significant computational advantages in terms of accuracy of the solution and computa- tional efficiency.This work attempts to exploit the characteristics of IGA for the numerical solution of coupled hydro–mechanical (HM) problems in saturated and partially saturated second gradient poro–elastoplastic geomaterials. On one hand, the second gradient model belonging to the theory of continua with microstructure ensures the objectivity of the results in presence of strain localization phenomena in terms of mesh independence of the numerical solution, which cannot be achieved with classical constitutive models without an internal length scale. On the other hand, the C1–continuity achievable by means of IGA basis functions allows a straightforward implementation of such higher order constitutive models, within a HM formulation derived from the classical mixture approach. In addition, the smoothness of the IGA basis functions proved to be very efficient in the modeling of coupled processes characterized by strong hydraulic gradients – such as the simulation of the downward propagation of a saturation front in a partially saturated slope subject to rainfall infiltration. Last but not least, it is worth noting that, as compared to the existing approaches based on Lagrange multipliers, the IGA approach to the solution of coupled hydro-mechanical (HM) problems in saturated and partially saturated second gradient materials allows a dramatic reduction in the number of degrees of freedoms required to achieve the same level of accuracy. This not only results in a significant increase of the computational efficiency, but also allows to extend the complete second gradient formulation to the analysis of realistic 3D problems, the solution of which has been presented in this work for the first time.The local second gradient poro–elastoplastic formulation developed in this work is implemented in the research-oriented code GeoPDEs, a Matlab open source IGA–FEM code developed at the University of Pavia. Based on the results obtained in a large series of representative 2D and 3D initial–boundary value problems analyzed in this work, it can be concluded that the combination of IGA and the second gradient elastoplasticity represents a powerful tool for the numerical simulation of geotechnical problems characterized by strong multiphysics couplings, highly nonlinear behavior of the soil, and strongly localized displacement and pore pressure gradients.

L’Analyse isogéométrique (AIG) est une méthode innovante de résolution numérique des équations différentielles, proposée à l’origine par Thomas Hughes, Austin Cottrell et Yuri Bazilevs en 2005. Cette technique de discrétisation est une généralisation de l’analyse par éléments finis classiques (AEF), conçue pour intégrer la conception assistée par ordinateur (CAO), afin de combler l’écart entre la description géométrique et l’analyse des problèmes d’ingénierie. Ceci est réalisé en utilisant des B-splines ou des B-splines rationnelles non uniformes (NURBS), pour la description des géométries ainsi que pour la représentation de champs de solutions inconnus.L’objet de cette thèse est d’étudier la méthode isogéométrique dans le contexte des problèmes hyperboliques en utilisant les fonctions B-splines comme fonctions de base. Nous proposons également une méthode combinant l’AIG avec la méthode de Galerkin discontinue (GD) pour résoudre les problèmes hyperboliques. Plus précisément, la méthodologie de GD est adoptée à travers les interfaces de patches, tandis que l’AIG traditionnelle est utilisée dans chaque patch. Notre méthode tire parti de la méthode de l’AIG et la méthode de GD.Les résultats numériques sont présentés jusqu’à l’ordre polynomial p= 4 à la fois pour une méthode deGalerkin continue et discontinue. Ces résultats numériques sont comparés pour un ensemble de problèmes de complexité croissante en 1D et 2D. / Isogeometric Analysis (IGA) is a modern strategy for numerical solution of partial differential equations, originally proposed by Thomas Hughes, Austin Cottrell and Yuri Bazilevs in 2005. This discretization technique is a generalization of classical finite element analysis (FEA), designed to integrate Computer Aided Design (CAD) and FEA, to close the gap between the geometrical description and the analysis of engineering problems. This is achieved by using B-splines or non-uniform rational B-splines (NURBS), for the description of geometries as well as for the representation of unknown solution fields.The purpose of this thesis is to study isogeometric methods in the context of hyperbolic problems usingB-splines as basis functions. We also propose a method that combines IGA with the discontinuous Galerkin(DG)method for solving hyperbolic problems. More precisely, DG methodology is adopted across the patchinterfaces, while the traditional IGA is employed within each patch. The proposed method takes advantageof both IGA and the DG method.Numerical results are presented up to polynomial order p= 4 both for a continuous and discontinuousGalerkin method. These numerical results are compared for a range of problems of increasing complexity,in 1D and 2D.

<div>As physical testing does not always yield insight into the mechanistic cause of failures, computational modeling is often used to develop an understanding of the goodness of a design and to shorten the product development time. One common, and widely used analysis technique is the Finite Element Method. A significant difficulty with the finite element method is the effort required to generate an analysis-suitable mesh due to the difference in the mathematical representation of geometry CAD and CAE systems. CAD systems commonly use Non-Uniform Rational B-Splines (NURBS) while the CAE tools rely on the finite element mesh. Efforts to unify CAD and CAE by carrying out analysis directly using NURBS models termed Isogeometric Analysis reduces the gap between CAD and CAE phases of product development. However, several challenges still remain in the field of isogeometric analysis. A critical challenge relates to the output of commercial CAD systems. B-rep CAD models generated by commercial CAD systems contain uncoupled NURBS patches and are therefore not suitable for analysis directly. Existing literature is largely missing methods to smoothly couple NURBS patches. This is the first topic of research in this thesis. Fracture-caused failures are a critical concern for the reliability of engineered structures in general and semiconductor chips in particular. The back-end of the line structures in modern semiconductor chips contain multi-material junctions that are sites of singular stress, and locations where cracks originate during fabrication or testing. Techniques to accurately model the singular stress fields at interfacial corners are relatively limited. This is the second topic addressed in this thesis. Thus, the overall objective of this dissertation is to develop an isogeometric framework for parametric domain decomposition and analysis of singular stresses using enriched isogeometric analysis.</div><div><br></div><div>Geometrically speaking, multi-material junctions, sub-domain interfaces and crack surfaces are lower-dimensional features relative to the two- or three-dimensional domain. The enriched isogeometric analysis described in this research builds enriching approximations directly on the lower-dimensional geometric features that then couple sub-domains or describe cracks. Since the interface or crack geometry is explicitly represented, it is easy to apply boundary conditions in a strong sense and to directly calculate geometric quantities such as normals or curvatures at any point on the geometry. These advantages contrast against those of implicit geometry methods including level set or phase-field methods. In the enriched isogeometric analysis, the base approximations in the domain/subdomains are enriched by the interfacial fields constructed as a function of distance from the interfaces. To circumvent the challenges of measuring distance and point of influence from the interface using iterative operations, algebraic level sets and algebraic point projection are utilized. The developed techniques are implemented as a program in the MATLAB environment named as <i>Hierarchical Design and Analysis Code</i>. The code is carefully designed to ensure simplicity and maintainability, to facilitate geometry creation, pre-processing, analysis and post-processing with optimal efficiency. </div><div><br></div><div>To couple NURBS patches, a parametric stitching strategy that assures arbitrary smoothness across subdomains with non-matching discretization is developed. The key concept used to accomplish the coupling is the insertion of a “parametric stitching” or p-stitching interface between the incompatible patches. In the present work, NURBS is chosen for discretizing the parametric subdomains. The developed procedure though is valid for other representations of subdomains whose basis functions obey partition of unity. The proposed method is validated through patch tests from which near-optimal rate of convergence is demonstrated. Several two- and three-dimensional elastostatic as well as heat conduction numerical examples are presented.</div><div><br></div><div>An enriched field approximation is then developed for characterizing stress singularities at junctions of general multi-material corners including crack tips. Using enriched isogeometric analysis, the developed method explicitly tracks the singular points and interfaces embedded in a non-conforming mesh. Solution convergence to those of linear elastic fracture mechanics is verified through several examples. More importantly, the proposed method enables direct extraction of generalized stress intensity factors upon solution of the problems without the need to use <i>a posteriori</i> path-independent integral such as the J-integral. Next, the analysis of crack initiation and propagation is carried out using the alternative concept of configurational force. The configurational force is first shown to result from a configurational optimization problem, which yields a configurational derivative as a necessary condition. For specific velocities imposed on the heterogeneities corresponding to translation, rotation or scaling, the configurational derivative is shown to yield the configurational force. The use of configurational force to analyze crack propagation is demonstrated through examples.</div><div><br></div><div>The developed methods are lastly applied to investigate the risk of ratcheting-induced fracture in the back end of line structure during thermal cycle test of a epoxy molded microelectronic package. The first principal stress and the opening mode stress intensity factor are proposed as the failure descriptors. A finite element analysis sub-modeling and load decomposition procedure is proposed to study the accumulation of plastic deformation in the metal line and to identify the critical loading mode. Enriched isogeometric analysis with singular stress enrichment is carried out to identify the interfacial corners most vulnerable to stress concentration and crack initiation. Correlation is made between the failure descriptors and the design parameters of the structure. Crack path from the identified critical corner is predicted using both linear elastic fracture mechanics criterion and configurational force criterion. </div>