Il est bien connu que les alliages à mémoire de forme (AMF) sont considérés comme une classe particulière de matériaux qui peuvent retrouver une forme préalablement définie par simple chauffage. Cette propriété remarquable, appelée effet mémoire de forme, peut être exploitée dans la conception d'applications originales afin de trouver les solutions intéressantes aux problèmes rencontrés dans divers champs industriels. Dans notre travail, nous proposons une loi constitutive tridimensionnelle thermomécanique adaptée aux alliages à mémoire de forme à base de fer. Elle tient compte de l'effet de la transformation martensitique et des mécanismes de glissement plastique et leurs interactions. La formulation adoptée est basée sur une description micromécanique simplifiée. La comparaison entre les forces motrices et les forces critiques d'activation des mécanismes mis en jeu nous ont permis de déterminer le type de comportement induit pour un niveau de chargement donné. Nous avons adopté le schéma implicite d'intégration de Newton-Raphson pour la résolution de ce système. Les résultats obtenus pour des chargements thermomécaniques sont comparés à ceux obtenus expérimentalement. / It is well known that Shape Memory Alloys (SMA) are a particular class of materials that can recover a memorized shape by simple heating. This remarkable property, called the Shape Memory Effect (SME), can be exploited in the design of original applications in order to find attractive solutions to problems encountered in various industrial fields. We propose a thermo-mechanical three-dimensional constitutive law adapted to Fe-based shape memory alloys. It takes into account the effect of the martensitic transformation and the plastic slip mechanisms and their interaction. The adopted formulation is based on a simplified micromechanical description. The macroscopic behaviour is derived by considering the equivalent homogeneous effect on a representative volume element. The Gibbs free energy expression is defined. Thermodynamic driving forces are then derived and compared to critical forces leading to the constitutive equations solved by Newton–Raphson numerical scheme. Obtained results for thermo-mechanical loadings are compared to experimental ones.

Les phénomènes de localisation de déformation et de coalescence de microfissures sont considérés comme les deux principaux précurseurs du processus de rupture des structures en matériaux cimentaires. Il est aujourd’hui largement reconnu que la rupture survient suivant les trois étapes successives que sont la nucléation des microfissures, la propagation des microfissures et leur coalescence à la présence d’une zone fortement localisée. La modélisation de ce processus de rupture présente un travail complexe et délicat. Cette thèse porte sur la modélisation multi-échelle de l’endommagement et de la rupture des matériaux cimentaires sous sollicitation hygrométrique avec des effets d’agrégat. Les travaux présentés dans cette thèse visent à décrire le comportement mécanique et le processus de rupture des matériaux cimentaires en prenant en compte des effets d’agrégats à une échelle mésoscopique. La pâte cimentaire sera considérée comme un milieu homogène qui sera décrit par un modèle de comportement élastoplastique homogénéisé, permettant de prendre en compte la localisation de déformation et le processus de fissuration. Une approche numérique, à l’aide de la méthode XFEM, sera proposée afin de prendre en compte les effets des agrégats et des interfaces. / The phenomena of strain localization and coalescence of microcracks are regarded as two principle precursors of the rupture process structures in cementitious materials. It is now widely recognized that the rupture occurs following by three stages that are microcracks nucleation, the propagation of microcracks and their coalescence in a highly localized area. The modeling of this process of rupture presents a complex and delicate work. This thesis focuses on the multi-scale modeling of damage and fracture of cementitious materials under hydrometric solicitation with effects of aggregates. The work presented in this thesis aims to describe the mechanical behavior and the rupture process of cementitious materials taking into account the effects of aggregates at mesoscopic scale. The cement paste is considered as a homogeneous medium which is described by an elastoplastic homogenized model, which takes into account the strain localization and cracking process. A numerical approach by using the extended finite element method(XFEM) will be proposed to take into account the effects of aggregates and interfaces.

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.

Le pliage est un procédé a priori élémentaire de mise en forme des tôles métalliques. La complexité des pièces techniques actuelles rend cependant pointues les opérations de pliage industriel. Le retour élastique apparaissant sur les pièces complexes après retrait des outils ne peut plus être anticipé simplement. Le problème posé est alors typique en mécanique : trouver la géométrie finale du flan en fonction des conditions aux limites imposées en forces et déplacements. Formellement, cela revient à résoudre un système d'équations différentielles dépendant du temps et des variables d'espace dont les conditions initiales sont données par les conditions aux frontières du volume. Le modèle de plasticité choisi pour le matériau et son identification sont déterminants pour la qualité de la solution. On réduit dans la mesure du possible la complexité du problème général par des hypothèses simplificatrices. L'intégration des équations réduites s'effectue alors analytiquement ou par une intégration numérique simplifiée. Toutefois les hypothèses nécessaires sont ici trop réductrices et limitent l'application à des cas d'école. Le pliage industriel ne peut finalement être modélisé qu'avec un code éléments finis. Celui-ci reproduit les phénomènes rencontrés lors d'un pliage industriel : courbure anticlastique, décalage de la fibre neutre, estampage... Une discussion est menée pour voir sous quelles conditions une simulation 2D diminuant les temps de calculs s'applique. Ces modèles fiabilisés sont utilisés couplés avec des méthodes d'analyses modales des défauts de forme pour optimiser les outillages.

Les précipités ductiles de niobium compris dans l'alliage à mémoire de forme Ni47Ti44Nb9 élargissent l'hystérésis de transformation après un traitement particulier. L'augmentation des températures de transformation inverse peut atteindre 80 °C. Cette caractéristique représente un grand intérêt pour des applications industrielles, notamment pour des bagues de serrage. Un modèle de comportement thermomécanique est présenté. Il s'appuie sur des observations expérimentales permettant de caractériser l'alliage, considéré comme un composite dont la matrice est composée de NiTi avec un comportement à mémoire de forme, et d'inclusions de niobium avec un comportement élastoplastique ductile. La technique de transition d'échelle de Mori-Tanaka entraîne la formulation d'une loi de comportement macrohomogène. Ce modèle est implémenté dans un code de calcul par éléments finis. Des dispositifs de serrage sont réalisés expérimentalement et modélisés numériquement afin de valider le modèle. / Niobium ductile precipitates included in Ni47Ti44Nb9 shape memory alloy broaden transformation hysteresis after a particular treatment. The increase of reverse transformation temperature can reach 80 °C. This feature is attractive for industrial applications, including tightening rings. A thermomechanical model is presented. It is based on experimental observations to characterize the alloy, considered as a composite whose matrix is composed of NiTi with shape memory alloy behavior and inclusions of niobium with elastic-plastic behavior. The Mori-Tanaka transition scale technique results in a macro-homogeneous behavior law formulation. This model is implemented in a finite element code. Clamping devices are realized experimentally and modelled numerically to validate the model.

Un modèle de comportement thermomécanique pour les alliages à mémoire de forme est présenté. Il prend en compte la transformation de phase martensitique, l'orientation des variantes de martensite ainsi que l'accommodation inélastique des macles au sein de la martensite formée sous une structure auto-accommodée. Un potentiel thermodynamique pour un volume élémentaire représentatif est proposé. Il est décrit à l'aide de trois variables internes définies à l'échelle macroscopique. Des forces thermodynamiques sont dérivées de ce potentiel et équilibrées en faisant intervenir des phénomènes dissipatifs. Le modèle est ensuite implanté dans un code de calcul par élément finis afin de dimensionner des structures en deux et trois dimensions. Ce modèle a servi par la suite de base à une description non locale du comportement superélastique permettant de prendre en compte les phénomènes de localisation observés dans les fils et les films minces d'AMF. Des éléments finis spécifiques sont développés afin de pouvoir prendre en compte ce type d'approche dans le cadre d'un calcul de structures. / A constitutive thermomechanical model for the behavior of shape memory alloys is presented. It takes into account the martensitic phase transformation, the orientation of martensite variants and the inelastic accommodation of twins inside self-accommodated martensite. A thermodynamical potential is built using three internal variables described at macroscopic scale. Driving forces are derived from this potential and the equilibrium is reached by considering dissipative phenomena. The model is then implemented into a finite element code in order to design two or three-dimensional structures. It was adopted as a fundamental for a non-local description of the superelastic behavior in order to take into account the localization phenomenon observed in SMA wires and thin films. Specific finite elements are developed to account with this type of approach in the framework of structures computation.

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