L’adaptation de maillages anisotrope permet d’améliorer grandement la précision de méthodes numériques telles que la méthode des éléments finis. L’adaptation peut être très coûteuse en temps de calcul, mais montre un grand potentiel de parallélisation. Ce projet consiste à développer une implémentation parallèle de l’adaptation de maillages anisotrope basée sur des estimateurs métrique et hiérarchique. On introduit les notions de base de l’adaptation et de la programmation parallèle. On développe une stratégie de parallélisation applicable à la fois en deux dimensions et en trois dimensions. Enfin, on applique cette stratégie dans le cadre de la méthode des éléments finis dans le but d’en évaluer l’efficacité et la rapidité.

Le comportement du mur de soutènement soumis à une excitation dynamique a été l'objet de nombreuses recherches. La première méthode est l'approche de Mononobe-Okabe (1929). Après cette méthode, plusieurs méthodes et approches ont été développées, notamment les approches pseudo-statiques et pseudo-dynamiques pour les murs libres, et l'approche de Wood (1973) pour les murs fixes. Ce mémoire présente une étude (2D) statique et dynamique d'un mur de soutènement qui soutient un remblai sableux. Le mur est considéré immobile à l'extrémité extérieure, tandis qu'à l'extrémité intérieure, il se déforme pour créer un domaine élastique dans le massif sableux. La réalisation de cette étude est assurée par des séries de simulations numériques avec les méthodes des éléments finis sous le logiciel commercialisé COMSOL. L'hypothèse sur le sable est un sol sec dont le contenu d'eau est nul. La validation du modèle est faite avec la comparaison des résultats obtenus avec ceux de Wood (1973), dans la littérature. Les analyses des domaines stationnaires, temporaires et fréquentielles ont été utilisées dans les différentes simulations dans COMSOL. Les principaux paramètres, qui ont été étudiés, sont : le coefficient de poisson, la longueur du remblai, la fréquence d'excitation, la rigidité verticale de l'élément d'interface, l'amortissement, les forces de gravité et l'amplitude d'excitation. Les résultats obtenus ont démontré que les pressions de tensions de Von Mises croissent avec l'augmentation du coefficient de poisson et de la longueur du remblai. La rigidité verticale d'élément d'interface engendre des tensions dans le domaine élastique principalement dans la partie supérieure du mur. L'absence de force de gravité constitue le cas le plus défavorable, surtout dans le volet dynamique pour le mur rugueux. C'est toutefois le contraire pour le mur lisse. L'étude fréquentielle a permis de déterminer que le système a un contenu fréquentiel bas. L'accélération spectrale a mis le point sur les caractéristiques de notre système, soit la période et la fréquence naturelles du système. / The behavior of the retaining wall under dynamic excitation has been the subject of many researches. The first method is the approach of Mononobe-Okabe (1929). After this method, several methods and approaches have been developed, including pseudo-static and pseudo-dynamic approaches for free walls, and Wood's (1973) approach for a fixed wall. This work presents a static and dynamic (2D) study of a retaining wall that supports a dry sand embankment. The wall is considered motionless at the outer end, while at the inner end, it deforms to create an elastic domain in the sandy mass. The realization of this study is ensured by a series of numerical simulations with finite element methods under the COMSOL software. The hypothesis on the sand is a dry soil whose water content is zero. The validation of the model is done with the comparison of the results obtained with those of Wood (1973) in the literature. Stationary, temporary and frequency domain analysis were used in the different simulations in COMSOL. The main parameters that have been studied are the poisson coefficient, the length of the backfill, the excitation frequency, the vertical stiffness of the interaction element, soil damping, the gravitational forces and the amplitude of the excitation. The results obtained showed that Von Mises stress pressures rise with the increase in the Poisson coefficient and the length of the embankment. The vertical rigidity generates tensions in the elastic domain mainly in the upper part of the wall. The absence of force of gravity constitutes the most unfavorable case especially in the dynamic shutter for the rough wall. It is, however, the opposite case for the smooth wall. The frequency study showed that the system has a low frequency content. The spectral acceleration has given the characteristics of our system; the period and the natural frequency of the system.

Tableau d’honneur de la Faculté des études supérieures et postdoctorales, 2014-2015 / Ce mémoire porte sur les transformations isoparamétriques quadratiques dans le contexte de la méthode des éléments finis, mais en gardant un point de vue théorique. Le but est de trouver des critères d’inversibilité afin de pouvoir appliquer des techniques d’adaptation de maillage avec des éléments quadratiques (courbes). Le premier chapitre introduit le sujet en présentant les difficultés qui surviennent avec les éléments quadratiques. Le deuxième chapitre présente une revue des principaux résultats sur le sujet et une caractérisation des conditions d’inversibilité en deux dimensions. Finalement, le troisième chapitre étudie ce qui se passe en trois dimensions et tente d’apporter des idées pour l’inversibilité dans un cas particulier.

Une conséquence du lemme de Céa permet de déterminer l’ordre de convergence optimal d’une solution estimée par la méthode des éléments finis. L’hypothèse importante qui est faite pour obtenir cette convergence est que la solution doit être suffisamment régulière. Dans ce mémoire, nous montrerons l’impact de l’adaptation sur la convergence de problèmes de diffusion ayant une singularité dans leur gradient. Afin de retrouver l’ordre optimal de la méthode sur des problèmes où cette hypothèse est fausse, nous développerons une méthode d’adaptation de maillage faite a posteriori, c’est-à-dire que l’on doit déjà avoir une approximation initiale de la solution pour l’utiliser. L’adaptation est dite hiérarchique et utilise le fait qu’une meilleure approximation peut être trouvée à partir de l’approximation initiale. Cette méthode est en développement depuis quelques années et a été étudiée en détail dans la thèse de Bois [2].

Les équations de Saint-Venant sont un système aux dérivées partielles jouant un rôle central dans la modélisation des écoulements océaniques. La méthode des éléments finis est particulièrement adaptée pour résoudre les équations de Saint-Venant car elle offre une grande flexibilité sur les domaines irréguliers ainsi qu’une variété d’espaces pour l’approximation de la solution. Or, la qualité de la solution numérique dépend de l’interaction entre ces espaces. Pour certaines combinaisons ou paires d’élément finis la solution numérique peut présenter des oscillations articiellement introduites par la discrétisation. Cette thèse porte sur le comportement numérique des solutions aux équations de Saint-Venant obtenues par différentes paires d’éléments finis. Tout d’abord, une étude sur la dispersion des ondes d’inertie-gravité est présentée pour une sélection de neuf paires d’éléments finis. Un ensemble de trois propriétés est ensuite mis en évidence afin que la discrétisation respecte le comportement des équations analytiques. Une méthode basée sur le calcul des noyaux est utilisée pour caractériser les modes stationnaires correspondant aux écoulements géostrophiques. Finalement, les espaces vectoriels de Raviart-Thomas et Brezzi-Douglas-Marini sont analysés. / The shallow-water equations system plays a central role in numerical oceanic models. The finite element method is particularly well suited to solve the shallow-water equations as it works on irregular meshes with a variety of approximation spaces. However, the behavior of the numerical solution highly depends on the interaction between these approximation spaces. For specific finite element pairs the solution may exhibit spurious oscillations induced by the discretization scheme. In this thesis, we analyze these oscillations for a wide selection of finite element pairs. The numerical dispersion of inertia-gravity waves is quantified with dispersion analyses. A constructive linear algebra approach is developed to compute the kernels of the discretized operators. The results are used to characterize the smallest representable vortices on both structured and unstructured meshes. A special attention is given to the Raviart-Thomas and Brezzi-Douglas-Marini approximation spaces.

No description available.

Cette thèse s'intéresse à la génération, à l'utilisation et à l'évaluation des maillages hex-dominants, composés d'hexaèdres et de tétraèdres, dans la cadre de la simulation numérique par la méthode des éléments finis. Les éléments finis hexaédriques sont souvent préférés aux éléments tétraédriques car ils offrent un meilleur ratio entre précision et temps de calcul dans un certain nombre de situations. Cependant, si la génération automatique de maillages tétraédriques est aujourd'hui un domaine bien maîtrisé, ce n'est pas le cas de la génération de maillages hexaédriques alignés avec le bord, qui reste un problème largement ouvert. En l'absence de progrès significatifs, les approches actuelles se contentent de maillages hex-dominants afin de tirer parti des performances supérieures des hexaèdres et de la flexibilité géométrique des tétraèdres, qui rend possible le maillage automatique. Dans une première partie, nous développons des algorithmes robustes pour la génération de maillages hex-dominants à partir de champs de directions, notamment pour l'isolement et le remplissage des régions difficiles à mailler (singularités et autres dégénérescences). Dans la seconde partie, nous essayons de déterminer dans quelles situations et dans quelle mesure les maillages hexaédriques, et hex-dominants générés précédemment, sont plus intéressants que les maillages tétraédriques. Ceci implique spécifiquement d'étudier plusieurs manières d'effectuer des simulations par éléments finis avec les maillages hybrides, dont une approche où nous utilisons des contraintes de continuité pour maillages non-conformes. Pour mesurer l'influence du maillage sur l'approximation des solutions, nous proposons une nouvelle méthode d'échantillonnage pour calculer très efficacement des distances globales entre solutions éléments finis définies sur des domaines compliqués / This thesis focuses on generation, usage and evaluation of hex-dominant meshes, which are made of hexaehedra and tetrahedra, in the context of the finite element method. Hexahedron finite elements are often preferred to tetrahedron elements because they offer a better compromise between accuracy and computation time in certain situations. However, if tetrahedral meshing is a well mastered subject, it is not the case of hexahedral meshing. Generating hexahedral meshes with elements aligned to the borders is still an open and difficult problem. Meanwhile, current automated approaches can use hex-dominant meshes in order to take advantage of both hexahedron accuracy and geometrical flexibility of tetrahedra. In the first part, we develop robust algorithms for the generation of hex-dominant meshes with elements aligned with the borders. Specifically, we propose a method to extract and fill the areas where hexahedral meshing is difficult (singularities and degeneracies). In the second part, we try to identify and to quantify the advantages of hexahedral and hex-dominant meshes over tetrehedral ones. This requires to study various ways to apply the finite element method on hybrid meshes, including one in which we propose to use continuity constraints on hexahedral-tetrahedral non-conforming meshes. To measure the impact of meshes on the finite element accuracy, we develop a new sampling method which allows to compute efficiently global distances between finite element solutions defined on complicated 3D domains

Cette thèse vise à développer et à évaluer une méthode de décomposition de domaine avec recouvrement dans le cadre de la modélisation du contrôle non destructif (CND) par courants de Foucault (CF). L'objectif d'une telle approche consiste à éviter le remaillage systématique de l'intégralité du domaine d'étude lors du déplacement de l'un de ses éléments constitutifs(par exemple, déplacement de la sonde CF au dessus de la pièce contrôlée). Plus précisément, il s'agit de concevoir une méthode de décomposition de domaine avec recouvrement qui s'appuie sur la théorie apportée par la méthode des éléments finis avec joints. En plus de s'affranchir de la contrainte d'une interface d'échange invariante avec le mouvement, la technique décrite dans ce travail réalise des transferts d'information réciproques entre les domaines. Cette étude présente les résultats théoriques ainsi que numériques liés à la simulation magnétodynamique. Par ailleurs, l'intérêt d'une telle méthode est illustré par des applications sur des configurations bidimensionnelles de CND par CF. / This thesis aims at studying and developing a domain decomposition method with overlapping subdomains for the modeling in eddy current (EC) non-Destructive testing (NDT). The idea behind such an approach is the possibility to avoid the systematic remeshes of the whole studying domain when some of its components are modified (for example the displacement of the coil above the conductor). More precisely, this work aims at designing a domain decomposition method with overlapping based on the theory of the mortar finite element method. In addition to remove the constraint owing to an coupling interface which is invariant with the displacement, the technique described, in this work, realizes reciprocal transfers of information between subdomains. This study presents the theoretical and numerical results attached to the magnetodynamic simulation. Moreover, the interest of such a method is illustrated by applications in some 2D modeling cases of EC NDT.

Nous construisons une méthode numérique fiable pour simuler un écoulement dans un milieu poreux modélisé par une équation elliptique. La simulation est rendue difficile par les hétérogénéités du milieu, la taille et la géométrie complexe du domaine de calcul. Un maillage d'hexaèdres réguliers ne permet pas de représenter fidèlement les couches géologiques du domaine. Par conséquent, nous sommes amenés à travailler avec un maillage de cubes déformés. Il existe différentes méthodes de volumes finis ou d'éléments finis qui résolvent ce problème avec plus ou moins de succès. Pour la méthode que nous proposons, nous nous imposons d'avoir seulement un degré de liberté par maille pour la pression et un degré de liberté par face pour la vitesse de Darcy, pour rester au plus près des habitudes des codes industriels. Comme les méthodes d'éléments finis mixtes standards ne convergent pas, notre méthode est basée sur un élément fini mixte composite. En deux dimensions, une maille polygonale est découpée en triangles en ajoutant un point au barycentre des sommets, et une expression explicite des fonctions de base a pu être obtenue. En dimension 3, la méthode s'étend naturellement au cas d'une maille pyramidale. Dans le cas d'un hexaèdre ou d'un cube déformé quelconque, la maille est divisée en 24 tétraèdres en ajoutant un point au barycentre des sommets et en divisant les faces en 4 triangles. Les fonctions de base de l'élément sont alors construites en résolvant un problème discret. Les méthodes proposées ont été analysées théoriquement et complétées par des estimateurs a posteriori. Elles ont été expérimentées sur des exemples académiques et réalistes en utilisant le calcul parallèle. / We develop a reliable numerical method to approximate a flow in a porous media, modeled by an elliptic equation. The simulation is made difficult because of the strong heterogeneities of the medium, the size together with complex geometry of the domain. A regular hexahedral mesh does not allow to describe accurately the geological layers of the domain. Consequently, this leads us to work with a mesh made of deformed cubes. There exists several methods of type finite volumes or finite elements which solve this issue. For our method, we wish to have only one degree of freedom per element for the pressure and one degree of freedom per face for the Darcy velocity, to stay as close to the habits of industrial software. Since standard mixed finite element methods does not converge, our method is based on composite mixed finite element. In two dimensions, a polygonal mesh is split into triangles by adding a node to the vertices's barycenter, and explicit formulation of the basis functions was obtained. In dimension 3, the method extend naturally to the case of pyramidal mesh. In the case of a hexahedron or a deformed cube, the element is divided into 24 tetrahedra by adding a node to the vertices's barycenter and splitting the faces into 4 triangles. The basis functions are then built by solving a discrete problem. The proposed methods have been theoretically analyzed and completed by a posteriori estimators. They have been tested on academical and realistic examples by using parallel computation.

Cette étude porte sur l'intégration de l'aluminium et de matériaux composites de type Polymère Renforcé de Fibre de Verre (PRFV) dans les ponts transportables ferroviaires. Ces matériaux ont un intérêt particulier en raison de leur résistance aux conditions environnementales rudes et leur légèreté. Cependant, la combinaison de ces matériaux légers avec les matériaux traditionnels dont l'acier est peu fréquente dans les structures de génie civil. De plus les références de conception fournissent très peu d'informations pour le dimensionnement de ce type d'assemblage. Il existe néanmoins des références de conception des assemblages simplement en aluminium ou simplement en PRFV. Cependant, en comparaison avec l'acier, les recommandations proposées pour l'aluminium et les PRFV sont basées sur un nombre restreint d'essais et le dimensionnement optimal d'un assemblage boulonné entre des pièces d'acier et d'aluminium ou de PRFV peut s'avérer difficile. Aussi, l'objectif principal de la thèse est d'identifier les paramètres géométriques optimum des assemblages boulonnés aluminium-acier et PRFV-acier et d'effectuer une analyse critique des recommandations de conception existantes et si nécessaire les améliorer et enfin, d'évaluer la contribution de la colle dans les assemblages boulonnés/collés. Sur la base de l'état des connaissances disponibles, des études expérimentales et/ou numériques des assemblages boulonnés, collés, et boulonnés/collés des plaques aluminium-acier et PRFV-acier ont été menées. Pour les assemblages boulonnés, des configurations à un boulon et deux boulons ont été expérimentalement et analytiquement analysées. Les résultats obtenus ont permis d'évaluer l'effet des paramètres géométriques tels que l'épaisseur des plaques, les pinces longitudinale et transversale, et les pas longitudinal et transversal sur la résistance de l'assemblage. Par la suite, les paramètres géométriques optimums des assemblages boulonnés en PRFV-acier et aluminium-acier ont été identifiés et l'analyse critique des équations de dimensionnement des assemblages en aluminium et en PRFV recommandés dans les références de conception a été effectuée. Dans le cas des assemblages aluminium-acier, une amélioration des équations de résistance à la pression diamétrale et de résistance en traction et cisaillement combinés ont été proposées. Pour les assemblages en PRFV-acier, une étude semi-empirique a permis de proposer une équation de résistance en section nette adaptée aux PRFV pultrudés. Pour les assemblages collés et boulonnés/collés, seules les configurations à un boulon ont été testées. Des adhésifs méthacrylate et époxyde ont été utilisés et l'effet de ces derniers sur les assemblages multi-matériaux a été évalué. Il est apparu que le type d'adhésif a un effet important tant sur la résistance de l'assemblage que sur son mode de rupture. En outre la rigidité des plaques a un grand impact dans la capacité portante de l'assemblage collé. La contribution de la colle dans un assemblage boulonné a été trouvée avantageuse seulement dans le cas de l'assemblage PRFV-acier. Dans le cas de l'assemblages aluminium-acier, pas d'accroissement significatif n'a été observé.