Le développement de nouvelles technologies permettant d’obtenir des procédés de fabrication novateurs, de même que l’utilisation de nouveaux matériaux et lubrifiants, ont rendu possible une amélioration globale de l’efficacité des systèmes mécaniques et la résolution de nombreuses défaillances. Néanmoins, certains types de détérioration comme le grippage restent difficiles à appréhender car ils résultent de nombreux paramètres connexes tels que le régime de lubrification, la température du contact et le coefficient de frottement. Le grippage est un phénomène qui apparait de manière brutale et qui se caractérise par des arrachements et des transferts de matières entre les surfaces en contact pouvant aller jusqu’à la destruction complète du système. Ce dernier est généralement associé à une mauvaise dissipation des calories entrainant des températures de contact élevées. De nombreux critères basés sur la dissipation d’énergie dans le contact ont été établis sans pour autant que l’un d’eux ne soit unanimement reconnu par la communauté scientifique. La première partie de cette étude a consisté à caractériser le comportement en frottement d’un couple matériau-lubrifiant grâce à l’utilisation d’une machine bi-disque. En parallèle, une modélisation thermique de ce banc d’essais a été réalisée en utilisant la méthode des réseaux thermiques. Les corrélations entre les formules analytiques et les mesures expérimentales ont permis d’aboutir à l’établissement d’une loi de frottement liant conditions de contact et température de masse des éprouvettes. Il a aussi pu être démontré que la température de masse des disques pouvait être considérablement différente de la température d’injection du lubrifiant en fonction des conditions opératoires. Ce résultat permet une analyse nouvelle des critères de grippage dont la plupart assimilent la température des éprouvettes à celle d’injection du lubrifiant. Malgré cela, les différents essais de grippage réalisés montrent que l’unique considération de la température de contact pour établir un critère de grippage s’avère insuffisante. / The development of new technologies, which allow innovative manufacturing processes, as well as the use of new materials and lubricants have led to an overall improvement of mechanical systems efficiency and reliability. However some failures, like scuffing, remain difficult to understand since they depend on many related parameters such as the lubrication regime, the contact temperature and the friction coefficient. Scuffing is a critical damage that appears suddenly and which is characterised by local welds and scratches between the sliding surfaces. It can lead to complete destruction of the mechanical system. This phenomenon is in general related to poor heat dissipation and overheating resulting in high contact temperature. Numerous studies were conducted in order to establish a scuffing criterion based on energetic approach. But none of them is unanimously recognized by the scientific community. The first part of this study consisted in characterizing the frictional behaviour of a given couple of lubricant and material thanks to the use of a twin-disc machine. In a second part, a thermal modelling of this test bench has been realised using the thermal network methodology. Correlations between analytical formulas and experimental measurements allowed to establish a friction law relating contact conditions and disc bulk temperature. It has also been demonstrated that the disc bulk temperature could be considerably different from the lubricant injection temperature depending on the operating conditions. This result allows a new interpretation of existing scuffing criteria, which for most of them consider the disc temperature to equal the oil temperature. Despite this analysis, the scuffing experiments performed reveal that accounting solely for the contact temperature is not sufficient to establish a reliable scuffing criterion.

Chaque fois que se produit un mouvement relatif entre deux systèmes, avec une interface à contact sec, le contact frottant induit des vibrations. La dynamique locale au contact (ruptures et la génération d'ondes) se couple avec la dynamique du système, donnant origine à des vibrations et affectant le comportement frictionnel macroscopique du système. Dans cette thèse, afin de développer une approche globale pour l'investigation des phénomènes multi-physiques, l'énergie a été utilisée comme une caractéristique physique universelle du couplage. La formulation de un bilan énergétique mécanique est utilisé pour identifier deux termes dissipatifs différents, i.e. la dissipation par amortissement matériel/système et la dissipation au contact. Les flux d'énergie, provenant des surfaces en contact et dus aux vibrations induites par frottement, excitent la réponse dynamique du système et, vice versa, l'influence de la réponse dynamique du système sur la dissipation d'énergie locale à l'interface de contact affecte les phénomènes tribologiques connexes. Dans cette thèse, les vibrations induites par le frottement ont été analysées en utilisant: l'approche par éléments finis pour étudier, par l'analyse des flux d'énergie, le couplage entre le contact et la dynamique du système; l'approche expérimentale pour valider les résultats numériques et observer l'influence des phénomènes pas encore inclus dans les modèles numériques; une approche avec une modèle à paramètres concentrés pour évaluer rapidement les effets des paramètres du système. L'analyse numérique par le modèle éléments finis 2D permet une répartition de l’énergie introduite dans le système mécanique entre les deux termes dissipatifs (amortissement matériau et contact), au cours de la réponse transitoire aussi bien en conditions stables qu’instables. En particulier, les vibrations induites par frottement modifient la capacité globale du système à absorber et dissiper l’énergie; une estimation de la puissance dissipée au contact, sans prendre en compte le comportement dynamique du système (flux d’énergie par les vibrations induites par frottement) peut conduire à des erreurs significatives dans la quantification de l’énergie dissipée au contact, ce qui affecte directement plusieurs phénomènes tribologiques. Les mesures expérimentales de crissement montrent comment les mêmes modes instables sont reproduits soit expérimentalement soit numériquement, validant l’utilisation de la simulation 2D transitoires pour l’analyse des vibrations instables induites par le frottement. L’équilibre énergétique a été utilisé sur le modèle à paramétrés concentrés, pour approcher le problème de la surestimation d’instabilité, qui est caractéristique d’une analyse des valeurs propres complexes. Un nouvel indice d’instabilité (MAI) a été défini, par des considérations énergétiques, pour comparer les différents modes instables et pour sélectionner le mode qui devient effectivement instable pendant le crissement. / Whenever relative motion between two system components occurs, through a dry contact interface, vibrations are induced by the frictional contact. The local dynamics at the contact (ruptures and wave generation) couples with the system dynamics, giving origin to vibrations and affecting the macroscopic frictional behavior of the system. In this thesis, in order to develop an overall approach to the investigation of the multi-physic phenomenon, the energy has been pointed out as a coupling physical characteristic among the several phenomena at the different scales. The formulation of a mechanical energy balance is used for distinguishing between two different dissipative terms, i.e. the dissipation by material/system damping and the dissipation at the contact. The energy flows coming from the frictional surfaces, by friction induced vibrations, excites the dynamic response of the system, and vice versa the influence of the system dynamic response on the local energy dissipation at the contact interface affects the related tribological phenomena. The friction-induced vibrations have been analyzed using three different approaches: the finite element approach, to investigate the coupling between the contact and system dynamics by the analysis of the energy flows; the experimental approach to validate the numerical results and observe the influence of phenomena not still included into the numerical model; a lumped parameter model approach to quickly investigate the effects of the system parameters. The numerical analysis by the 2D finite element model allowed investigating the repartition of the energy introduced into the mechanical system between the two dissipative terms (material damping and contact) during both stable and unstable friction-induced vibrations. In particular, it has been shown how the friction-induced vibrations modify the overall capacity of the system to absorb and dissipate energy; an estimation of the power dissipated at the contact, without considering the dynamic behavior of the system (energy flows by friction induced vibrations) can lead to significant error in the quantification of the dissipated energy at the contact, which affects directly several tribological phenomena. The experimental squeal measurements show how the same unstable modes are recovered both experimentally and numerically, validating the use of the 2D transient simulations for the reproduction of the unstable friction-induced vibrations. Once the energy balance formulated, it has been used on the lumped model to approach the instability over-prediction issue characteristic of the complex eigenvalue analysis. By energy considerations, a newer instability index (MAI) has been defined to compare the different unstable modes and to select the mode that becomes effectively unstable during the transient response. The Modal Absorption Index allows quantifying the capability of each mode to exchange energy with the external environment.

Les pneus sont un organe déterminant dans la tenue de route des véhicules. Cette thèse porte sur la modélisation du contact entre la bande de roulement d'un pneumatique et une route sèche, afin de comprendre les différents phénomènes physiques mis en jeu ainsi que leurs rôles relatifs dans le frottement. La rugosité multi échelle des sols routiers les rendent difficiles à modéliser avec une simulation par élements finis standard. En utilisant l'hypothèse que la gomme de la bande de roulement est très grande devant la taille des rugosités, elle peut être considérée comme un massif semi-infini. Il est alors possible de résoudre efficacement le problème de contact en ne discrétisant que la surface du massif de gomme. Cette résolution est faite à l'aide d'un algorithme de Gradient Conjugué, au cours duquel les calculs matriciels sont effectués par Transformée de Fourier Rapide (FFT). La viscoélasticité de la gomme est prise en compte en régime transitoire. Les interactions à l'interface entre la gomme et le sol sont modélisés par une loi de frottement ainsi que par une énergie d'adhésion. Les résultats montrent le rôle primordial de la viscoélasticité qui, couplée à la rugosité multi-échelle du sol, modifie la surface du contact au cours des différents étapes de mise en glissement d'un pneumatique, faisant ainsi varier le frottement. / Tires are a key component for the handling and safety of personal vehicles. In this thesis a model of the contact between the tire tread and a dry road is described. It aims at understanding the different physical phenomena taking place in such a contact and their relative role in tire friction. Modeling the multiple scales of road roughness is difficult using a standard Finite Element Method. The hypothesis that the rubber of the tire tread is very large compared to the largest scale of surface roughness is made, so that it can be considered as a semi-infinite half-space. This way, the contact problem can be solved by discretizing the rubber surface only. The solver is a specific Conjugate Gradient iterative method, in which the matrix-vector products are performed with Fast Fourier Transforms. Transient viscoelasticity is accounted for with a step-by-step approach. The algorithm is able to model surface interactions such as Coulomb friction and adhesion. Results show the crucial role played by viscoelasticity. Coupled with the road roughness, it changes the contact surface during the different steps of tire sliding, which in turns impacts friction.

L’étude du frottement dans les contacts lubrifiés fortement chargés est un sujet complexe. En effet, le frottement dépend fortement du comportement rhéologique du lubrifiant dans le centre du contact et, celui-ci n'est pas bien connu à pression et taux de cisaillement élevé. Diverses méthodes expérimentales ont été développées pour expliquer le comportement plateau dans les courbes de frottement, associé à l’existence d’une contrainte limite de cisaillement, mais aucune d'elles n'a fourni une image complète des mécanismes réels impliqués. Dans la continuité de ces efforts, des études sont présentées dans ce mémoire. Le premier défi dans ce travail est d'effectuer des mesures de frottement dans des conditions isothermes nominales, ce qui signifie que même si des effets thermiques doivent se produire dans toute mesure de frottement, il est possible de les minimiser et de rendre les résultats insensibles à une faible dissipation d'énergie dans le volume expérimental d'intérêt. La minimisation de l’échauffement du lubrifiant aide à se focaliser sur l’origine mécanique de la contrainte limite de cisaillement et de mieux caractériser sa dépendance à la pression et à la température. C'est pourquoi, tout d'abord, une série d'expériences a été réalisée sur deux lubrifiants, un diester pur (benzyl benzoate) et une huile minérale de turbine (Shell T9) avec des vitesses d'entraînement variables. Cela nous permet d'abord d'observer directement l'influence des effets thermiques sur les valeurs de la contrainte limite de cisaillement et ensuite, de déterminer les conditions expérimentales qui limitent ces effets tout en assurant un régime en film complet. Le deuxième objectif est de caractériser le comportement en frottement des lubrifiants sous conditions isothermes nominales et sur une large gamme de pression (jusqu'à 3 GPa) et de température (jusqu’à 80°C) afin d'établir un nouveau modèle découplé permettant de décrire la dépendance à la température et à la pression de la contrainte limite sous fortes charges. Enfin, l'étude se focalise sur la compréhension du comportement microscopique des lubrifiants dans des conditions extrêmes de cisaillement et de pression. Des mesures in situ de spectroscopie Raman et Brillouin ont été également effectuées sous conditions statiques, afin d'étudier le changement de phase du lubrifiant. / Friction in highly loaded lubricated contacts is a complex issue. Indeed, it highly depends on the lubricant rheological behaviour in the Hertzian region, which is not well known under such high pressure and high shear stress. Various experimental methods have been developed to explain the plateau-like behaviour in friction curves referred to as the limiting shear stress (LSS), but none of them provided a full picture of the real mechanisms involved. In a continuation of these efforts, some investigations are presented in this manuscript. The first challenge in this work is to carry out friction measurements under nominal isothermal conditions, meaning that even if thermal effects must occur in any friction measurement, it is possible to minimize them and to make the results almost insensitive to a weak energy dissipation within the experimental volume of interest. Minimizing shear heating of the lubricant help us to focus on the mechanical origin of the LSS and to better characterize its dependence to pressure and temperature. That’s why, first of all, a series of experiments was performed on two lubricants, a pure diester fluid (benzyl benzoate), and a commercial turbine mineral oil (Shell T9) with varying entrainment velocities. This allow us first to directly observe the influence of the lubricant shear heating on the LSS values and then to determine the experimental conditions which limit this thermal effect while ensuring a full film regime. The second objective is to characterize the frictional behavior of both lubricants under nominal isothermal conditions and over a wide range of pressure (up to 3 GPa) and temperature (up to 80°C) in order to establish a new uncoupled model to describe the temperature and pressure dependence of the limiting shear stress under highly loaded conditions. Finally, the study focuses on the understanding of the microscopic behavior of lubricants under extreme shear and pressure conditions. In situ Raman and Brillouin spectroscopy investigations were also conducted under static conditions, in order to study the lubricant phase changes under various pressure and temperature conditions.

Cette thèse est dédiée à l’étude des contacts lubrifiés tore-plan sous diverses conditions. Ces contacts se situent à l’interface entre l’extrémité torique des rouleaux et le collet de la bague dans les roulements à rouleaux. La première complexité de cette étude provient de la géométrie particulière des solides concernés. La deuxième est générée par la cinématique complexe qui règne dans ces contacts. Afin de comprendre les mécanismes physiques à l’œuvre, une approche duale (expérimentale et numérique) est adoptée. Le banc d’essai Jérotrib permet une première étude basée sur l’hypothèse que le contact élastohydrodynamique tore-plan est similaire à un contact elliptique équivalent. Grâce à une méthode d’interférométrie optique en lumière blanche qui a été adaptée aux spécificités du contact en question, des mesures précises de l’épaisseur de film ont été effectuées dans un nombre significatif de conditions. Sur cette base, un modèle numérique thermo-élastohydrodynamique a été validé avec précision. Ce dernier a permis d’étudier les écoulements de fluide à l’entrée du contact afin de mettre en évidence leur influence sur le champ d’épaisseur de film. Le modèle numérique a ensuite été amélioré afin de prendre en compte la vraie forme des solides. Il a été validé en épaisseur de film par le banc d’essai Tribogyr, dans des conditions similaires à celles rencontrées dans les vrais roulements. Il a été montré que le cisaillement du fluide est responsable de l’échauffement des solides, qui diminue par suite l’épaisseur de film : ceci souligne la nécessité de modéliser cet échauffement global pour prédire la séparation des surfaces. Par ailleurs, lors de l’étude, le champ de pression et d’épaisseur de film ont perdu leurs symétries à cause de la cinématique et de la forme des solides. Toutefois, le comportement du contact est resté similaire à celui d’un contact elliptique, en dehors de certains cas limites. / This thesis is dedicated to the study of torus on plane contacts under various operating conditions. They can be found at the interface between the torus roller-end and the flange in roller bearings. The first challenge of this thesis is to deal with unusual mating geometries. The other challenge is the presence of a complex kinematic which operates in these contacts. In order to further develop the understanding of such a contact, a dual approach (experimental and numerical) is adopted. The Jérotrib test-rig enables a first study, by considering that the élastohydrodynamic torus on plane contact can be modelled by an elliptical equivalent contact. Thanks to a differential colorimetric interferometry method which was improved and adapted during the thesis, precise film thickness measurements are carried out under a rather wide range of operating conditions. A thermo-elastohydrodynamic numerical model is developed and validated by comparing its results to the ones of the test-rig. A numerical study on film forming is then proposed and the role of the contact ellipticity is investigated. The numerical model is improved in order to take into account the actual shape of the solids. A film thickness validation of the model is proposed, thanks to measurements performed on the Tribogyr test-rig. The operating conditions are very similar to the one encountered in actual bearings, and the mating solids have representative geometries: it is an actual torus-on-plane contact. It is demonstrated that the lubricant shearing is responsible for the solids temperature rise, which in its turn, reduces the film thickness. It appears mandatory to be able to predict this global warming of the bodies. It is also demonstrated that the pressure and film thickness distributions lose their symmetry because of the spinning kinematic and the solids shape. However, the behaviour of the torus-on-plane contact appears very similar to the one of an elliptical equivalent contact, apart from some limit cases.

La diminution de l’impact environnemental des moteurs électriques conduit à des améliorations du rendement et de la compacité. Puisque la taille de la transmission est réduite, des problèmes d’échauffements peuvent apparaître car la surface d’échange avec l’environnement est fortement réduite, diminuant les possibilités de refroidissement. L’ensemble du système doit être caractérisé en termes de pertes de puissance et de capacité de transfert thermique pour déterminer si la température maximale en fonctionnement correspond aux exigences de fiabilité tout au long de la durée de vie du système. Cette thèse propose un modèle thermique d’une transmission avec des engrenages spiraux coniques, lubrifiée par bain d’huile. Le modèle thermique utilise la méthode des réseaux thermiques. Par conséquent, une estimation précise des pertes et des transferts thermiques entre les différents éléments est nécessaire. Une revue des méthodes existantes est présentée, indiquant que les pertes par barbotage des engrenages spiraux coniques ne sont pas totalement maîtrisées, entraînant la conduite d’une série d’essais sur un banc d’essai dédié. Cette étude se conclue sur l’application du modèle à la transmission ainsi que la proposition de points d’améliorations. / Environmental concern leads to improvements in both efficiency and power density of electrical motors. For compact transmissions, thermal issues can become critical: the surface available for thermal exchange with the environment is reduced, minimizing the possibility to cool down the elements. To ensure the system reliability, the operating temperature for each operating condition must not exceed the required maximum oil temperature. In order to evaluate its thermal behavior, the whole system should be therefore evaluated in terms of power losses and heat transfer. This PhD proposes a thermal model of an existing oil bath lubricated gear unit with spiral bevel gears. The thermal model uses the thermal network method, which needs an accurate prediction of power losses of the different elements and heat transfer between them. A review of the existing methods is presented. Since the churning losses of spiral bevel gears has never been fully characterized, specific tests have been done on a dedicated test bench. The last part of this work suggests possible improvements on the gearbox design.

La corrosion sous contrainte (CSC) est un phénomène synergique d'endommagement qui résulte d'un processus de corrosion (dissolution, adsorption) et d'une rupture mécanique (fissuration). Les mécanismes de couplage mécano-électrochimique en pointe de fissure nécessaires à la compréhension du phénomène sont encore mal connus puisqu'ils dépendent du système d'étude (métal/milieu agressif) et font intervenir de nombreux facteurs mécaniques et électrochimiques. Dans cette thèse, nous nous proposons d'étudier les interactions réciproques entre la dissolution et l'état de contrainte mécanique en pointe de fissure (facteurs d'intensité des contraintes) pour le cas du Zircaloy-4 en milieu aqueux halogéné. Les éprouvettes sont d’abord pré-fissurées par fatigue à l’air avec la technique du Load-Shedding, ce qui permet alors de maîtriser le facteur d’intensité des contraintes résiduel de pré-fissuration par fatigue. Ensuite, une pré-oxydation thermique est réalisée pour produire une couche de film passif sur les surfaces de l’éprouvette. Les réactions électrochimiques sont alors concentrées en pointe de fissure qui induit également une concentration des effets mécaniques. Des techniques de corrélation d'images sont développées dans le but d'identifier les facteurs d'intensité des contraintes et d'estimer la longueur de fissure en temps réel. Ainsi des essais originaux de CSC, permettant d'imposer les facteurs d'intensité en pointe de fissure, sont conduits et les résultats montrent alors clairement les effets synergiques entre l'évolution des mesures mécaniques et des courants de dissolution. D'autre part, l'existence d'un facteur d'intensité des contraintes seuil de propagation en CSC, présentant une forte dépendance à l'histoire de chargement et à la variation de l'état de contrainte locale en pointe de fissure (variation des facteurs d'intensité des contraintes), est mise en évidence. Ceci montre alors que la plasticité doit évoluer pour que la dissolution se produise. Ainsi, le temps caractéristique d'arriver de nouvelles dislocations en pointe de fissure ne doit pas dépasser la durée caractéristique des réactions de passivation. Enfin sur la base des résultats expérimentaux, un modèle de loi de propagation phénoménologique de fissure en CSC est proposé et ses paramètres sont identifiés et validés à partir de différents essais de CSC. / Stress corrosion cracking (SCC) is a damage phenomenon which results from the synergy between corrosion process (dissolution, adsorption) and mechanical fracture (crack propagation). Although this phenomenon is well known, its modelling is still a challenging issue, especially concerning mechano-electrochemical coupling mechanisms at crack tip, because it depends on model system (metal/aggressive media) and large number of mechanical and electrochemical factors. In this thesis, mutual interactions between dissolution and the stress state around the crack tip (stress intensity factor) are studied in the case of Zircaloy-4 in aqueous halide solution. Samples are first pre-cracked in air by using fatigue load-shedding procedure to control the stress intensity factor. Then, pre-oxidization is used to produce a thin protective passive layer on their surface. The electro-chemical reactions are thus concentrated at the crack tip which also induces a concentration of the mechanical effect. During the test, digital images of the sample surface are acquired. Digital Image Correlation is performed a posteriori in order to obtain the evolution of the crack length and the stress intensity factors. Further, a specific procedure is developed in order to perform the DIC analysis while the test is running. This allows to control the load so that a given value of the stress intensity factor is prescribed. With this innovative experimental technique, we perform experimental tests that allow to discriminate the effects between different stress corrosion cracking mechanisms. It is suggested that once a critical anodic polarization is exceeded, the crack growth rate depends on the stress intensity factor but also on its time derivative. Indeed, a threshold effect is obtained on the stress intensity factor, meaning that plasticity must increase for the dissolution reaction to occur, but also on its rate meaning that time for plasticity to produce new dislocations must not exceed the characteristic duration of the passivation reaction. A phenomenological crack propagation model is proposed and its parameters are then identified and validated from the experimental measurements during SCC tests.

Cette thèse concerne l’étude des contacts pivotants rencontrés à la conjonction collet-rouleau, entre la bague d’un roulement et le flanc d’un rouleau. L’orientation principale de l’écoulement du lubrifiant peut changer lorsque le contact est mis à l’oblique. Cette cinématique complexe influe sur le comportement du contact. Elle est donc étudiée par une approche duale, expérimentale et numérique. Le banc d’essai Tribogyr permet l’expérimentation du contact à l’échelle 1:1. Une méthode pour la mesure de l’épaisseur du film lubrifiant par interférométrie optique en lumière blanche a été développée sur le banc d’essai et rend possible la mesure d’épaisseurs entre 0 et 800 nm, avec une résolution de quelques nanomètres. La mesure des efforts dans le sens de l’écoulement montre des similitudes avec les contacts de type roulement-glissement bien que le coefficient de frottement soit globalement plus faible. Les efforts transverses ont des valeurs du même ordre de grandeur que les efforts longitudinaux. Ils sont dus au cisaillement transverse induit par le pivotement. Un modèle numérique a été développé dans le but de simuler ces contacts pivotants. Le modèle inclut le calcul des températures et la rhéologie non Newtonienne du lubrifiant dans une stratégie de résolution par éléments finis, totalement couplée. La validation avec des résultats expérimentaux issus de Tribogyr, en épaisseur de film et en frottement, a été effectuée. Il est montré que l’épaisseur de film chute lorsque le pivotement et l’obliquité cisaillent le fluide, entraînant des effets rhéo-fluidifiants et thermo-fluidifiants. En cas de fort pivotement, le lubrifiant sortant peut être réinjecté à nouveau vers l’intérieur du contact et les transferts de chaleur entre lubrifiant et solides en sont fortement perturbés. Une forte obliquité entraîne à la fois la formation d’une augmentation locale de l’épaisseur de film et peut aussi provoquer la sous-alimentation du contact. Plusieurs campagnes expérimentales couplées à l’utilisation intensive du modèle numérique ont permis de comprendre les phénomènes physiques entrant en jeu et de prévoir l’efficacité, en terme de pertes de puissance, d’un contact pivotant. / This thesis is devoted to the study of spinning contacts located in bearing between the roller-end and the ring flange. The main direction of the lubricant flow may change when the contact is subjected to skew. This complex kinematics influences the contact behaviour. A dual experimental-numerical approach is proposed to study this problem. The Tribogyr test-rig allows the experimentation of the contact at the 1:1 scale. A film thickness measurement method, based on white light interferometry, was developed on Tribogyr. This method allows the measurement of film thickness between 0 and 800 nm with an accuracy of a few nanometres. The measurement of forces in the main flow direction shows similarities with classical rolling-sliding contacts. However, the friction coefficient is globally lower as soon as spin is involved. Transverse forces are of the same order of magnitude as the longitudinal forces. This is due to transverse shearing caused by the spin. A numerical model has been developed for the simulation of these spinning contacts. The finite element model, which is based on a fully-coupled solving strategy, takes into account the temperature calculation and the lubricant non-Newtonian rheology. Its validation with Tribogyr experimental results in terms of film thickness and friction has been conducted. Spin and skew effects induce high shear-thinning and thermal-thinning of the lubricant that lead to a decrease of the film thickness. Under high spinning condition, the lubricant exiting the contact may be re-injected to the contact inlet. Consequently, the heat transfers between the lubricant and the solids in contact are modified. In contact subjected to high skew, a local increase (dimple) of the film thickness may occur. Important skew may also lead to starvation conditions. Many experimental campaigns, coupled with an intensive use of the numerical model, allowed to understand the physical phenomena involved as well as to predict the efficiency, in terms of power losses, of the spinning contacts.

En tribologie, la modélisation numérique est aujourd'hui un outil indispensable pour étudier un contact afin de pallier les limites expérimentales. Pour comprendre de mieux en mieux les phénomènes mis en jeu, les modèles ne se situent plus à une seule échelle, mais en font intervenir plusieurs, rendant plus que jamais le concept de triplet tribologique incontournable. Travaillant avec cette philosophie et en se basant sur l'approche Non Smooth Contact Dynamics, dont nous rappelons les grandes lignes, nous proposons de franchir deux cas: proposer des modèles offrant des résultats quantitatifs et mettre en place les premières pièces d'une homogénéisation au niveau du contact (VER). Dans le premier cas, l'étude du couplage éléments finis/éléments discrets au sein d'une même simulation a pour but de proposer des modèles plus "réalistes". Même si l’interface utilisée est déjà présente au coeur du contact et ne va pas évoluer, elle permet de mettre en évidence l’utilisation d’outil de mesure permettant de lier le mouvement des particules aux instabilités dynamiques et permet d’avoir des résultats qualitatifs mais aussi quantitatifs puisque la comparaison avec les taux de contraintes expérimentaux sont en très bonne adéquation. Dans le second cas, le VER sous sollicitations tribologiques est étudié afin d'étendre les techniques d'homogénéisation aux problèmes de contact afin de s'affranchir de la description des interfaces aux grandes échelles en trouvant un moyen d'homogénéiser le comportement hétérogène de l'interface et de le faire dialoguer avec le comportement continu des corps en contact en faisant remonter, dans un sens, des grandeurs moyennées à l'échelle microscopique à l'échelle macroscopique des premiers corps et dans l'autre sens, se servir des données locales à l'échelle macroscopique comme conditions limites à l'échelle microscopique. / In tribology, the numerical modeling has become an indispensable tool for studying a contact to overcome the experimental limitations. To have a better understanding of the phenomena involved, the models are no longer located at a single scale, but involve several ones, more than ever, making the concept of tribological triplet as a unavoidable concept. Working with this philosophy and approach based on the Non Smooth Contact Dynamics framework, which we remind some outlines, we propose to cross two steps~: model that can offer quantitative results and that implement the first ingredient to perform a homogenization at a contact level. In the first case, the study of coupling finite elements/discrete elements within the same simulation aims to propose models that are more "realistic". Even if the interface is already present in the contact and not going to evolves, it can highlight the use of measurement tool of spot particles via dynamic instabilities and allows to have not only qualitative results but also quantitative ones since the comparison with the experimental strain rates are in very good agreement. In the second case, the study of VER in tribological charges is performed to extend the homogenization techniques to contact problems in order to overcome the interface description on large scales by finding a way to homogenize the heterogeneous behavior of the interface and make a dialogue with the continue behavior of bodies in contact by send up, in a sense, average values of the microscopic scale to the macroscopic scale and in the other sense, use local data of the macroscopic scale as boundary conditions at the microscopic scale.

Les aubes de turboréacteurs sont aujourd’hui sujets à de nouvelles technologiques novatrices : les matériaux composites tissées 3D. Ces modifications matériaux permettent des réductions de masse significatives, et conduisent à des performances énergétiques et environnementales accrues. Ces structures tissées possèdent généralement des architectures complexes et requièrent une bonne caractérisation mécanique. Des outils de modélisation textile et d’homogénéisation ont été développés, la géométrie et les propriétés physiques de l’architecture tissée y sont analysées. Toutefois, il est important de prévoir les endommagements engendrés dans ces matériaux au niveau des zones de contact soumises à du fretting. En effet, les pieds d’aubes de turboréacteurs font face à des sollicitations de type fretting. Il en résulte deux types d’endommagements : l’amorçage et propagation de fissures et l’usure des surfaces en contact. Quantifier ces endommagements requiert une analyse fine du contact à partir d’un outil de calcul à la fois robuste et rapide. Pour mener à bien cet objectif, un code semi-analytique permettant aujourd’hui de traiter des problèmes hétérogènes de contacts élasto-plastiques est développé. La méthode d’inclusion équivalente proposée au sens d’Eshelby est utilisée dans le solveur de contact pour décrire l’effet des hétérogénéités (cavités, inclusions, fibres ou torons). L’un des corps en contact peut contenir de multiples hétérogénéités de formes parallélépipédiques, ellipsoïdales, et ses formes dégénérées (ellipsoïde oblate, ellipsoïde prolate, sphère, cylindre, disque plat,...). Cette méthode est modifiée et améliorée afin de prendre en compte les influences mutuelles entre les inclusions hétérogènes voisines, et la décohésion à l’interface hétérogénéité/matrice. Un premier couplage avec un modèle de contact aube/disque Eléments Finis permet de prendre en compte les effets de structures. Cette méthode est analogue à un zoom structural lorsqu’un maillage fin est nécessaire dans des zones de contact à fort gradient de contraintes soumises à un endommagement de type fretting. Un second couplage avec le logiciel WiseTex permet de décrire la géométrie réelle du tissage du composite, les propriétés matériaux des fibres et de la matrice. Une segmentation fine du modèle permet de discrétiser les mèches en de multiples hétérogénéités ellipsoïdales équivalentes et d’appliquer la méthode d’inclusion équivalente modifiée. A partir de cette mise en donnée d’un matériau composite tissé revêtu, des calculs de contact tridimensionnel multiéchelle en situation de fretting sont réalisés avec succès. / New baseline turbofan engines feature advanced blade technology made of 3D woven composites for a significant reduction in weight and an increase in energy and environmental performance. Woven structures generally have complex architectures wich require high level of mechanical model. Modeling textile and homogenization tools have been developed to provide precisely the geometry and physical properties of the woven architecture. However, it is important to predict the damage mode under contact zone which are subjected to fretting. Indeed, the blade roots of turbofan engines are damaged by fretting phenomenon. Fretting modes generate two kinds of damage : (i) the initiation and propagation cracks, (ii) wear on contact surfaces. Quantify these two kinds of damage requires a detailed contact analysis from a robust and fast contact solver. In order to reach this goal, a semi-analytical solver allows to solve heterogeneous elasto-plastic contact problems. Equivalent inclusion method in the sense of Eshelby allows to describe accurately the effect of inhomogeneities (cavities, inclusions, fibers or strands). Only one of contact bodies contains multiple heterogeneous of cuboïdal and ellipsoidal shapes, and their degenerated forms (oblate spheroid, prolate ellipsoid, sphere, cylinder, flat disk,...). This method is modified and improved in order to take into account the mutual influence between neighboring heterogeneous inclusions and decohesion at the interface heterogeneity/matrix. A first coupling with a finite element model of blade/disk contact allows to take into account the effects of structure. This method is analogous to a structural zoom when structural fine mesh is necessary within areas of high contact stress gradient subjected to fretting. A second coupling with the software WiseTex allows to describe the geometry of the actual weaving of the composite, the material properties of the fibers and the matrix. A fine segmentation of the numerical model allows to discretize reinforcements to multiple equivalent ellipsoidal heterogeneities and to apply the modified equivalent inclusion method. Three-dimensional simulation of frictional contact model are successfully completed from these data of a coating woven composite.