Les garnitures mécaniques sont utilisées dans de multiples applications pour réaliser l'étanchéité autour d'arbres en rotation. Ces composants peuvent fonctionner efficacement pendant plusieurs années en conditions stables, mais leur durée de vie est significativement réduite lorsque les conditions varient. L'objectif de ce travail de recherche est de développer et d’utiliser un banc d'essais et code de calcul pour étudier l'impact de pulsations de pression, d’inversions de pression et du chargement dynamique résultant sur les performances de garnitures mécaniques ayant des faces mésalignées et présentant des défauts de planéité.Le solveur fluide d'un modèle numérique de garnitures mécaniques a été étendu aux conditions transitoires. Un module résolvant la dynamique des forces et des moments a été ajouté afin de prédire le déplacement axial et les déplacements angulaires de la face montée de manière flexible. Afin de caractériser les performances de garnitures, un banc d'essais générant des pulses de pression a été instrumenté et des méthodes de mesure de perte de volume d'huile et d'entrée d'eau mises en place.Des garnitures mécaniques à faces parallèles puis mésalignées, fonctionnant sous pulsations et inversions de pression, ont été testées expérimentalement et simulées. Seules de très faibles augmentations d'eau dans l'huile ont été observées, sans augmentation au cours du temps, et sans fuite d'huile mesurable. Les faibles valeurs d'entrées d'eau sont dues à la faible épaisseur de film et à la courte durée des inversions de pression. Une garniture mécanique expérimentale à fort défaut de planéité a aussi été testée. Contrairement aux autres paramètres, le défaut de planéité semble augmenter significativement la fuite et promouvoir les entrées d'eau et pourrait ainsi être à l'origine de certaines défaillances. / Face seals are mechanical devices used to seal rotating shafts in numerous applications. While they can operate efficiently under steady conditions for years, they tend to fail prematurely when operating in severe, or rapidly varying conditions. The focus of this research work is the development and use of an experimental and a numerical method to investigate the impact of pressure pulses, pressure inversions and induced dynamic loading on the performance of mechanical face seals exhibiting face misalignment and waviness.The fluid solver of a state-of-the art face seal numerical model was extended to transient conditions and a module solving the dynamics for the axial and angular degrees of freedom of the flexibly-mounted stator added. A system-level experimental setup generating pressure pulses was instrumented and methods to characterise face seal performance in terms of oil volume loss and ingression of water outer-fluid selected and implemented.Face seals, with flat and misaligned faces, operating under pressure pulses and pressure inversions were experimentally tested and simulated. They show only slight increase of water in the oil, no increase over time, and no measurable oil leakage. The low water ingression is due to low film thickness combined with the short duration of pressure inversions. An exploratory face seal of high waviness was also experimentally tested. Contrary to the other parameters, the waviness appears to significantly increase the leakage and promote water ingression and could thus be at the origin of some seal failures.

La plupart des problèmes d'ingénierie des deux derniers siècles ont été résolus grâce à des modèles structuraux pour poutres, plaques et coques. Les théories classiques, tels que Euler-Bernoulli, Navier et de Saint-Venant pour les poutres, et Kirchhoff-Love et Mindlin-Reissner pour plaques et coques, ont permis de réduire le problème générique 3-D, dans le problème unidimensionnel pour les poutres et deux dimensionnelle pour les coques et les plaques. Théories raffinés d'ordre supérieur ont été proposées au cours du temps, comme les modèles classiques ne consentez pas à d'obtenir une complète domaine des contraintes et des déformations. La Carrera Unified Formulation (UF) a été proposé au cours de la dernière décennie, et permet de développer un grand nombre de théories structurelles avec un nombre variable d'inconnues principales au moyen d'une notation compacte et se référant à des nuclei fondamentales. Cette formulation unifiée permet de dériver carrément des modèles structurels d'ordre supérieur, pour les poutres, plaques et coques. Dans ce cadre, cette thèse vise à étendre la formulation pour l'analyse des structures fonctionnellement gradués (FGM), en introduisant aussi le problème thermo-mécanique, dans le cas des poutres fonctionnellement gradués. Suite à la formulation unifiée, les variables génériques déplacements sont écrits en termes de fonctions de base, qui multiplie les inconnues. Dans la deuxième partie de la thèse, de nouvelles fonctions de bases pour la modélisation des coques, qui représentent une approximation trigonométrique des variables déplacements, sont pris en compte / Most of the engineering problems of the last two centuries have been solved thanks to structural models for both beams, and for plates and shells. Classical theories, such as Euler-Bernoulli, Navier and De Saint-Venant for beams, and Kirchhoff-Love and Mindlin- Reissner for plates and shells, permitted to reduce the generic 3-D problem, in onedimensional one for beams and two-dimensional for shells and plates. Refined higher order theories have been proposed in the course of time, as the classical models do not consent to obtain a complete stress/strain field. Carrera Unified Formulation (UF) has been proposed during the last decade, and allows to develop a large number of structural theories with a variable number of main unknowns by means of a compact notation and referring to few fundamental nuclei. This Unified Formulation allows to derive straightforwardly higher-order structural models, for beams, plates and shells. In this framework, this thesis aims to extend the formulation for the analysis of Functionally Graded structures, introducing also the thermo-mechanical problem, in the case of functionally graded beams. Following the Unified Formulation, the generic displacements variables are written in terms of a base functions, which multiplies the unknowns. In the second part of the thesis, new bases functions for shells modelling, accounting for trigonometric approximation of the displacements variables, are considered.

Les alliages de titane β métastables sont des matériaux de structure essentiels pour les applications aéronautiques de part leurs très bonnes propriétés mécaniques. En effet, ils présentent une résistance spécifique élevée, une bonne ductilité et forgeabilité et une excellente réponse aux traitements thermiques. Toutefois, il existe encore aujourd'hui à leur sujet des controverses et des questions ouvertes et ce, malgré les efforts pour comprendre les mécanismes d'évolution microstructurale au cours de traitements thermo-mécaniques et pour déterminer les phases en présence et leur contribution aux les propriétés mécaniques résultantes. Ce travail de thèse a pour objectif de déterminer la nature de la phase β et de caractériser la transformation β→α à haute et basse températures par des caractérisations fines en microscopie électronique à balayage et à transmission couplées à des mesures d'orientations cristallographiques et de composition chimique. L'alliage étudié est un Ti-5553 avec une microstructure initiale 100% β obtenue par mise en solution et trempe. Il a été démontré expérimentalement que la structure de la phase β métastable n'est pas purement cubique centrée. Les points de la phase β dans les clichés de diffraction présentent un allongement (streaking) et des points supplémentaires sont visibles aux positions de diffraction 1/2, 1/3 et 2/3. Par ailleurs, les images MET ont un aspect en moiré. A partir de ces résultats et de calculs crystallographiques, il a été prouvé que des déplacements atomiques sur les plans {110}β et {112}β forment une structure intermédiaire entre celle de la phase β parente et celles des phases α et ω, prouvant que la phase β a intrinsèquement initié une transformation. L'étude de la précipitation au cours du procédé thermomécanique dans le domaine α+β a révélé que des précipités α discontinus, équiaxes ou légèrement allongés (1~2μm) se forment aux joints β de forte et de faible désorientation mais rarement au coeur des grains β produisant ainsi une microstructure en "collier". La relation d'orientation de Burgers (ROB) entre les phases α et β est progressivement détruite par la déformation. L'écart à la ROB est plus marqué pour les précipités α qui se forment au joint de grains qu'à l'intérieur des grains. L'écart à la ROB augmente aussi avec la déformation, mais diminue avec la vitesse de déformation. Au cours des déformations en bas du domaine α+β, les précipités α ont une morphologie qui dépend de leur position. Au coeur des bandes de glissement, les grains α/β sont équiaxes et ne respectent pas la ROB. Entre les bandes de glissement, la microstructure est lamellaire où les phases α/β alternent et respectent la ROB. Dans ce dernier cas, une forte sélection de variantes a été observée: Seuls les deux ou trois variants favorisant l'accommodation de la déformation se sont formés. A titre de comparaison, dans l'état non déformé, les 12 précipités sont présents. La transformation β→α est retardée en cours de compression à haute température. Ceci est attribué à une compétition entre adoucissement et transformation de phase. Au contraire, celle-ci est favorisée au cours de la compression à plus basse température du fait que les défauts cristallins induits par la déformation jouent le role de sites de germination et que la croissance des précipités soit accéléré alors que l'adoucissement soit ralenti. Dans le Ti-5553, le mécanisme de déformation dominant est le glissement des dislocations. Dans les déformations en bas du domaine α+β, du glissement simple ou multiple avec deux ou trois systèmes de glissement activés. L'identification de ces systèmes a pu être effectuée par des analyses de traces. Cette thèse a résolu la nature de la phase β métastable et constitue un travail de référence pour l'étude de la transformation β→α au cours de traitement thermomécanique / Metastable β titanium alloys are important structural materials for aeronautical applications due to their high strength to density ratio, good ductility and workability and excellent hardenability. Despite the efforts in resolving the complex microstructural evolution related to thermomechanical processes and in gaining knowledge on the produced phases and their contribution to the resultant mechanical properties, there are still some controversial and unresolved issues. The aim of the present PhD work is to determine precisely the metastable nature of β phase and to characterize finely the characteristics of the β→α transformation during high and low temperature thermomechanical treatments. Investigations were performed on a Ti-5553 alloy with the single β phase initial microstructure obtained by solution treatment followed by quenching using scanning and transmission electron microscopy (SEM/TEM) coupled to crystallographic orientation measurements and chemical analyses. It was demonstrated experimentally that the structure of the β phase in the metastable titanium alloy is not “pure” body centered cubic. Diffraction diagrams presents streaking of the β diffraction spots and additional spots at the 1/2, the 1/3 and 2/3 diffraction positions. Also, striations are observed in TEM images. From this experimental evidence and crystallographic calculations, it was proved that atomic displacements on the {110}β and {112}β planes formed a structure between that of the parent β phase and that of the α or ω phase, demonstrating pre-phase transformation tendency. The study of the precipitation during thermomechanical processing at higher temperature in the α+β region revealed that discontinuous equiaxed or short rod shaped α precipitates (1~2μm) mainly form on the high angle and low angle β grain boundaries but seldom in β grain interiors, forming the “necklace” microstructure. The Burgers orientation relationship (BOR) between the α and β phases is destroyed gradually by the deformation. The BOR deviation of grain boundary α is larger than that of intragranular α. The deviation from the BOR increases both with the increasing strain and decreasing strain rate. During the deformation at the lower temperature in the α+β region, the α precipitates exhibit different morphologies: such as lamellar α, equiaxed α and irregular α depending on their localization. Within the slip bands, equiaxed α/β grains which do not respect the BOR are present. However, between the bands, lamellar α and β phases maintaining the BOR are distributed alternately. In that last case a strong variant selection is observed as only the two or three variants that form are those which can accommodate the macroscopic deformation. Comparatively, in absence of compression all 12 variants are formed. The β→α phase transformation is retarded during the hot compression at higher temperature region, which is attributed to the competition between softening and phase transformation. On the contrary, it is promoted during compression at lower temperature region due to the more inducted deformation defects acting as α phase nucleation sites and due to accelerating growth of α precipitates and retarded softening. Dislocation slip is the leading deformation mechanism for the Ti-5553 alloy. Under the lower temperature deformation condition, single or multiple-slip bands with two or three different activated slip systems would form during the hot deformation process. Identification of these slip systems have been done by trace analysis. These results provide new insights into the structural nature of β metastable phase and valuable reference for β→α phase transformation during thermo-mechanical treatment in metastable β titanium alloys

Orogenic belts are fundamental features of plate tectonics. The crustal structure of orogens around the world shows a wide range of deformation styles from narrow, asymmetric wedges like the Pyrenees to wide, plateau-like orogens such as the Zagros or the Himalaya. The primary controlling factor on the size and structure of an orogen is the amount of convergence between the colliding plates. However, there are important additional factors providing major controls on the structural development of orogens. Among the potential parameters that can affect the style of deformation are the crustal strength, inherited weaknesses, and the surface processes. These parameters have been studied extensively in the past but their relative importance remains unclear. The aim of this thesis is to improve our understanding of: (1) How surface processes affect mountain building, with a special focus on the relationship between thin-skinned foreland and thick-skinned internal deformation of orogens. (2) How inherited extensional structures affect mountain building. The study was carried out using the Pyrenees as a special reference case. To answer our research questions we have used a wide range of state-of-the-art numerical modelling tools. In paper 1 we present a new method where we couple a structural-kinematic model and a thermo-kinematic model to evaluate the consistency of existing balanced section reconstructions with independent thermochronology data. In papers 2 and 3 we use 2D lithospheric scale thermo-mechanical models with surface process algorithms. Using the above toolset, we demonstrate that syntectonic sedimentation results in longer basement thrust sheets as well as longer thin-skinned thrust sheets and a generally wider orogen. Conversely erosion tends to narrow the wedge and reduce the orogenic loading of the colliding plates, limiting the space available for deposition in the flexural foreland deeps. We also demonstrate that inherited extensional structures play a crucial role in mountain building as they facilitate the migration of deformation into the undeformed basement of the overriding plate. Moreover, a significant amount of lower-crustal/mantle-lithospheric material is preserved at shallow depths only in the presence of extensional inheritance, but significant erosion is needed in order to bring this material to the surface. Our models also show that thin-skinned thrust sheets are generally rooted in the footwall of basement thrusts as they form outward-propagating sequences. As soon as a new basement thrust forms, the thin-skinned sequence situated on top of the new basement thrust-sheet is abandoned in favour of starting a new sequence in the footwall of the new thrust. Regarding our case study, it was possible to reproduce the section restoration using a structural-kinematic model with high accuracy up to the 36-Ma time slice and with limited accuracy up to the 50-Ma time slice. The thermochronometric ages predicted by the thermo-kinematic modelling are generally in good agreement with both the high- and low-temperature thermochronology data available in the Central Pyrenees; hence we conclude that the restoration is to a first order consistent with these datasets. The predicted thermochronological ages approximate the available low-temperature thermochronology data better by taking into account the late-stage burial and re-excavation scenario affecting the southern flank of the Pyrenean wedge presented by Coney et al. (1996), and quantified by Fillon and van der Beek (2012). In conclusion, our model experiments suggest, that extensional inheritance played a prime role in the structural evolution of the Pyrenees, with the major characteristics of the North Pyrenean Unit, including the presence of steep, inverted normal faults, the relative tectonic quiescence of the area after the early inversion and the presence of a lower-crustal body at shallow depth below the unit, best recaptured by our accordion models. / Orogenic belts are fundamental features of plate tectonics. The crustal structure of orogens around the world shows a wide range of deformation styles from narrow, asymmetric wedges like the Pyrenees to wide, plateau-like orogens such as the Zagros or the Himalaya. The primary controlling factor on the size and structure of an orogen is the amount of convergence between the colliding plates. However, there are important additional factors providing major controls on the structural development of orogens. Among the potential parameters that can affect the style of deformation are the crustal strength, inherited weaknesses, and the surface processes. These parameters have been studied extensively in the past but their relative importance remains unclear. The aim of this thesis is to improve our understanding of: (1) How surface processes affect mountain building, with a special focus on the relationship between thin-skinned foreland and thick-skinned internal deformation of orogens. (2) How inherited extensional structures affect mountain building. The study was carried out using the Pyrenees as a special reference case. To answer our research questions we have used a wide range of state-of-the-art numerical modelling tools. In paper 1 we present a new method where we couple a structural-kinematic model and a thermo-kinematic model to evaluate the consistency of existing balanced section reconstructions with independent thermochronology data. In papers 2 and 3 we use 2D lithospheric scale thermo-mechanical models with surface process algorithms. Using the above toolset, we demonstrate that syntectonic sedimentation results in longer basement thrust sheets as well as longer thin-skinned thrust sheets and a generally wider orogen. Conversely erosion tends to narrow the wedge and reduce the orogenic loading of the colliding plates, limiting the space available for deposition in the flexural foreland deeps. We also demonstrate that inherited extensional structures play a crucial role in mountain building as they facilitate the migration of deformation into the undeformed basement of the overriding plate. Moreover, a significant amount of lower-crustal/mantle-lithospheric material is preserved at shallow depths only in the presence of extensional inheritance, but significant erosion is needed in order to bring this material to the surface. Our models also show that thin-skinned thrust sheets are generally rooted in the footwall of basement thrusts as they form outward-propagating sequences. As soon as a new basement thrust forms, the thin-skinned sequence situated on top of the new basement thrust-sheet is abandoned in favour of starting a new sequence in the footwall of the new thrust. Regarding our case study, it was possible to reproduce the section restoration using a structural-kinematic model with high accuracy up to the 36-Ma time slice and with limited accuracy up to the 50-Ma time slice. The thermochronometric ages predicted by the thermo-kinematic modelling are generally in good agreement with both the high- and low-temperature thermochronology data available in the Central Pyrenees; hence we conclude that the restoration is to a first order consistent with these datasets. The predicted thermochronological ages approximate the available low-temperature thermochronology data better by taking into account the late-stage burial and re-excavation scenario affecting the southern flank of the Pyrenean wedge presented by Coney et al. (1996), and quantified by Fillon and van der Beek (2012). In conclusion, our model experiments suggest, that extensional inheritance played a prime role in the structural evolution of the Pyrenees, with the major characteristics of the North Pyrenean Unit, including the presence of steep, inverted normal faults, the relative tectonic quiescence of the area after the early inversion and the presence of a lower-crustal body at shallow depth below the unit, best recaptured by our accordion models.

Cette thèse CIFRE a été réalisée dans le cadre d'un partenariat entre l'ENTPE Université de Lyon, AREA et EIFFAGE Travaux Publics. Les méthodes actuelles de détermination et de prévision de l'endommagement des structures de chaussées autoroutières sont peu fiables et nécessitent d'être améliorées. L'objectif de la thèse est de mettre en place une procédure d'analyse inverse donnant accès aux paramètres de comportement des matériaux de chacune des couches de chaussées à partir des mesures de déflexion. Une première étape vise à définir plus précisément les états de contraintes et de déformations existants dans les différentes couches de chaussées. Pour cela, une instrumentation de chaussée autoroutière située proche d'Aix les Bains a été réalisée. 47 capteurs de déformation et de température ont été implantés dans la structure à trois niveaux différents lors de travaux d'entretien effectués en août 2012. En parallèle, des matériaux de chaussée ont été prélevés sur autoroutes pour caractériser leur comportement en laboratoire. Des essais de module complexe ont permis de déterminer les propriétés thermo-visco-élastiques des matériaux. Ces résultats ont permis de modéliser la déformation de la chaussée et de la comparer avec les mesures. Des mesures de déflexion ont également été réalisées sur la chaussée instrumentée à l'aide d'un curviamètre. Les résultats ont pu être comparés à la modélisation. Les différences observées ont permis de proposer une analyse critique des mesures de déflexion. / This thesis was conducted into a parnership between ENTPE university of Lyon, AREA and Eiffage Travaux Publics. Current methods of determination and prediction of damage to highway pavement structures are unreliable and need to be improved. The aim of the thesis is to establish an back analysis procedure from the deflection measurements fiving access to material behavior parameters of each pavement layers. A first step is to define more precisely the states of stress and deformation existing in the different layers of pavement. For this reason, a motorway pavement instrumentation was performed near Aix les Bains in France. 47 strain sensors and temperature were implanted into the structure in three layers during maintenance works in August 2012.In parallel, pavement materials were collected on motorways to characterize their behavior. Complex modulus tests were used to determine the thermo-viscoelalstic properties of materials. These results were used to modeling the deformation of the pavement and compared with measurements. Deflection measurements were also performed on the instrumented pavement with a curviameter device. The reseults were compared to modeling. The observed differences have helped to provide a critical analysis of deflection measurements.

Neste trabalho, propõe-se um modelo numérico-computacional representativo dos processos de tratamento térmico, que seja uma ferramenta eficiente e forneça meios para um entendimento efetivo do mecanismo de geração de tensões residuais durante a têmpera de aços. Foram investigados os fenômenos térmicos, mecânicos e de transformação de fase observados na têmpera, bem como o acoplamento entre esses três fenômenos. O modelo utiliza o Método dos Elementos Finitos (MEF) e o programa ABAQUS®, além de rotinas numéricas em FORTRAN responsáveis pela resolução do problema termo-mecânico-microestrutural acoplado. A utilização de sub-rotinas, que implementam uma alteração na formulação (matemática e numérica) do programa de Elementos Finitos, permite incluir no modelo as informações presentes em uma curva CRC (curva de resfriamento contínuo) do aço SAE 4140, implementando o cálculo de deformações da peça simulada de modo incremental e cumulativo. Os resultados mostram que a utilização das sub-rotinas desenvolvidas neste trabalho permitiu implementar, conjuntamente com o programa ABAQUS®, o cálculo das frações volumétricas, durezas, distorções e tensões que surgem em um tratamento térmico de têmpera, simulando as transformações martensítica, perlítica, bainítica e ferrítica. Os resultados dos modelos foram equivalentes aos relatados pela literatura, principalmente no que se refere às durezas e tensões associadas a cada transformação de fase. Em particular, os resultados indicam que a transformação martensítica está sempre associada à formação de tensões compressivas. Ensaios experimentais foram realizados a fim de validar os modelos computacionais propostos, utilizando-se um teste Jominy adaptado e instrumentado, de modo a permitir a amostragem da variação de temperaturas no material. Ensaios metalográficos permitiram correlacionar as frações volumétricas transformadas durante a têmpera do corpo de prova Jominy aos valores calculados pelo modelo numérico acoplado. / The objective of this work is to analyze residual strains and stresses and volumetric expansion due to phase transformations that occur during quenching of a steel body, as well as to predict these phase transformations. The coupled thermo-mechanical-phase transformation problem was analyzed, specifically in terms of the quenching process. Different computational models were presented, based on the finite element software ABAQUS® and on the use of FORTRAN subroutines. The continuous-cooling-transformation (CCT) diagrams of SAE 4140 steel are represented differently in each model, depending on the transformed phases and correspondent volumetric expansion. These subroutines include information from the CCT diagrams of SAE 4140 into a FORTRAN code. The subroutine calculates all the microstructures resulting from quenching (ferrite, pearlite, bainite, and martensite), depending on cooling rate. The numerical analysis conducted in this work provided results in terms of the temperature and stresses developed during quenching. The properties determined in this work are hardness, yield strength, volumetric fraction and distortion. Hardness has been predicted by the use of analytical equations. The finite element analyses were able to explain and reproduce phenomena observed during quenching of a steel cylinder. In particular, numerical results indicated that martensite formation is always related to a compressive stress field. The results of the models are in qualitative agreement with data provided by literature, particularly, in relation to the stresses originated by each different phase transformation during quenching process. Experimental testing was conducted, based on the analysis of the quenching of a Jominy probe, in order to validate the computational model developed in this work.

Le mémoire porte sur l’étude de la fiabilité des assemblages utilisant des alliages de brasure sans-plomb en environnement sévère pour des applications électroniques. Une méthode de préparation métallographique fiable et reproductible en vue de l’analyse microstructurale est présentée. L'effet du vieillissement thermique statique sur l’évolution microstructurale de billes de brasure en SAC, la tenue à la fatigue mécanique des assemblages brasés et l'étude de l'endommagement thermomécanique sont développés. Le seul effet notable constaté après vieillissement thermique est l'augmentation des épaisseurs des intermétalliques de contact. A la suite de sollicitations mécaniques, les ruptures sont plus souvent observées dans les pistes que dans les brasures, et ne permettent pas d'incriminer des éléments de la microstructure brasée comme facteurs déterminant de propagation des fissures. Les sollicitations thermomécaniques conduisent à une recristallisation, et les fissures se propagent alors le long des joints de grains recristallisés. / The thesis focuses on the study of the reliability of assemblies using Pb-free alloys under' harsh environment for electronic applications. A method for preparing metallographic reliable and reproducible for the microstructural analysis is presented. The effect of thermal aging on the static microstructural evolution of solder balls in SAC, the resistance to mechanical fatigue of solder joints and the study of thermomechanical damage are developed. The only significant effect observed after thermal aging is the increase in thickness of the intermetallic contact. Following mechanical, ruptures are most often seen in the copper tracks of the assembly, and do not blame the solder microstructure elements as determinants of crack propagation. The thermomechanical stresses lead to recrystallization, and then cracks propagate along the grain boundaries of recrystallized.

Cette thèse porte sur le thème général de la fiabilité des circuits microélectroniques. Le but de notre travail fut de développer un outil de simulation multi-physiques pour la conception des circuits intégrés fiables qui possède les caractéristiques innovatrices suivantes : • (i) L’intégration dans un environnement de conception microélectronique standard, tel que l’environnement Cadence® ; • (ii) La possibilité de simulation, sur de longues durées, du comportement des circuits CMOS analogiques en tenant compte du phénomène de vieillissement ; • (iii) La simulation de plusieurs physiques (électrique-thermique-mécanique) couplées dans ce même environnement de CAO en utilisant la méthode de simulation directe. Ce travail de thèse a été réalisé en passant par trois grandes étapes traduites par les trois parties de ce manuscrit. / This thesis was carried out under the theme of the microelectronics Integrated Circuits Reliability. The aim of our work was to develop a multi-physics simulation tool for the design of reliable integrated circuits. This tool has the following innovative features : • (i) The integration in a standard microelectronics design environment, such as the Cadence® environment ;• (ii) The possibility of efficient simulation, over long periods, of analog CMOS circuits taking into account the aging henomenon ; • (iii) The simulation of multiple physical behaviours of ICs (electrical-thermalmechanical) coupled in the same environment using the direct simulation method. This work was carried out through three main stages detailed in the three parts of this Manuscript.

Dans ce projet, le matériau est l'alliage de magnésium AZ31-O en tôle. L'épaisseur de tôles est de 1,2 mm. Les essais de traction à chaud sont réalisés afin d'étudier la ductilité de l'alliage de magnésium AZ31-O, la température et l'influence la vitesse de déformation sont incluses dans tous les tests. Le résultat d'analyse montre que la ductilité est renforcée avec une température croissante et une vitesse de déformation décroissante, le phénomène d'adoucissement est évident à la température élevée. La propriété anisotrope n'est pas considérée dans ce projet. Les essais Nakazima à chaud avec le poinçon d'hémisphère sont réalisés pour étudier la formabilité de l'alliage de magnésium AZ31-O. Enfin, la FLD (Forming Limit Diagram) est identifiée et les comparaisons montrent que la formabilité est préférable à une température plus élevée. En outre, les prédictions des limites de formage sont effectuées dans le modèle M-K. La comparaison montre clairement avec la prédiction théorique ne convient pas avec l'expérience. Les simulations des éléments finis sont effectuées pour un emboutissage par poinçon hémisphérique et un emboutissage en croix. Tout d'abord, les simulations d'emboutissage de poinçon hémisphérique sont réalisés sur FORGE® et sur ABAQUS®. Les résultats des simulations de FORGE et de ABAQUS sont comparés afin d'étudier la différence de divers codes de simulation des éléments finis. Deuxièmement, le comportement de d'endommagent est étudié dans FORGE par modèle d'endommagement Lemaitre. Enfin, la simulation d'emboutissage en croix qui est un benchmark de la conférence 2011 NUMISHEET est réalisée avec FORGE. La charge de poinçon, l'épaisseur et la distribution de température sont obtenues et comparées pour chaque simulation. En outre, ces résultats de la simulation de benchmark (FORGE) sont également comparés à d'autres logiciels de simulation en conférence. Les résultats des analyses détaillées sont présentés dans cette thèse. / In this project, the material is AZ31-O magnesium alloy sheet. The sheet thickness is 1.2mm. Warm tensile tests are performed to study ductility of AZ31-O magnesium alloy, the temperature and strain rate influence are included in all tests. The analysis result shows the ductility is enhanced with temperature increasing and strain rate decreasing, and the softening phenomenon is obvious at high temperature. The anisotropic property is not considered in this project. Warm Nakazima tests with hemisphere punch are performed to study formability of AZ31-O magnesium alloy. Finally, the FLD (Forming Limit Diagram) is identified and the comparisons distinctly show that the formability is better at higher temperature. Moreover, the forming limits predictions are performed in M-K model. The comparison clearly shows the theoretical prediction do not fit well with experiment. Finite element simulations are performed for a hemisphere punch deep drawing and a cross-shaped deep drawing. Firstly, the hemisphere punch deep drawing simulations are performed in FORGE® and ABAQUS®. The simulation result from FORGE and ABAQUS are compared in order to study the difference of various finite element simulation codes. Secondly, the damage behavior is studied in FORGE by Lemaitre damage model. Finally, the cross-shaped deep drawing simulation which is a benchmark of NUMISHEET 2011 conference is performed with FORGE. The punch load, thickness and temperature distribution are obtained and compared for each simulation. Furthermore, this benchmark simulation results (FORGE) are also compared with other various simulation software in conference. The detailed analysis results are presented in this thesis.

Electrical machines in electrified heavy-duty vehicles are subjected to dynamic temperature loadings during normal operation due to the different driving conditions. The Electrical Insulation System (EIS) in a stator winding is aged as an effect of these dynamic thermal loads. The thermal loads are usually high constant temperatures and thermal cycling. The high average constant thermal load is well-known in the electrical machine industry but little is known about the effect of temperature cycling. In this project, the ageing of the EIS in stator windings due to temperature cycling is examined. In this project, computational simulations of different simplified models that represent the electrical insulation system are made to analyse the thermo-mechanical stresses that is induced due to thermal cycling. Furthermore, a test object was designed and simulated to replicate the stress levels obtained from the simulations. The test object is to ease the physical testing of electrical insulation system. Testing a complete stator takes time and has the disadvantage of having a high mass, therefore a test object is designed and a test method is provided. The results from the finite element analysis indicate that the mechanical stresses induced will affect the lifetime of the electrical insulation system. A sensitivity study of several thermal cycling parameters was performed, the stator core length, the cycle rate and the temperature cycle amplitude. The results obtained indicate that the stator core length is too short to have a significant effect on the thermo-mechanical stresses induced. The results of the sensitivity study of the temperature cycle rate and the temperature cycle amplitude showed that these parameters increase the thermo-mechanical stresses induced. The results from the simulations of the test object is similar to the results from the simulations of the stator windings, which means that the tests object is valid for testing. The test method that is most appropriate is the power cycling test method, because it replicates the actual application of stator windings. The thermally induced stresses exposing the slot insulation exceeds the yield strength of the material, therefore plastic deformation may occur only after one thermal cycle. The other components in the stator are exposed to stresses below the yield strength. The thermally induced stresses exposing the slot insulation are high enough to low cycle fatigue the electrical insulation system, thus thermo-mechanical fatigue is an ageing factor of the electrical insulation system.