Les problèmes de mécanique du contact sont des problèmes multi-échelles mettant en jeux de nombreux phénomènes physiques. Les premières études concernant ce domaine datent de la fin du XIXème siècle et les développements majeurs ont été réalisés au cours du XXème siècle en parallèle du besoin croissant des ingénieurs de prévoir le comportement des matériaux sous sollicitations tribologiques. L'évolution des besoins industriels et les avancées technologiques réalisées dans le domaine du numérique conduisent à réaliser des simulations tridimensionnelles ayant pour objectif la prédiction du comportement de pièces sous sollicitations thermo-mécaniques transitoires, pour, soit alléger les structures, augmenter le niveau de sollicitations, étendre la durée de vie... Ces simulations se révèlent très souvent coûteuses en termes de temps de calcul et d'espace mémoire et nécessitent par conséquent l'utilisation de super calculateurs. Dans ce contexte, cette thèse propose un modèle innovant basé sur les techniques multigrilles avec raffinement local afin de réaliser ces simulations pour des coûts numériques faibles. Ce modèle est basé sur les équations de Lamé généralisées et l'équation de la chaleur de Fourier discrétisée à l'aide des différences finies. Le système linéaire obtenu est résolu à l'aide de la méthode itérative de Gauss-Seidel couplée avec les techniques multigrilles. Ces techniques permettent d'accélérer la convergence d'un problème en utilisant plusieurs grilles et des opérateurs de transfert. Afin de garantir une convergence optimale et de minimiser la taille mémoire dans le cas de variations de propriétés importantes, des techniques numériques de localisation et d'optimisation ont été mises en place. Les applications visées ici sont centrées sur l'utilisation de revêtements ou de matériaux innovants pour permettre les gains attendus. Des validations du modèle ont été effectuées en comparant nos résultats avec ceux issus de la littérature. Des études paramétriques ont permis d'étudier l'influence de l'épaisseur du revêtement, de la valeur du module de Young mais aussi d'une couche à gradient de propriété sur le champ de contrainte et la tenue du système revêtement/substrat sous sollicitation de contact. Des études similaires ont été conduites sous sollicitations thermiques. L'intérêt porté aux variations de propriétés des matériaux selon toutes les directions de l'espace a conduit à étudier l'effet de la microstructure, qui est constituée de grains ayant chacun leurs propriétés propres, sur les champs de contraintes. La mise en évidence de cet effet est explicitement montrée au travers de calculs de durée de vie utilisant des descriptions statistiques de type Weibull. La dispersion observée sur les résultats est conforme aux observations expérimentales.

Dans ce travail de thèse, on s’intéresse au traitement endovasculaire de l’anévrisme de l’aorte abdominale (EVAR). Cette technique mini-invasive couramment utilisée et connaît une croissance importante depuis 10 ans. Elle repose sur l’exclusion du sac anévrismal par le déploiement au niveau de l’anévrisme d’une ou plusieurs endoprothèses introduites par voies fémorales. Au cours de l’intervention, plusieurs types d’outils de rigidité variable sont introduits pour permettre la navigation de l’endoprothèse. La structure vasculaire subit alors des déformations importantes. Ces déformations sont en général sans incidence sur le bon déroulement de l’intervention. Cependant dans certains cas, notamment pour les patients présentant des anatomies défavorables (fortes tortuosités ou angulations, important degré de calcification, longueur importante des artères iliaques communes et externes) les déformations produites par l’insertion des guides rigides peuvent avoir des conséquences sur le déroulement de l’intervention. Actuellement leur anticipation repose principalement sur l’expérience du chirurgien. La simulation mécanique semble être un outil adapté pour fournir des indicateurs plus objectifs et utiles au praticien lors du planning de son intervention : cette pratique permettrait en guidant et sécurisant le geste chirurgical de diminuer potentiellement les risques de complications peropératoires et postopératoires. La première étape du travail a consisté à développer un modèle mécanique de la structure aorto-iliaque et une méthode de simulation permettant de répondre au problème mécanique posé. Ce modèle a été paramétré de façon patient-spécifique à partir des données préopératoires disponibles. Puis la deuxième étape du travail a consisté à valider la modélisation développée en la confrontant à des données peropératoires réelles obtenus sur 28 cas de patients opérés au CHU Rennes. L’ensemble des méthodes développées à enfin été intégré à un module de démonstration du logiciel EndoSize® (Therenva, France). / Endovascular Aneurysm Repair (EVAR) is a mini-invasive technique that is commonly used to treat Abdominal Aortic Aneurysms (AAA). It relies on the exclusion of the aneurysm sac by introducing one or more stent-grafts through the femoral arteries and deploying them inside the aneurysm. During the procedure, several tools of varying stiffness are introduced to enable the delivery of the stent graft to its deployment site. During this process, the vascular structure undergoes major deformations. Usually, these have no consequence on the smooth progress of the procedure. However, in some instances, particularly when the patient presents an unfavorable anatomical profile (major tortuousness or angulation, deep calcification, long length of the common and external iliac arteries), the deformation caused by the insertion of stiff guidewires can have major consequences. Today, their prediction relies mainly on the surgeon’s experience. Numerical simulation appears to be an appropriate tool to give the practitioner more objective and more useful indicators when planning the procedure: guiding the surgical act and making it safer using such an approach would potentially reduce the risks of intraoperative and postoperative complications. In the first step of the work, we developed a mechanical model of the aorto-iliac vascular structure and a simulation methodology to answer the mechanical problem. This patient-specific model has been parametrized based on available preoperative data. Then the second step of the work consisted in the validation of this model by confronting the simulation results to real intraoperative 3D data that were collected on 28 cases of patients operated at the University Hospital of Rennes. All the methods that were developed during this PhD were integrated in demonstration module of EndoSize® software (Therenva, France).

La durée de vie en fonctionnement des pièces mécaniques peut être augmentée par la présence de contraintes résiduelles de compression. Inhérentes à la plupart des procédés de fabrication, contraintes résiduelles jouent un rôle sur la tenue en service des éléments mécaniques. La connaissance et la maîtrise de ces contraintes résiduelles et des procédés associés sont donc fondamentales. La mise en place de méthodes numériques performantes pour prédire ces contraintes résiduelles permettra à terme d'éviter de nombreux et coûteux essais et d'étudier l'influence des principaux paramètres. Cette thèse présente le développement et l'application de méthodes de calcul semi analytiques (SA) à la modélisation de procédés mécaniques de mise en compression des surfaces. Les méthodes SA, initialement développées pour la simulation des contacts élasto-plastiques, ont l'avantage de temps de calcul très significativement réduits par rapport aux méthodes numériques conventionnelles. Cette méthode est d'abord utilisée pour la simulation d'un procédé connu sous le nom de galetage, avec un chargement de type contact roulant. Enfin elle est utilisée pour la simulation d'impacts uniques puis répétés, et les développements associés présentés. Le contact roulant entre deux massifs élasto-plastiques, sans considération du frottement, est d'abord étudié. L'influence de la plasticité et du type d'écrouissage (isotrope ou cinématique), de la géométrie des massifs en contact (sphéres ou ellipsoides) et du type de chargement (indentation ou roulement) sur la pression de contact et les déformations plastique sont analysés. La simulation d'impacts est ensuite abordée. La méthode développée est tout d'abord validée numériquement, puis confrontée à l'expérience. Trois matériaux ont été plus particulièrement étudiés : 316L, AA 7010, et l'Inconel 600. Les dimensions des impacts ainsi que les déformations générées, mesurées par stéréo corrélation sont utilisées pour valider expérimentalement la méthode. Le procédé de billage ultrasonore a été tout particulièrement étudié. Dans un premier temps nous nous sommes attachés à la décrire la cinématique des billes dans une chambre fermée, mises et maintenues en mouvement par une sonotrode. L'utilisation de formules analytiques pour l'estimation des coefficients de restitution, lors des nombreux chocs entre les billes ou avec les parois de la chambre, a permis d'affiner le calcul des vitesses moyennes d'impact en fonction des paramètres du procédé. La méthode SA est ensuite utilisée pour déterminer le champ de déformations plastiques induits par les impacts. Enfin une méthode de projection est proposée pour permettre in fine la prédiction des contraintes résiduelles pour des structures fines ou épaisses. / The life time of mechanical components can be increased by the presence of compressive residual stresses. Inherent to most production processes, residual stresses play a critical role on the mechanical parts behaviour. The knowledge and mastering of residual stresses and linked processes are thus fundamental. The development of efficient numerical methods to predict these residual stresses will allow to save costly experiments and to study the influence of the main parameters. This PhD presents the development and application of semi analytical methods (SAM) to the modelling of mechanical processes of compressive residual stress generation. The SAMs, initially developed for the simulation of elasto-plastic contacts, have the advantage of significantly low computation time compared to classical numerical methods. This method is first used to simulate the low plasticity burnishing process, with a rolling loading. Then, it is used for the simulation of impacts, first unique and then repeated. The frictionless rolling contact between two elasto-plastic bodies is first studied. The influence of plasticity, of the hardening model (isotropic or kinematic), of the geometry of the bodies in contact (spheres or ellipsoids) and of the loading type (indentation or rolling) on the contact pressure and plastic strains are analysed. Impacts simulation is then addressed. The developed method is first validated numerically then confronted to experimentations. Three materials have been particularly studied: 316L, AA 7010 and Inconel 600. The impacts dimensions and the generated strains, measured by digital image correlation, are used to validate experimentally the method. The ultrasonic shotpeening process has been specifically studied. The description of the kinematics of the shots put in movement by a sonotrode in a closed peening chamber has first been studied. The use of analytical formulae for the estimation of the coefficients of restitution, during the numerous impacts between shots and with the chamber’s walls, allowed refining the calculation of the average impact velocity as a function of the process parameters. The SAM is the used to determine the plastic strain field induced by the impacts. At last a projection method is proposed to finally determine the residual stress field in thick or thin structures.

Le récent développement des menaces terroristes renforce l'effort de recherche du CEA et d'EDF pour la protection des citoyens et des installations. De nombreux scénarios doivent être envisagés comme, par exemple, la chute d'un avion de ligne sur une structure de génie civil. La dispersion du carburant dans la structure, son embrasement sous forme de boule de feu et les effets thermiques associés sont des éléments essentiels du problème. L'utilisation de modèles numériques est indispensable car des expériences seraient difficiles à mettre en œuvre, coûteuses et dangereuses. Le problème type que l'on cherche à modéliser est donc l'impact d'un réservoir rempli de fluide, sa déchirure et la dispersion de son contenu. C'est un problème complexe qui fait intervenir une structure mince avec un comportement fortement non-linéaire allant jusqu'à rupture, un fluide dont la surface libre peut varier drastiquement et des interactions fluide-structure non permanentes. L'utilisation des méthodes numériques traditionnelles pour résoudre ce problème semble difficile, essentiellement parce qu'elles reposent sur un maillage. Cela complique la gestion des grandes déformations, la modélisation des interfaces variables et l'introduction de discontinuités telles que les fissures. Afin de s'affranchir de ces problèmes, la méthode sans maillage SPH (\og Smoothed Particle Hydrodynamics \fg) a été utilisée pour modéliser le fluide et la structure. Ce travail, inscrit dans la continuité de recherches précédentes, a permis d'étendre un modèle de coque SPH à la modélisation des ruptures. Un algorithme de gestion des interactions fluide-structure a également été adapté à la topologie particulière des coques. Afin de réduire les coûts de calcul importants liés à ce modèle, un couplage avec la méthode des éléments finis a également été élaboré. Il permet de n'utiliser les SPH que dans les zones d'intérêt où la rupture est attendue. Finalement, des essais réalisés par l'ONERA sont étudiés pour valider la méthode. Ces travaux ont permis de doter le logiciel de dynamique rapide Europlexus d'un outil original et efficace pour la simulation des impacts de structures minces en interaction avec un fluide. Un calcul démonstratif montre enfin la pertinence de l'approche et sa mise en œuvre dans un cadre industriel. / The recent development of terrorist threats increases the research effort of the french Atomic Energy Commission (CEA) and the French Electricity company (EDF) for the protection of citizens and facilities. Many scenarios should be considered as, for example, the fall of an airliner on a civil engineering structure. The dispersion of fuel in the structure, the formation of a fireball and associated thermal effects are essential elements of the problem. The use of numerical models is essential because experiences would be difficult to organize, costly and dangerous. The typical problem that we want to model is the impact of a tank filled with fluid, its rupture and the dispersion of its contents. It is a complex problem which involves a thin structure with a highly non-linear behavior up to rupture, a fluid with a free surface that can vary drastically and non permanent fluid-structure interactions. The use of traditional numerical methods to solve this kind of problems is difficult, mainly because they rely on a mesh. This complicates the management of large deformations, the modeling of moving interfaces and the introduction of discontinuities such as cracks. To overcome these problems, the meshfree method SPH (Smoothed Particle Hydrodynamics) was used to model both the fluid and the structure. This work, which is a continuation of previous research, has extended a model of SPH shell to the modeling of ruptures. An algorithm for managing fluid-structure interactions has also been adapted to the particular topology of shells. To reduce the important computational costs associated with this model, a coupling with the finite element method was also developed. It allows the use of SPH in areas of interest where the rupture is expected. Finally, tests performed by the french Aerospace Lab (ONERA) are studied to validate the method. This work helped to provide fast dynamic software Europlexus an original and effective tool for the simulation of the impact of thin structures interacting with fluid. A demonstrative calculation finally shows the relevance of the approach and its use within an industrial framework.

Grâce à leur géométrie spécifique, faite d’une superposition de couches unidirectionnelles de fibres continues, liées par des coutures plus ou moins complexes, les renforts NCF (Non-Crimp Fabric) possèdent une combinaison de propriétés intéressantes. Ces propriétés permettent l'utilisation de procédés de production tels que le procédé RTM (Resin Transfer Moulding) et permettent l’utilisation des caractéristiques de rigidité élevée des fibres dans la pièce finale. Ainsi, ces renforts trouvent des applications dans différents secteurs de l'industrie des composites et notamment en aéronautique. Les couches, formées de fibres parallèles juxtaposées, sont liées par différents types de coutures. Ce lien conduit à des propriétés différentes de celles observées sur un renfort tissé où les mèches de chaîne et de trame sont liées par entrelacement. Par conséquent, la cinématique de déformation des renforts NCF secs lors de la mise en forme peut être différente. Dans ce travail de recherche, nous menons d'abord une étude expérimentale sur les mécanismes de déformation se produisant lors de tests élémentaires et lors d’une expérience d’emboutissage. Plusieurs renforts NCF et un renfort interlock sont étudiés. Ensuite, nous introduisons un nouveau modèle pour la simulation de la mise en forme du renfort de composite NCF sec. Cette approche implique l’utilisation d’éléments finis lagrangien semi-discrets pour représenter les couches et d’éléments de barres pour représenter les coutures. Ce modèle nécessite un algorithme de contact spécifique pour gérer la relation entre les éléments finis. Une loi de frottement de Coulomb avec seuil de glissement est utilisée. Enfin, le modèle est validé par l’intermédiaire de l’essai de bias test puis la mise en forme sur un hémisphère est simulée et comparée à l'expérience. / Due to their specific geometry, made of overlapping layers of unidirectional plies of continuous fibres, linked by more or less complex stitching, non-crimp fabric (NCF) reinforcements have advantageous combination of properties. These properties enable the use of process in production such as RTM process and allow the full use of stiffness characteristics of the fibres in the final part. Thus, these reinforcements find applications in different areas of composites industry. The plies of juxtaposed parallel fibres are bound by a tricot or a chain stitching with a lead. This link leads to different properties of the NCF compared to those of a woven reinforcement where the warp and weft yarns are linked by the interweaving. Therefore, deformation kinematics of dry NCF reinforcements might be different. In this research work, we first carry out an experimental study on deformation mechanisms occurring during elementary tests and a drawing experiment. Several non-crimp fabrics and an interlock fabric are studied. Then, we introduce a new model for the simulation of the dry composite reinforcement preforming. This approach involves a semi discrete Lagrangian finite element to represent the plies and a truss element to represent the stitching. This model requires a specific contact algorithm to handle the relation between the finite elements. A Coulomb law with a friction threshold is used. Finally, the model is validated using the bias extension test and the preforming on a hemisphere is simulated and compared to the experiment.

Dans le cadre de la maintenance conditionnelle, l'analyse vibratoire temporelle repose sur 50 ans d'expérience et jouit, de fait, d'une représentativité importante dans le secteur éolien. Cette analyse, appliquée sous des conditions de fonctionnement non stationnaires, fait néanmoins preuve d'un manque de précision dans la détection de défauts mécaniques. Afin de synchroniser l'échantillonnage avec les défauts des éléments tournants de la ligne d'arbre, il est alors nécessaire de s'appuyer sur une discrétisation angulaire du vecteur d'information. L'objet de cette étude est le développement d'un outil de surveillance disposant de cette qualité et s'appuyant sur les variations que présente la vitesse instantanée de la machine tournante pour évaluer l'état des éléments tournants. La vitesse angulaire instantanée, associée à l'analyse spectrale, est assujettie à des perturbations dont les origines sont identifiées et les influences quantifiées. Ainsi, il est envisageable de dimensionner un système de mesure en fonction des caractéristiques de la chaîne cinématique sous surveillance. Ce document propose également plusieurs outils de traitement originaux qui permettent, sous des conditions de fonctionnement non stationnaires, d'améliorer suffisamment l'observation pour que l'utilisation d'un seul capteur soit envisageable pour assurer la surveillance de la ligne d'arbre dans sa totalité. Ces outils mettent en exergue la dualité temporelle angulaire de la surveillance de machines tournantes et tirent parti des conditions de fonctionnement non stationnaires par la connaissance de la vitesse instantanée. Les procédés d'acquisition combinés aux outils de traitement sont mis à l'épreuve sur plusieurs banc d'essais. Sensible au degrés de désalignement, au type de défaut d'engrènement et à la gravité d'un écaillage, la vitesse angulaire instantanée est dorénavant une alternative confirmée à la surveillance vibratoire, même angulaire. Une vaste campagne de mesures a été lancée en parallèle et en continu sur une éolienne dans l'objectif ambitieux de mettre en place un procédé de pronostic. Les diverses conditions de fonctionnement sont prises en compte par une approche statistique permettant d'utiliser avantageusement leurs influences sur les mesures qui en sont issues. Ces travaux constituent une première étape dans la réalisation d'un modèle de surveillance multi-niveau adapté aux machines à fonctionnement variable dont les éoliennes sont un parfait exemple de complexité. / The overall background for the thesis is the need to develop new methods for monitoring machines subject to nonstationary operating conditions of speed and load, as typified by wind turbines. This has become particularly important because of the increasing prevalence of wind turbines for the sustainable supply of electrical power. Because of the cost of installing and maintaining wind turbines, often situated in remote areas and increasingly used off-shore, it has become imperative to use the most up-to-date monitoring methods to avoid unforeseen failure. For the same reason there has been a tendency to increase the size of individual units, but this gives problems for the transmissions used to convert the low speed input at a fraction of a Hz to the speed of conventional generators in the range 16-30 Hz or so. The power delivered by a turbine is proportional to the square of its diameter, but the input torque is proportional to the cube, and a large part of the cost of the gearbox is related to the torque. Anecdotally, this appears to have resulted in a situation where many transmissions are failing in a shorter time than their design life, and this means not only that their design must be improved, but also that they must be monitored more closely than many other machines, to detect and diagnose incipient faults at the earliest possible stage. The thesis seeks to help to redress the situation by proposing a number of innovative methods to improve the detection and diagnosis of faults in machines with variable speed and load, and especially wind turbine transmissions, not only to make the monitoring more efficient, but also more economical.

Lors de la réparation d’un défaut, le processus consiste à éliminer la matière dans la zone englobant le défaut, puis à recharger cet affouillement par un dépôt de cordons de soudure, afin de remettre en conformité la pièce en vérifiant les critères de conception. Pourtant la zone rechargée constitue en elle-même un site sensible où le niveau des contraintes résiduelles peut être préjudiciable de par son importance, notamment localement, s’il porte des contraintes de traction. Les risques de récidive d’endommagement pourraient être réduits par certains procédés de parachèvement. Ces procédés permettent de générer une zone de compression en peau externe et en sous couche de la pièce traitée. Par conséquent, la durée de vie d’une pièce ainsi parachevée pourrait être augmentée. Parmi les procédés existants, nous nous intéresserons dans cette étude à la mise en compression par billage ultrasonore. L’objectif de cette contribution est de prédire le niveau des contraintes résiduelles induites par le billage sur une pièce rechargée par soudage multi-passes. Pour mener à bien cet objectif, une étude expérimentale du billage a été réalisée sur des éprouvettes minces et épaisses afin d’identifier les paramètres influents et de disposer d’une base de résultats expérimentaux – comme par exemple la répartition quantitative des contraintes résiduelles dans la matière traitée – à confronter à la simulation. L’originalité de ce travail est de mettre en œuvre, en termes d’évaluation, une modélisation tridimensionnelle représentant correctement les phénomènes locaux et basée sur un modèle semi-analytique. Ce modèle numérique présente l’avantage de temps de calcul considérablement réduits par rapport aux méthodes éléments finis plus conventionnelles, sans pour autant sacrifier à la qualité des résultats. Pour cette application, ce modèle semi-analytique a été validé par des benchmarks aux éléments finis dans le cas de mono et multi impacts. Puis cette approche prédictive a été validée par l’expérience. Enfin, une simulation numérique du soudage a été réalisée afin de disposer d’un champ de contraintes résiduelles de soudage. Les résultats des calculs, obtenus par chaînage de simulations numériques des procédés de soudage puis de billage, ont ensuite été confrontés aux grandeurs mesurées sur des éprouvettes initialement précontraintes par flexions alternées puis billées selon des paramètres du processus identifiés antérieurement. L’état précontraint (qu’il provienne de la thermique du soudage ou de l’écrouissage de flexion alternée) a été introduit en données d’entrée, en tant que pré-chargement, dans le modèle numérique simulant les impacts. / During a repair welding process, the first step is to remove the material in the defected area, then refill this chamfer by welding rods. Since, the nickel based alloy, such as IN600, is sensitive to the phenomenon of stress-corrosion cracking, it must be ensured that the level of residual stress after welding is lower than the triggering threshold of stress-corrosion cracking. At this stage, the ultrasonic shot peening (USP) process can be used as the final finishing step. The aim is to introduce compressive residual stresses on the surface of the welded area in order to prevent occurrence of cracking. The objective of this study is to predict the level of residual stresses, induced by the ultrasonic shot peening, on structure made of IN600 repaired by multipass welding. To achieve this objective, an experimental study of the ultrasonic shot peening was carried out. Experiments similar to the “Almen” test were performed on thin and thick specimens. Influential parameters were identified. Experimental results, such as the level of residual stresses, were obtained. The originality of this work is that a three-dimensional modeling of multiple impacts, based on a Semi-Analytical Method (SAM), was implemented. This method can reduce significantly the computing time compared to finite element method, without sacrificing the quality of results. This was first validated by benchmarks for both single and multiple impacts, then proved by experiments. Moreover, a numerical simulation of welding was performed. The field of welding residual stresses should be used as the input data for the impact model. In fact, in this study，the residual state，was generated by alternatif bending test. These specimens were then treated by USP process with the parameters identified previously. As the final step, numerical results in terms of distortions and residual stresses were compared with the experimental data in order to validate the chaining of welding and peening simulation.

La fatigue à nombre de cycles élevé (HCF) est un mode de défaillance courant et dangereux pour les aubages de turbomachines. Elle est induite par les efforts dynamiques élevés générés lors de résonances présentes dans la plage de fonctionnement de ces machines. Les dispositifs amortisseurs basés sur l’utilisation du frottement sec, tels que les nageoires ou les frotteurs sous-plateformes, permettent de réduire les amplitudes vibratoires, voire de repousser les fréquences de résonance hors des zones de fonctionnement. Cependant la conception de ces dispositifs reste encore largement basée sur l’empirisme et ils peuvent être la source d’un effet de désaccordage potentiellement nuisible. L’objectif ici est de développer des modélisations adaptées au traitement du problème de vibration des aubages en présence de frottement sec, ceci afin de mieux maîtriser les comportements physiques mis en jeu et donc, de mieux maîtriser leur processus de conception. Plusieurs modélisations numériques sont testées et confrontées à des résultats de référence. Une comparaison entre procédures de résolution temporelle et fréquentielle est menée et montre l’efficacité des méthodes fréquentielles. La méthode fréquentielle de la balance harmonique à plusieurs harmoniques est adaptée au problème et exploitée dans le cadre d’une étude énergétique. Cette étude conduit à une meilleure compréhension des phénomènes mis en jeu lors de l’aplatissement des pics en fonctionnement ; elle permet de démontrer que l’alternance des états de contact glissant et bloqué est à l’origine de cet aplatissement et non la dissipation d’énergie comme souvent avancé dans la littérature. Enfin, la méthode est exploitée pour décrire le comportement des disques aubés désaccordés. La méthode de Monte Carlo est utilisée pour obtenir les caractéristiques statistiques de la réponse forcée d’un système discret, en tenant compte des variations stochastiques des paramètres du contact notamment, la charge normale, la raideur du contact et le coefficient de frottement. Les résultats obtenus permettent de mieux comprendre les effets de la nature variable de ces paramètres fondamentaux sur la dynamique d’ensemble du système non linéaire. / Fatigue with high number of cycles (HCF) is a current and dangerous mode of failure for the blades of turbo shaft engines. It is induced by the high dynamic stresses generated at resonance in the operating range of these machines. The rubbing devices based on use of dry friction, such as shrouds or under-platform dampers, make it possible to reduce the vibratory amplitudes, to even push back the resonance frequencies out of the operation zones. However the design of these devices remains still largely based on empiricism and they can be the source of a potentially harmful effect of mistuning. The goal is to develop adapted modelling for the treatment of blades vibration problem in the presence of dry friction, this in order to better control the concerned physical behaviours and thus, to better control their process of design. Several numerical modelling are tested and confronted with reference results. A comparison between procedures of time and frequency domains resolution is carried out and shows the effectiveness of frequential methods. The frequential method of the Harmonic Balance including several harmonics is adapted to the problem and is used within an energy study. This study leads to a better comprehension of the phenomena of peaks flattening and it shows that the alternation of the states of slipping and sticking contact is the real cause of this flattening and not energy dissipation like often advanced in the literature. Lastly, the method is exploited to describe the behaviour of mistuned bladed disks. The Monte Carlo method is used to obtain the statistical characteristics of the forced response for a lumped system, by taking account of the stochastic variations of parameters in the contact, specifically the normal load, the stiffness of the contact and the coefficient of friction. Results obtained make it possible to better include/understand the effects of the variable nature of these fundamental parameters on dynamics of nonlinear system.

Un modèle de trains planétaires est proposé afin de tenir compte de l’influence d’erreurs de position et d’excentricités en lien avec d’éventuels montages ‘flottants’ sur le comportement dynamique d’une transmission. La formulation théorique repose sur le formalisme des torseurs de déplacements infinitésimaux pour simuler à la fois les erreurs géométriques et les degrés de liberté du modèle. Une des propriétés principales de cette approche est que la géométrie des engrènements et les excitations correspondantes sont couplées aux degrés de liberté, conduisant ainsi à des excitations complexes présentant des modulations d’amplitude et de phase. Les résultats de simulation sont comparés avec des mesures sur banc d’essai et un très bon accord est obtenu en terme de partage de charge entre les satellites, validant ainsi le modèle de contact développé. Enfin, des résultats d’études paramétriques portant sur le rôle de certaines erreurs ainsi que sur l’apport éventuel de solaire et/ou satellites flottants dans des applications grandes vitesses concluent ce travail de thèse. / A dynamic model of planetary gears is presented which accounts for planet position errors and eccentricities for either rigid mounts or floating members. The theoretical formulation relies on infinitesimal generalised displacement screws which can simulate both errors and deflections. A unique feature of this model is that mesh properties (geometry and excitations) are coupled with the degrees-of-freedom thus leading to complex frequency and amplitude modulated excitation sources. For a number of planetary gears, it is found that the simulated load sharing between the planets compare well with the experimental evidence thus validating the contact modelling strategy. Finally, the results of extensive parameter analyses are displayed which illustrate the role of certain errors along with the interest and drawbacks of floating sun-gears or planets in high-speed applications.

La détection et le diagnostic des défauts d’engrenages ont été traditionnellement basés sur l'analyse des signaux vibratoires et acoustiques. Mais, ces méthodes peuvent être coûteuses suite aux difficultés techniques de mise en œuvre de capteurs sur des pièces en rotation. Dans ces travaux de thèse, une méthode originale de détection de défauts locaux dans des engrenages entraînés par des moteurs électriques est proposée en se basant sur le suivi des courants statoriques. Le système électromécanique est simulé numériquement en combinant un modèle électrique dynamique de moteurs asynchrones (modèle de Kron et modèle des réseaux de perméances) et un modèle mécanique dynamique d’engrenages droits et hélicoïdaux incluant les arbres et les paliers. Certains défauts classiques du moteur (excentricité du rotor) et des défauts d'engrenages (écaillages sur les flancs actifs de dentures) sont considérés et une technique itérative de résolution est présentée qui permet de résoudre le problème électromécanique couplé. De nombreux résultats sont présentés qui illustrent la possibilité et la sensibilité de détection de défauts d’engrenages par mesure de courants électriques sur le moteur. Sur la base de ces simulations, il est conclu que l’analyse de la signature des courants du moteur apparaît comme une méthode suffisamment précise pour détecter et surveiller les défauts locaux d’engrenage comme les écaillages. / The fault diagnosis in geared transmissions has been traditionally based on vibration monitoring but the implementation of sensors can be expensive and, in some applications, present technical difficulties. In this thesis, an original way of detecting faulty gears driven by electrical motors is proposed based on the monitoring of the motor currents. The corresponding electromechanical coupled system is numerically simulated by combining i) an electrical dynamic model of asynchronous motors (Kron’s model and permeance network model) and ii) a mechanical dynamic model of spur and helical gears (3D lumped parameter model) including the shafts and bearings. Some classic motor defects (rotor eccentricity) and gear faults (tooth spalls) are considered and an iterative solution technique is presented which makes it possible to solve the coupled electromechanical problem. A number of results are presented, which illustrate the possibility of tooth fault detection by stator current measurements with regard to the position and dimensions of the defect. Based on these simulations, it is concluded that motor current signature analysis can be a viable tool for the detection of localized gear faults such as tooth pitting.