Le planage est un procédé de mise en forme utilisé dans l’industrie de l’aluminium dans le but de corriger les défauts de planéité et de réduire les contraintes résiduelles dans les bandes minces par des flexions sous traction. Après avoir étudié les possibilités pour modéliser analytiquement le passage d’une tôle dans une planeuse multi-rouleaux et établi les limites de cette approche, on propose dans ce travail de thèse un modèle éléments finis en trois dimensions pour simuler le fonctionnement d’une planeuse sous traction pure. Les distributions de déformations plastiques et de contraintes résiduelles, dans la largeur et dans l’épaisseur, sont calculées. Des défauts initiaux de planéité sont introduits dans la bande en entrée de la machine et on vérifie que ceux-ci sont bien corrigés à la sortie.Ensuite, pour limiter la taille du modèle et diminuer le temps nécessaire à la simulation, une seconde approche avec un chaînage de configurations élémentaires et un transfert de données est comparée au modèle complet d’étireuse. Une application de ce modèle alternatif est effectuée en étudiant l’influence du profil des rouleaux. Enfin, la déformée de la bande, après retour élastique et flambement éventuel, est prédite / Leveling is a forming process used in the aluminum industry in order to correct flatness defects and minimize residual stresses in thin metallic strips thanks to bending under tension. After testing an analytical model of multiroll leveling and raising the limits of this approach, this thesis manuscript introduces a three-dimensionnal Finite Element model to simulate the sheet conveying in an industrial configuration of stretcher. It can compute plastic strains and residual stresses through width and thickness. Initial flatness defetcs are taken into account in the entering strip and we verify how they are corrected at the exit. Then, to reduce model size and computation time, a second approach with two elementary configurations coupling and data transfer is compared to whole model of stretcher. This alternative model is applied to study the profiled rolls effect. Finally, it predicts the deformed strip after springback and potential buckling phenomena

Cette thèse concerne l’analyse d’écoulements laminaires et turbulents en régime permanent par la méthode des éléments finis. Le but est de présenter des nouvelles méthodologies qui permettent de contrôler efficacement et automatiquement certains paramètres du problème physique et numérique afin d’obtenir les résultats désirés par l’ingénieur. Ces résultats sont des besoins tangibles telle que la performance du système mesurée en termes de coefficients aérodynamiques, de pertes d’énergie, etc. Pour contrôler la précision de l’anticipation numérique, la simulation de l’écoulement n’est pas suffisante. Le lien entre les résultats et les paramètres est manquant. Une analyse de sensibilité permet d’obtenir cette information. Nous verrons que si cette analyse est bien faite, son coût devient marginal. Dans un premier temps, l’estimation numérique du résultat doit être fiable. Par conséquent, on présente une nouvelle façon de contrôler efficacement l’erreur liée à la discrétisation spatiale nécessaire à ce genre de méthodes computationnelles. La sensibilité du résultat par rapport aux résidus des équations de la physique est alors utilisée. Nous validons nos nouveaux outils sur des écoulements laminaires. Ensuite, nous continuons avec des écoulements turbulents. Une fois que ces quantités sont bien estimées numériquement, le problème peut être optimisé virtuellement. Ainsi, l’amélioration de la conception réelle peut être fortement accélérée. Dans ce travail, l’optimisation géométrique est abordée. Nous présentons de nouvelles façons de gérer les contraintes de conception. La sensibilité du résultat par rapport à des perturbations de la paroi est alors utilisée. / This thesis concerns the analysis of fluid mechanics by computer. The goal is to present new methodologies that allow controlling effectively and automatically parameters of the physical and numerical problem to obtain the results the engineer desires. These results are concrete needs as the performance of the system measured in terms of aerodynamic coefficients, energy losses, etc. To control the result, the flow simulation is not sufficient. The link between results and parameters is missing. A sensitivity analysis provides this information. If we do this analysis appropriately, its cost becomes marginal. First, the numerical estimate of the result must be reliable. Consequently, we give an new way to control effectively the error related to the spatial discretization necessary for this kind of computational methods. We use the sensitivity of the result with respect to the residuals of the physics equations. We are validating our new tools on laminar flows. Then we continue with turbulent flows. Once these quantities are well estimated numerically, the problem can be optimized virtually. Thus, the actual design improvement can be accelerated. In this work, geometric optimization is achieved. We present new ways to manage design constraints. We use the sensitivity of the result to wall perturbations.

Ce mémoire présente les travaux de mise en œuvre d'un modèle numérique par éléments finis de la cellule unitaire représentative d'un composite interlock 3D, à partir de tomographies. Ce matériau se compose de torons tissés dans les directions chaine et trame, ainsi qu'à travers l'épaisseur. Ce type d'interlock 3X a été largement étudié ces dix dernières années, sans pour autant que des modèles numériques performants soient réalisés, ce qui rend nécessaire une modélisation numérique de sa cellule unitaire. Il présente pourtant des avantages vis-à-vis des autres types de composites, notamment grâce à l'absence du phénomène de délamination, cause principale de dégradation des composites. Ce travail présente la démarche et les moyens mis en œuvre pour obtenir la cellule unitaire numérique du matériau. Des tomographies de différents états du tissu sec, non compacté et compacté, ont été réalisées afin de générer un modèle éléments finis à partir du logiciel de traitement d'image AVIZO ainsi que du logiciel de génération de textile TEXGEN. Un modèle simplifié d'une tranche de cellule unitaire a été choisi afin de simuler et valider le procédé de compaction numérique avec le logiciel élément finis Abaqus. Les résultats du modèle de compaction ont été comparés aux tomographies du tissu sec compacté. Les résultats de cette simulation simplifiée sont prometteurs. La déformation et les déplacements des torons obtenus sur la cellule unitaire du tissu sec simulé sont très proches de ceux observés sur les tomographies du tissu compacté de la sorte. / This thesis presents the work conducted toward the development of a finite element model of the representative unit cell of a 3D interlock composite, based on tomographies. This material is composed of yarns woven in the warp and weft directions, as well as through the thickness. This kind of 3X interlock has been very little studied to date, which makes numerical modelling of its unit cell necessary. However, it present advantages over other types of composites, in particular the absence of delamination phenomenon, the main cause of composite degradation. This work presents the approach and the means employed to obtain the numerical unit cell of the material. Tomographies of different states of the dry, uncompacted and compacted fabric were carried out in order to generate a finite element model from the image processing software AVIZO and the textile generation software TEXGEN. A simplified model of a unit cell slice was chosen in order to simulate and validate the numerical compaction process with the Abaqus finite element software. The results of the compaction model were compared with tomographies of dry texile. The results of this simplified simulation are promising. The deformation and displacements of the yarns obtained on the unit cell of the simulated dry textile are very close to those observed on tomographies of the compacted textile.

Cette thèse porte sur la modélisation numérique des problèmes de changement de phase. Le but est de proposer un algorithme robuste, basé essentiellement sur la méthode CutFEM, pour résoudre le problème avec la plus grande précision possible. Dans notre démarche, la méthode level-set est employée pour une représentation implicite de l'interface entre différentes phases, tandis que l'imposition des conditions essentielles sur cette dernière se fait à l'aide de la variante symétrique de la méthode de Nitsche. Le premier résultat de la thèse porte sur le développement d'un post-traitement fiable, facile à implémenter, pour une meilleure estimation du flux à travers une frontière immergée. La solution proposée est inspirée de la méthode de l'intégrale de domaine, à la différence que nous introduisons une stabilisation de type Ghost Penalty pour éliminer la sensibilité de la qualité de l'approximation à l'interface par rapport aux nœuds du maillage. L'idée est validée à travers une série d'exemples 2D et 3D.Le deuxième résultat concerne le développement d'un solveur CutFEM pour la simulation des problèmes de Stefan à deux phases sans convection. Dans notre approche, nous nous servons du premier résultat pour une meilleure évaluation de la vitesse de l'interface. Quant à la stabilisation des formulations d'extension de la vitesse et du transport de la fonction level-set, nous avons opté pour la méthode CIP (Continuous Interior Penalty) qui a déjà fait preuve d'une grande efficacité pour des problèmes fortement dominés par la convection. Des exemples numériques 2D permettent de confirmer l'efficacité de notre algorithme puisque des taux de convergence optimaux sont obtenus. Une simulation 3D est également proposée pour démontrer la validité du solveur dans une configuration plus réaliste. Pour le dernier résultat, nous reconsidérons le problème de Stefan, mais cette fois en prenant en compte la convection naturelle dans la phase liquide. L'équation de conduction de la chaleur est alors modifiée par l'ajout d'un terme de convection. De plus, cette équation est couplée aux équations de Navier-Stokes, introduisant de fait des non-linéarités supplémentaires. Un schéma itératif basé sur la méthode de Newton est alors proposé pour une résolution efficace du système. Les autres étapes de résolution sont assez similaires à celles de l'algorithme développé pour les problèmes de Stefan sans convection. Une étude numérique approfondie du solveur est présentée. Premièrement, l'efficacité de la discrétisation en temps et en espace est étudiée pour des problèmes (relativement simples) avec des solutions manufacturées. Ensuite, pour un problème beaucoup plus complexe modélisant la solidification de l'eau, une comparaison des résultats obtenus numériquement avec les résultats expérimentaux connus dans la littérature est présentée. Enfin, d'autres problèmes de changement de phase encore plus complexes tels que la fusion du gallium dans une cavité rectangulaire, ou encore celle du n-octadécane dans une cavité cubique sont considérés, et une comparaison des résultats de simulation avec les données expérimentales est également présentée. / This thesis focuses on the numerical modeling of phase change problems. The aim is to propose a robust CutFEM-based algorithm to solve the problem with the highest possible accuracy. In our approach, the level-set technique is used for an implicit representation of the interface, whereas the imposition of the essential conditions on the interface is done using the symmetric variant of the Nitsche method. The first result of the thesis is the development of a reliable, easy-to-implement post-processing for an accurate estimate of the flux across an interface. The proposed solution is inspired by the domain integral method, with the difference that we introduce Ghost Penalty stabilization to eliminate the sensitivity to the interface location relative to the mesh nodes. The accuracy of the method is validated through a series of two- and three-dimensional examples. The second result concerns the development of a CutFEM method for the numerical simulation of two-phase Stefan's problems without convection. In our approach, we use the idea of the first result for an accurate evaluation of the normal speed of the interface. Regarding the stabilization of the velocity extension and front transportation problems, we opted for the Continuous Interior Penalty method which has shown great effectiveness for convection-dominated problems. Two-dimensional numerical examples are used to confirm the effectiveness of our algorithm since optimal convergence rates are achieved. A three-dimensional simulation is also provided to demonstrate the validity of the solver in a more realistic configuration. For the last result, we reconsider the Stefan problem, but this time taking into account natural convection in the liquid phase. In this phase, the heat conduction equation is then modified by the addition of a convection term. Furthermore, this equation is coupled to the Navier-Stokes equations, thereby introducing additional nonlinearities. An iterative scheme based on Newton's method is then proposed for an efficient solution of the system. The other steps of the proposed solver are quite similar to those of the solver developed for Stefan problems without convection. An in-depth numerical study of the solver is presented. First, the efficiency of discretization in time and space is studied for (relatively simple) problems with manufactured solutions. Then, for a much more complex problem modeling the solidification of water, a comparison of the results obtained numerically with the experimental results known in the literature is presented. Finally, other even more complex phase change problems such as the fusion of gallium in a rectangular cavity, or that of n-octadecane in a cubical test cell are considered, and a comparison of the simulation results with experimental data is also provided.

Cette recherche vise à améliorer la compréhension et l'optimisation des mousses de caoutchouc, en particulier des mousses de caoutchouc naturel. Cela se fait par une intégration complète des résultats de laboratoire et des modèles par éléments finis (EF). L'objectif principal est de simuler, analyser et concevoir des mousses de caoutchouc avec des propriétés mécaniques sur mesure en utilisant des modèles EF pour diverses applications, telles que dans l'industrie du pneu. La nouveauté réside dans l'incorporation de la géométrie 3D réelle dans les modèles EF, abordant une lacune critique dans la littérature actuelle qui repose souvent sur des méthodes mathématiques ou analytiques simplifiées. L'approche interdisciplinaire du projet combine des données expérimentales et des méthodes numériques, offrant une stratégie efficace et économique pour le développement de produits en mousse de caoutchouc optimisés. Tout au long de la recherche, différents sous-objectifs sont poursuivis, notamment la modélisation des matériaux hyperélastiques non moussés, l'exploration de la modélisation géométrique 3D avancée, l'incorporation de la modélisation multi phase, et l'analyse de la concentration des contraintes et des points de rupture. Pour atteindre ces objectifs, différents facteurs affectant les propriétés mécaniques des mousses de caoutchouc ont été sélectionnés et étudiés. Ces facteurs incluent la densité relative de la mousse (0,3 à 0,5), la teneur en agent moussant dans les mousses de caoutchouc naturel (2 phr à 10 phr) et dans les mousses d'élastomère polyoléfinique (2 phr à 5 phr). La conclusion résume les principales résultats, souligne l'impact potentiel de la recherche et décrit les perspectives d'exploration dans le domaine des mousses de caoutchouc. / This work aims to advance the understanding and optimization of rubber foams, especially natural rubber (NR) foams. This is done through a comprehensive integration of laboratory results and finite element (FE) models. The man goal is to simulate, analyze, and design rubber foams with tailored mechanical properties using FE models for diverse applications, such as in the tire industry. The novelty lies in incorporating real 3D geometry into FE models, addressing a critical gap in the current literature that often relies on simplified mathematical or analytical methods. The project's interdisciplinary approach combines experimental data and numerical methods, providing an efficient and cost-effective strategy for the development of optimized rubber foam products. Throughout the research, different sub-objectives are pursued, including modeling hyperelastic unfoamed materials, exploring advanced 3D geometric modeling, incorporating multiphase modeling, and analyzing stress concentration and failure points. To achieve these objectives, different factors affecting the mechanical properties of rubber foams have been selected and investigated. These factors include the relative foam density (0.3 to 0.5), foaming agent content in natural rubber (NR) foams (2 phr to 10 phr) and in polyolefin elastomer (POE) foams (2 phr to 5 phr). The conclusion summarizes the key findings by underscoring the research's potential impact, and outlines avenues for future exploration in the field of rubber foams.

Une technique de renforcement des structures en béton consiste à lier rigidement à ces structures un matériau plus résistant à la traction, par exemple un tissu de fibres de carbone dans une matrice plastique. Le lien mécanique entre le matériau renforçant et la structure se fait actuellement par une colle de type réactive (époxy). Les différentes applications ont montré qu'un facteur limitant de ce procédé était la qualité de l'interface béton/colle à travers laquelle transitent les efforts mécaniques sous forme de contraintes de cisaillement. <br />On propose dans ce travail une modélisation de la phase de pénétration d'une colle fluide (colle à solvant) dans un matériau cimentaire (pâte de ciment). Une première partie consiste à définir les différents phénomènes physiques associés au processus de pénétration puis à établir le système d'équations différentielles décrivant analytiquement ces phénomènes. Une seconde partie décrit les moyens expérimentaux et les sources bibliographiques permettant de disposer des caractéristiques nécessaires à la résolution des équations. Le problème de la colle déposée en surface du matériau cimentaire est enfin étudié par une méthode numérique utilisant les concepts de volumes finis et d'éléments finis. Les résultats obtenus sont présentés et sont comparés avec différentes méthodes de validation expérimentales développées au cours de ce travail.

La notion de modules de Drinfeld est le centre de cette thèse, cette notion fut introduite par Drinfeld en 1973, comme étant des " modules elliptiques" appelés de nos jours modules de Drinfeld. Ceux sont des objets algèbriques analogues aux courbes elliptiques sur les corps des nombres et sur les corps finis, obtenus par la réduction modulo une place non-archimédiennene. Une étude de l'arithmétique de tels objet devient légitime, motivée par l'arithmétique des courbes définies sur un corps fini initiée par Artin, Hasse et Weil. Dans cette direction on pousse cette analogie, pour un module de Drinfeld de rang 2, à la majorité de points étudiés pour des courbes elliptiques sur un corps fini. On donne plus précisement un analogue du théorème de Weil, théorème de Deuring-Waterhouse, et un analogue du travail de S. Vladut sur la cyclicité de tel structure algébrique.

Dans cette these, nous etudions les schemas boite. Ils ont ete introduits par H.B. Keller en 1971. Dans un premier temps, on s'est interesse a des problemes elliptiques de type Poisson. Plusieurs schemas boite pour des domaines de $\mathbb(R)^2$ mailles par des triangles ou des rectangles ont ete introduits. Dans ce cas, la discretisation s'effectue sur la forme mixte du probleme en prenant la moyenne des deux equations (conservation et flux) sur les cellules du maillage. La methode peut etre qualifiee de ``methode volumes finis mixte de type Petrov-Galerkin ``. Une des difficultes du design de cette famille de schemas reside dans le choix des differents espaces de fonctions (approximation et test) qui doivent satisfaire des conditions de compatibilite de type Babuska-Brezzi. En revanche, cette methode de discretisation ne necessite qu'un seul maillage (le maillage du domaine). De plus, on montre dans la plupart des cas que le schema obtenu est equivalent a un probleme découplé : la résolution d'un probleme variationnel pour l'inconnue principale et une formule locale pour le gradient (le flux). Cette formulation facilite le calcul des inconnues discretes. Des resultats de stabilite et les calculs d'erreurs reposant sur la theorie des elements finis ont ete etablis. Une etude numérique valide ces resultats pour quelques cas tests. Dans le cadre du Groupement de Recherche MoMaS pour le stockage des dechets nucleaires dans la Meuse, j'ai ensuite etudie des problemes de convection-diffusion instationnaires. Un schéma boite permettant d'approcher ces equations dans le cas monodimensionnel a ete introduit. Des coefficients de decentrement propres a chaque maille permettent de controler le schema (precision, stabilite). Afin de generaliser rapidement ce schema au cas bidimensionnel, je me suis concentree sur une extension du schema boite monodimensionnel par la methode ADI (Alternating Direction Implicit).

Les méthodes numériques classiques sont inadaptées aux problèmes de diffusion au sein des matériaux cimentaires, du fait de l'écart entre l'échelle grossière de travail, de l'ordre du mètre, et l'échelle fine de la description du milieu, de l'ordre du micromètre, On présente dans cette thèse plusieurs méthodes de simulations multi-échelles couplant Volumes Finis et Éléments Finis, ainsi que leurs implémentations informatiques efficaces. En particulier, on a développé une chaîne de calcul multi-échelle utilisant la plate-forme SALOME (génération des maillages, post-traitement des solutions) et le code de calcul parallèle MPCube (résolution des problèmes) qui permet de réaliser automatiquement et efficacement des simulations multi-échelles. Un soin particulier a été apporté à la parallélisation des tâches et à l'optimisation des méthodes aux spécificités des matériaux cimentaires. On a appliqué cette chaîne de calcul à plusieurs échantillons de matériaux cimentaires, notamment des modèles de mortiers et de pâtes de ciment. Les résultats de ces simulations ont permis de déterminer une diffusivité numérique équivalente, ainsi que de reconstruire une solution à l'échelle fine.

Les modèles de turbulence de type Navier-Stokes en moyenne de Reynolds (RANS) au premier ordre sont étudiés dans cette thèse. Ils sont constitués des équations de Navier-Stokes, auxquelles on adjoint un système d'équations de bilan pour des échelles scalaires caractéristiques de la turbulence. L'évaluation de celles-ci permet, grâce à une relation algébrique, de calculer une viscosité additionnelle dite "turbulente", modélisant la contribution de l'agitation turbulente dans les équations de Navier-Stokes. Les problèmes d'analyse numérique abordés se placent dans le contexte d'un algorithme à pas fractionnaire constitué d'une approximation, sur un maillage régulier, des équations de Navier-Stokes par éléments finis non-conformes de Crouzeix-Raviart, ainsi que d'un ensemble d'équations de bilan de la turbulence de type convection-diffusion, discrétisées par la méthode de volumes finis standard. Un schéma numérique basé sur une discrétisation de volumes finis, permettant de préserver la positivité des échelles turbulentes telles que l'énergie cinétique turbulente (k) et son taux de dissipation (ε), est ainsi proposé dans le cas des modèles k − ε standard, k − ε RNG et leur extension k − ε − v2 − f. La convergence du schéma numérique proposé est ensuite étudiée sur un problème modèle constitué des équations de Stokes incompressibles et d'une équation de convection-diffusion stationnaires, couplées par les viscosités et le terme de production turbulente. Il permet d'aborder la difficulté principale de l'analyse d'un tel problème : l'expression du terme de production turbulente amène à considérer, pour les équations de bilan de la turbulence, un problème de convection-diffusion avec second membre appartenant à L1. Enfin, afin d'aborder le problème instationnaire, on montre la convergence du schéma de volumes finis pour une équation de convection-diffusion modèle avec second membre appartenant à L1. Les estimations a priori de la solution et de sa dérivée en temps sont obtenues dans des normes discrètes dont les espaces correspondants ne sont pas duaux. Un résultat de compacité plus général que le théorème de Kolmogorov usuel, qui se pose comme un équivalent discret du Lemme d'Aubin-Simon, est alors proposé et permet de conclure à la convergence dans L1 d'une suite de solutions discrètes.