Cette thèse a pour objectif le développement et l'implémentation d'un outil de résolution éléments-finis pour les problèmes d'extrusion de matériaux visco-plastiques fondus qui soit adapté à l'industrie du pneumatique. Dans ces travaux nous présenterons une méthodologie complète pour résoudre un modèle couplé comprenant le calcul des écoulements dans l'outillage, la répartition de température dans le domaine et la position d'interfaces entre plusieurs matériaux non miscibles. Pour cela, nous utiliserons les équations de Stokes non linéaires, l'équation de la chaleur, la méthode des interfaces diffuses et l'adaptation de maillage hiérarchique. Après avoir détaillé les différentes méthodes utilisées pour la résolution des différents problèmes, nous présenterons des résultats numériques dans plusieurs géométries d'outillages liés à l'extrusion des polymères. / The objective of this thesis is the development and the implementation of a finite element resolution tool for extrusion problems of molten visco-plastic materials which is suitable for the tire industry. In this works we will present a complete methodology to solve a coupled model including the calculation of the flow field in the equipment, the temperature distribution in the domain and the position of interfaces between several immiscible materials. To achieve this goal, we will solve the nonlinear Stokes equations, the heat equation, the diffuse interface method and we will use hierarchical mesh adaptation. After describing the various methods used for the resolution of the different problems, we present numerical results in several industrial geometries of tools related to the extrusion of polymers.

Ce mémoire de maitrise présente une approche innovante de résolution de systèmes multi échelles intitulée « Finite-Element-Self-consistent Clustering Analysis (FE_SCA) » qui peut être traduite par Élément Fini couplé à un modèle réduit par regroupement auto cohérent des volumes élémentaires représentatifs. Cette méthode appartient à la catégorie des méthodes de mécanique computationnelle en homogénéisation multi-échelle. Étant un bon substitut à la méthode EF², à chaque point d'intégration du modèle macroscopique, un volume élémentaire représentatif réduit de la microstructure est résolu à l'aide des équations de Lippmann Schwinger. Ces dernières, étant non-linéaires, un schéma de résolution implicite est considéré. Dans ce travail, la méthode itérative de Newton-Raphson a été utilisée. La réduction des volumes élémentaires représentatifs est faite avec la méthode de regroupement en K-moyennes (k-means clustering) formant des groupements ayant le même comportement mécanique. Cette première phase s'appelle la phase hors-ligne, et vise à préparer une base de données qui comporte les informations compressées de la microstructure. La phase hors-ligne se compose de trois étapes : 1-Construction de la géométrie et homogénéisation de la microstructure avec la méthode Fast Fourier Transform (FFT) pour obtenir les tenseurs de concentrations de contraintes Aᵢⱼ, 2- regroupement des zones ayant les tenseurs Aᵢⱼ identiques, 3- Créations des tenseurs d'interactions Dᴵᴶ entre les groupements formés. Ensuite, cette base de données composée des tenseurs d'interaction Dᴵᴶ et de la répartition volumique des groupements formés est utilisée dans la deuxième partie de résolution multi échelle nommée la phase en-ligne. L'implémentation de la méthode est faite à partir des deux logiciels commerciaux suivants : Python pour la partie hors-ligne et ABAQUS/explicit pour la création du VUMAT de la partie en-ligne. La méthode est testée sur un matériau composite unidirectionnel avec deux phases ayant un comportement élastique. / This master's thesis presents an innovative approach of multi-scale simulation entitled Finite-Element-Self-consistent Clustering Analysis (FE_SCA). This method fits into the category of numerical methods for multi-scale homogenization. Being a good substitute for FE², at each integration point of the macroscopic model, a Reduced Order Model (ROM) of the microstructure is solved using the Lippmann Schwinger equations. As the latter are non-linear, an implicit solution scheme is considered. In this work, the iterative Newton-Raphson method was used. The reduced order model of Representative Volume Element (RVE) is carried out using the k-means clustering method, forming groups with the same mechanical behavior. This first phase is called the Offline phase and aims to prepare a database containing compressed microstructure information. The Offline phase consists of three steps: 1-Construction of the geometry and homogenization of the microstructure using the Fast Fourier Transform (FFT) method to obtain stress concentration tensors Aᵢⱼ, 2- Clustering of voxels according to tensors Aᵢⱼ, 3- Creation of interaction tensors Dᴵᴶ between each pair of clusters I,J. This database of tensors Dᴵᴶ and the volume distribution of the clusters is then used in the second part of the multi-scale resolution process, called the Online phase. The method is implemented using the following two commercial software packages: Python for the OFFLINE phase and ABAQUS/explicit for the Online phase using a user material subroutine VUMAT. The method is tested on a unidirectional composite material with two elastic phases.

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é.

No description available.

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.

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.

Ce projet de maîtrise présente l'intégration d'une méthode d'optimisation topologique dans le processus de Conception et Fabrication Assistée par Ordinateur. Il fut réalisé dans le cadre d'un projet multidisciplinaire issu d'une collaboration entre le groupe de recherche en optimisation des structures de l'Université de Sherbrooke (OptiS) et l'Équipe de Recherche en Intégration CAO-Calcul de l'UQTR (ÉRICCA).Ce projet multidisciplinaire consiste à développer un gratuiciel multiplateforme d'optimisation des structures intégrant la CAO à l'optimisation afin de permettre le développement complet de pièces ou structures mécaniques en partant du modèle CAO initial, jusqu'au modèle CAO final optimisé. Deux objectifs principaux sont visés dans le cadre de ce projet de maîtrise. Premièrement, implanter la méthode d'optimisation topologique par homogénéisation (méthode SIMP) pour des structures quelconques en 3D. Deuxièmement, développer une méthode de lissage pour réduire le bruit présent sur le maillage optimisé résultant de l'optimisation topologique par la méthode SIMP. Nous avons atteint ces deux objectifs en développant un processus d'optimisation complètement automatique en sept étapes. Elles correspondent respectivement à la modélisation géométrique, l'entrée des données initiales du problème (conditions aux limites, matériau, etc.), la sous-division de la géométrie en sous-domaines de design et de non-design, le maillage automatique adapté aux sous-domaines multiples, l'optimisation topologique, le lissage du maillage de surface et finalement la reconstruction de la géométrie finale. Les résultats ont démontré que notre implantation de la méthode SIMP fonctionne et donne des résultats très intéressants qui s'apparentent aux résultats présentés dans la littérature. Néanmoins, le développement d'une méthode de lissage de triangulation basée sur les méthodes classiques a démontré que ces méthodes sont très mal adaptées à des maillages très bruités tels que ceux obtenus avec la méthode SIMP. En somme, ce projet a permis de faire un grand pas vers l'intégration complète de l'optimisation comme une étape à part entière du processus de CAO/FAO.