Les chaussées subissent des sollicitations liées au trafic et au climat conduisant à leur dégradation, par fissuration notamment. Il est nécessaire dans le contexte actuel de pouvoir modéliser le comportement de ces structures multicouches endommagées afin de prévoir leur durée de vie résiduelle ou dimensionner des solutions de renforcement. L’objectif de la thèse est ainsi de proposer un outil de calcul dédié à l'analyse 3D des chaussées fissurées ou comportant des discontinuités. L’approche retenue repose sur la modélisation simplifiée d’une chaussée par un empilement de plaques du Modèle Multi-particulaire des Matériaux Multicouches à 5n équations d’équilibre (M4-5n). Un outil de calcul rapide de référence de chaussées 2D fissurées et une méthode de résolution générale du M4-5n par Eléments Finis mixtes sont développés. Le point de départ de la méthode de résolution est l’écriture, pour le M4-5n, du principe variationnel basé sur le théorème de l'énergie complémentaire où la condition de contraintes statiquement admissibles est assurée à partir de multiplicateurs de Lagrange. La discrétisation des efforts généralisés utilise des espaces d’interpolation permettant le bon conditionnement du système d’équations algébriques à résoudre et garantissant la stabilité de la solution. La méthode est implémentée dans FreeFem++. Elle ramène le problème 3D initial à une modélisation EF 2D et conduit à des valeurs finies des efforts généralisés au niveau des fissures ou décollement d’interface. L’outil de calcul final ainsi développé est validé et appliqué à l’étude de la réponse d’une structure fissurée,représentative d’une chaussée testée en vraie grandeur sur le site de l’IFSTTAR. / Pavements are multilayer structures which undergo cracking distress due to traffic and climatic factors. It is important nowadays to be able to model the mechanical response of such damaged pavements in order to assess their residual lifetime or to design reinforcement solutions. In this context, the present thesis aims at developing a numerical tool dedicated to the analysis of pavements incorporating cracks or discontinuities. In the developed approach, the pavement structure is modeled as a stacking of “plate” elements of typeM4-5n (Multi-Particle Models of Multilayer Materials) which considers 5n equilibrium equations. A reference quick 2D calculation tool for cracked pavements and a general solving of M4-5n by the mixed Finite Element (FE) method was developed. The starting point for this method is the derivation for M4-5n of the variational principle based on the complementary energy theorem whose condition of statically admissible stress is taken into account using Lagrange multipliers. Discretization of the generalized stresses involves interpolation spaces, proposed to avoid ill-conditioned system of algebraic equations after discretization and to insure stability of the solution. The developed method is implemented in a FreeFem++ script. In this method, the initial 3D problem can be handled through FE simulations in 2D and finite values of the generalized stresses are obtained at crack and interlayer debonding locations. The developed numerical tool was validated and applied to the study of the mechanical response of a structure with cracks representative of a pavement tested underfull-scale conditions during an accelerated fatigue test performed at IFSTTAR.

Through enrichment of the elastic potential by the second-order gradient of deformation, gradient elasticity formulations are capable of taking nonlocal effects into account. Moreover, geometry-induced singularities, which may appear when using classical elasticity formulations, disappear due to the higher regularity of the solution. In this contribution, a mixed finite element discretization for finite strain gradient elasticity is investigated, in which instead of the displacements, the first-order gradient of the displacements is the solution variable. Thus, the C1 continuity condition of displacement-based finite elements for gradient elasticity is relaxed to C0. Contrary to existing mixed approaches, the proposed approach incorporates a rot-free constraint, through which the displacements are decoupled from the problem. This has the advantage of a reduction of the number of solution variables. Furthermore, the fulfillment of mathematical stability conditions is shown for the corresponding small strain setting. Numerical examples verify convergence in two and three dimensions and reveal a reduced computing cost compared to competitive formulations. Additionally, the gradient elasticity features of avoiding singularities and modeling size effects are demonstrated.

We connect the Mimetic Finite Difference method (MFD) with the finite-volume two-point flux scheme (TPFA) for Voronoi meshes. The main effect is reducing the saddle-point system to a much smaller symmetric-positive definite matrix. In addition, the generalization allows MFD to seamlessly integrate with existing porous media modeling technology. The generalization also imparts the monotonicity property of the TPFA method on MFD. The connection is achieved by altering the consistency condition of the velocity bilinear operator. First-order convergence theory is presented as well as numerical results that support the claims. We demonstrate a methodology for using MFD in modeling fluid flow in fractures coupled with a reservoir. The method can be used for nonplanar fractures. We use the method to demonstrate the effects of fracture curvature on single-phase and multi-phase flows. Standard benchmarks are used to demonstrate the accuracy of the method. The approach is coupled with existing reservoir simulation technology. / text

Nous développons une nouvelle approche éléments finis pour des équations aux dérivées partielles elliptiques de type élasticité linéaire ou Stokes sur une surface fermée de R3. La surface considérée est décrite par le zéro d'une fonction de niveau assez régulière. Le problème se ramène à la minimisation d'une fonctionnelle énergie pour le champ de vitesse sous contraintes. Les contraintes sont de deux types : (i) la vitesse est tangentielle à la surface, (ii) la surface est inextensible. Cette deuxième contrainte équivaut à l'incompressibilité surfacique du champ de vitesse. Nous abordons ce problème de deux façons : la pénalisation et l'introduction de deux multiplicateurs de Lagrange. Cette dernière méthode a l'avantage de traiter le cas de la limite incompressible d'un écoulement en surface dont nous présentons pour la première fois l'analyse théorique et numérique. Nous montrons des estimations d'erreurs sur la solution discrète et les tests numériques confirment l'optimalité des ces estimations. Pour cela, nous proposons plusieurs approches pour le calcul numérique de la normale et la courbure de la surface. L'implémentation utilise la librairie libre d'éléments finis Rheolef. Nous présentons aussi des résultats de simulations numériques pour une application en biologie : la morphogenèse de l'embryon de la drosophile, durant laquelle des déformations tangentielles d'une monocouche de cellules avec une faible variation d'aire. Ce phénomène est connu sous le nom de l'extension de la bande germinale. / We develop a novel finite element approach for linear elasticity or Stokes-type PDEs set on a closed surface of $mathbb{R}^3$. The surface we consider is described as the zero of a sufficiently smooth level-set function. The problem can be written as the minimisation of an energy function over a constrained velocity field. Constraints areof two different types: (i) the velocity field is tangential to the surface, (ii) the surface is inextensible. This second constraint is equivalent to surface incompressibility of the velocity field. We address thisproblem in two different ways: a penalty method and a mixed method involving two Lagrange multipliers. This latter method allows us to solve the limiting case of incompressible surface flow, for which we present a novel theoretical and numerical analysis. Error estimates for the discrete solution are given andnumerical tests shows the optimality of the estimates. For this purpose, several approaches for the numerical computation of the normal and curvature of the surface are proposed. The implementation relies on the Rheolef open-source finite element library. We present numerical simulations for a biological application: the morphogenesis of Drosophila embryos, duringwhich tangential flows of a cell monolayer take place with a low surface-area variation. This phenomenon is known as germ-band extension.

We address two computational challenges for numerical simulations of Darcy flow problems: rough and uncertain data. The rapidly varying and possibly high contrast permeability coefficient for the pressure equation in Darcy flow problems generally leads to irregular solutions, which in turn make standard solution techniques perform poorly. We study methods for numerical homogenization based on localized computations. Regarding the challenge of uncertain data, we consider the problem of forward propagation of uncertainty through a numerical model. More specifically, we consider methods for estimating the failure probability, or a point estimate of the cumulative distribution function (cdf) of a scalar output from the model. The issue of rough coefficients is discussed in Papers I–III by analyzing three aspects of the localized orthogonal decomposition (LOD) method. In Paper I, we define an interpolation operator that makes the localization error independent of the contrast of the coefficient. The conditions for its applicability are studied. In Paper II, we consider time-dependent coefficients and derive computable error indicators that are used to adaptively update the multiscale space. In Paper III, we derive a priori error bounds for the LOD method based on the Raviart–Thomas finite element. The topic of uncertain data is discussed in Papers IV–VI. The main contribution is the selective refinement algorithm, proposed in Paper IV for estimating quantiles, and further developed in Paper V for point evaluation of the cdf. Selective refinement makes use of a hierarchy of numerical approximations of the model and exploits computable error bounds for the random model output to reduce the cost complexity. It is applied in combination with Monte Carlo and multilevel Monte Carlo methods to reduce the overall cost. In Paper VI we quantify the gains from applying selective refinement to a two-phase Darcy flow problem.

We discuss a full-waveform based material profile reconstruction in two-dimensional heterogeneous semi-infinite domains. In particular, we try to image the spatial variation of shear moduli/wave velocities, directly in the time-domain, from scant surficial measurements of the domain's response to prescribed dynamic excitation. In addition, in one-dimensional media, we try to image the spatial variability of elastic and attenuation properties simultaneously. To deal with the semi-infinite extent of the physical domains, we introduce truncation boundaries, and adopt perfectly-matched-layers (PMLs) as the boundary wave absorbers. Within this framework we develop a new mixed displacement-stress (or stress memory) finite element formulation based on unsplit-field PMLs for transient scalar wave simulations in heterogeneous semi-infinite domains. We use, as is typically done, complex-coordinate stretching transformations in the frequency-domain, and recover the governing PDEs in the time-domain through the inverse Fourier transform. Upon spatial discretization, the resulting equations lead to a mixed semi-discrete form, where both displacements and stresses (or stress histories/memories) are treated as independent unknowns. We propose approximant pairs, which numerically, are shown to be stable. The resulting mixed finite element scheme is relatively simple and straightforward to implement, when compared against split-field PML techniques. It also bypasses the need for complicated time integration schemes that arise when recent displacement-based formulations are used. We report numerical results for 1D and 2D scalar wave propagation in semi-infinite domains truncated by PMLs. We also conduct parametric studies and report on the effect the various PML parameter choices have on the simulation error. To tackle the inversion, we adopt a PDE-constrained optimization approach, that formally leads to a classic KKT (Karush-Kuhn-Tucker) system comprising an initial-value state, a final-value adjoint, and a time-invariant control problem. We iteratively update the velocity profile by solving the KKT system via a reduced space approach. To narrow the feasibility space and alleviate the inherent solution multiplicity of the inverse problem, Tikhonov and Total Variation (TV) regularization schemes are used, endowed with a regularization factor continuation algorithm. We use a source frequency continuation scheme to make successive iterates remain within the basin of attraction of the global minimum. We also limit the total observation time to optimally account for the domain's heterogeneity during inversion iterations. We report on both one- and two-dimensional examples, including the Marmousi benchmark problem, that lead efficiently to the reconstruction of heterogeneous profiles involving both horizontal and inclined layers, as well as of inclusions within layered systems. / text

In der vorliegenden Arbeit wird ein gemischtes Element gegen Locking-Effekte untersucht. Dazu wird ein Fünf-Feld-Hu-Washizu-Funktional (CoFEM-Element) für lineare und quadratische Hexaeder-Elemente unter einer hyperelastischen, isotropen, polykonvexen sowie einer transversal-isotropen Materialformulierung implementiert. Die resultierenden nichtlinearen Gleichungen werden mithilfe eines Mehrebenen-NEWTON-RAPHSON-Verfahren unter Beachtung einer konsistenten Linearisierung gelöst. Als repräsentatives Beispiel der numerischen Untersuchungen dient der einseitig eingespannte Cook-Balken mit einer quadratischen Druckverteilung am Rand. Zur Beurteilung des CoFEM-Elements wird das räumliche Konvergenzverhalten für unterschiedliche Polynomgrade und für verschiedene Netze unter Beachtung der algorithmischen Effizienz untersucht. / This paper presents a mixed finite element formulation of Hu-Washizu type (CoFEM) designed to reduce locking effects with respect to a linear and quadratic approximation in space. We consider a hyperelastic, isotropic, polyconvex material formulation as well as transverse isotropy. The resulting nonlinear algebraic equations are solved with a multilevel NEWTON-RAPHSON method. As a numerical example serves a cook-like cantilever beam with a quadratic distribution of in-plane load on the Neumann boundary. We analyze the spatial convergence with respect to the polynomial degree of the underlying Lagrange polynomials and with respect to the level of mesh refinement in terms of algorithmic efficiency.

In der vorliegenden Arbeit wird ein gemischtes Element gegen Locking-Effekte untersucht. Dazu wird ein Fünf-Feld-Hu-Washizu-Funktional (CoFEM-Element) für lineare und quadratische Hexaeder-Elemente unter einer hyperelastischen, isotropen, polykonvexen sowie einer transversal-isotropen Materialformulierung implementiert. Die resultierenden nichtlinearen Gleichungen werden mithilfe eines Mehrebenen-NEWTON-RAPHSON-Verfahren unter Beachtung einer konsistenten Linearisierung gelöst. Als repräsentatives Beispiel der numerischen Untersuchungen dient der einseitig eingespannte Cook-Balken mit einer quadratischen Druckverteilung am Rand. Zur Beurteilung des CoFEM-Elements wird das räumliche Konvergenzverhalten für unterschiedliche Polynomgrade und für verschiedene Netze unter Beachtung der algorithmischen Effizienz untersucht. / This paper presents a mixed finite element formulation of Hu-Washizu type (CoFEM) designed to reduce locking effects with respect to a linear and quadratic approximation in space. We consider a hyperelastic, isotropic, polyconvex material formulation as well as transverse isotropy. The resulting nonlinear algebraic equations are solved with a multilevel NEWTON-RAPHSON method. As a numerical example serves a cook-like cantilever beam with a quadratic distribution of in-plane load on the Neumann boundary. We analyze the spatial convergence with respect to the polynomial degree of the underlying Lagrange polynomials and with respect to the level of mesh refinement in terms of algorithmic efficiency.

Cette étude a pour objectif principal d’établir une stratégie de modélisation robuste, fiable et performante pour décrire des propagations de fissures d’échelle centimétrique en régime ductile dans des composants industriels. Le modèle d’endommagement de GTN écrit en grandes déformations est utilisé pour modéliser l’endommagement ductile. Ce modèle conduit généralement à une localisation de la déformation, conformément à l’expérience. L’échelle caractéristique de ce phénomène est introduite dans les équations de comportement via l’adoption d’une formulation non locale.Sur le plan numérique, ce modèle non local rend bien compte de la localisation dans une bande d’épaisseur donnée lorsqu’on raffine suffisamment le maillage. Par ailleurs, le problème de verrouillage numérique associé au caractère initialement isochore de la déformation plastique est limité en utilisant une formulation à base d’éléments finis mixtes. Enfin, la distorsion des éléments totalement cassés (i.e. sans rigidité apparente), qui pourrait nuire à la bonne convergence des simulations numériques, est traitée par une régularisation viscoélastique.L’ensemble de ces ingrédients sont appliqués pour simuler la propagation de fissure dans un milieu infini plasticité confinée), de sorte à établir un lien avec les approches globales en J-Δa. L’émoussement, l’amorçage et la (grande) propagation de fissure sont bien prédits. Le modèle est également appliqué à une tuyauterie métallique testée en grandeur réelle dans le cadre du projet européen Atlas+. Après une phase d’identification des paramètres sur éprouvette, les réponses globales et locales d’autres éprouvettes et du tube sont confrontés aux résultats expérimentaux. Ces résultats illustrent le degré de robustesse, de fiabilité et de performance qu’on peut attendre du modèle. / The major goal of this work is to establish a robust, reliable and efficient modeling technique so as to describe ductile tearing over a distance of several centimeters in industrial cases. The GTN damage model expressed in the context of finite strains is chosen to model ductile damage. Generally, the model leads to strain localization in agreement with experimental observations. The characteristic length scale of this phenomenon is introduced into the constitutive equations through the use of a nonlocal formulation.On a numerical ground, the nonlocal model controls the width of the localization band as soon as the mesh is sufficiently refined. Besides, the issue of volumetric-locking associated with plastic incompressibility is handled using a mixed finite element formulation. Finally, the distortion of broken elements (i.e. without any stiffness), which may affect the computational convergence of numerical simulations, is treated using a viscoelastic regularization.The improved GTN model is applied to simulate crack propagation under small-scale yielding conditions, so as to establish a relation with the global (J-Δa) approach. Crack tip blunting, crack initiation and (large) crack propagation are well captured. The model is also applied to a full-scale metallic pipe in the framework of the UE project Atlas+. After a phase of parameter calibration based on the experimental results on some small specimens, the global and local responses of other small specimens and of the full-scale pre-cracked pipe are compared with the experimental results. The results illustrates the robustness, the reliability and the efficiency of the current model.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Neste trabalho é apresentado um novo método acurado com passo de tempo fracionário, baseado em uma técnica de decomposição de operadores, para a solução numérica de um sistema governante de equações diferenciais parciais que modela escoamento trifásico água-gás-óleo imiscível em reservatórios de petróleo heterogêneos no qual os efeitos de compressibilidade do gás não foram levados em conta. A técnica de decomposição de operadores em dois níveis permite o uso de passos de tempo distintos para os três problemas definidos pelo procedimento de decomposição: convecção, difusão e pressão-velocidade. Um sistema hiperbólico de leis de conservação que modela o transporte convectivo das fases fluidas é aproximado por um esquema central de diferenças finitas explícito, conservativo, não oscilatório e de segunda ordem. Este esquema é combinado com elementos finitos mistos, localmente conservativos, para a aproximação numérica dos sistemas de equações parabólico e elíptico associados aos problemas de transporte difusivo e de pressão-velocidade, respectivamente. O operador temporal associado ao sistema parabólico é resolvido fazendo-se uso de uma estratégia implícita de solução (Backward Euler). O modelo matemático para escoamento trifásico considerado neste trabalho leva em conta as forças de capilaridade e expressões gerais para as funções de permeabilidade relativa, campos variáveis de porosidade e de permeabilidade e os efeitos da gravidade. A escolha de expressões gerais para as funções de permeabilidade relativa pode levar à perda de hiperbolicidade escrita e, desta maneira, à existência de uma região elíptica ou de pontos umbílicos para o sistema não linear de leis de conservação hiperbólicas que descreve o transporte convectivo das fases fluidas. Como consequência, a perda de hiperbolicidade pode levar à existência de choques não clássicos (também chamados de choques transicionais ou choques subcompressivos) nas soluções de escoamentos trifásicos. O novo procedimento numérico foi usado para investigar a existência e a estabilidade de choques não clássicos, com respeito ao fenômeno de fingering viscoso, em problemas de escoamentos trifásicos bidimensionais em reservatórios heterogêneos, estendendo deste modo resultados disponíveis na literatura para problemas de escoamentos trifásicos unidimensionais. Experimentos numéricos, incluindo o estudo de estratégias de injeção alternada de água e gás (Water-Alternating-Gas (WAG)), indicam que o novo procedimento numérico proposto conduz com eficiência computacional a resultados numéricos com precisão. Perspectivas para trabalhos de pesquisa futuros são também discutidas, tomando como base os desenvolvimentos reportados nesta tese. / We present a new, accurate fractional time-step method based on an operator splitting technique for the numerical solution of a system of partial differential equations modeling three-phase immiscible water-gas-oil flow problems in heterogeneous petroleum reservoirs in which the compressibility effects of the gas was not take into account. A two-level operator splitting technique allows for the use of distinct time steps for the three problems defined by the splitting procedure: convection, diffusion and pressure-velocity. A system of hyperbolic conservation laws modelling the convective transport of the fluid phases is approximated by a high resolution, nonoscillatory, second-order, conservative central difference scheme in the convection step. This scheme is combined with locally conservative mixed finite elements for the numerical solution of the parabolic and elliptic problems associated with the diffusive transport of fluid phases and the pressure-velocity problem, respectively. The time discretization of the parabolic problem is performed by means of the implicit backward Euler method. The mathematical model for the three-phase flow considered in this work takes into account capillary forces and general expressions for the relative permeability functions, variable porosity and permeability fields, and the effect of gravity. The choice of general expressions for the relative permeability functions may lead to the loss of strict hyperbolicity and, therefore, to the existence of an elliptic region of umbilic points for the systems of nonlinear hyperbolic conservation laws describing the convective transport of the fluid phases. As a consequence, the loss of hyperbolicity may lead to the existence of nonclassical shocks (also called transitional shocks or undercompressive shocks) in three-phase flow solutions. The numerical procedure was used in an investigation of the existence and stability of nonclassical shocks with respect to viscous fingering in heterogeneous two-dimensional flows, thereby extending previous results for one-dimensional three-phase flow available in the literature. Numerical experiments, including the study of Water-Alternating-Gas (WAG) injection strategies, indicate that the proposed new numerical procedure leads to computational efficiency and accurate numerical results. Directions for further research are also discussed, based on the developments reported in this thesis.