Dans le cadre de la réduction des émissions polluantes et de la consommation des véhicules à moteur thermique, les constructeurs cherchent à diminuer la cylindrée et la vitesse de rotation des moteurs de chaines cinématiques. Ces évolutions conduisent, du fait du principe même du moteur à pistons, à une augmentation significative des irrégularités de rotation de celui-ci. Depuis quelques années, le système à pendule est apparu dans les groupes moto-propulseurs automobiles. Il agit à la manière d’un batteur, accordé sur l’ordre d’allumage du moteur thermique, et permet ainsi une réduction des vibrations. Cependant, les fortes non-linéarités intrinsèques aux pendules provoquent un désaccord du système à grande amplitude synonyme de perte de performances. Cette thèse a pour but d’améliorer la compréhension et le comportement du système en interaction avec la chaîne de traction automobile. En renfort des traditionnelles méthodes d’intégrations temporelles, le système non linéaire est résolu par la méthode asymptotique numérique couplée à la méthode de l’équilibrage harmonique. Une méthode originale de continuation d’antirésonance est également proposée ainsi que des règles de conception issues de développements analytiques. La validation par l’expérience montre une amélioration significative des performances du système. / In the context of the reduction of polluting emissions and fuel consumption of thermal engines of vehicles, automotive manufacturers try to reduce cylinder capacity and engine speed of rotation. These evolutions lead to significant increase of irregularities of rotation. The so-called centrifugal pendulum vibration absorber is a recent solution of mitigation of torsional vibrations in automotive powertrains. It acts as a mass damper tuned on the firing order of the engine and allows reduction of vibrations. However, strong non-linearities intrinsic to pendular systems cause a detuning of the device at large amplitude of motion resulting in a loss of performances. This thesis aims to improve the understanding and the behavior of the system in interaction with an automotive driveline. In support of classic time integration procedures, the nonlinear system is solved through the asymptotic numerical method coupled to the harmonic balance method. In addition, an original continuation of antiresonance method is proposed as well as some design rules derived from analytical developments. Experimental validation shows a significant enhancement of performances of the system.

Les vibrations qui sont en général source de dérangement, d’usure et même destruction des machines et structures mécaniques doivent être contrôlées ou éliminées. Pour cette raison, la lutte contre les vibrations est devenue depuis des années un enjeu majeur pour les chercheurs de laboratoire et de développement dans l’industrie afin de développer des solutions efficaces contre ces problèmes. De nombreuses technologies ont donc été développées. Parmi ces technologies, les absorbeurs de vibration non linéaires présentent des performances importantes dans l’atténuation de vibration sur une large bande de fréquences. C’est dans ce contexte que cette thèse se focalise sur l’analyse et l’optimisation des absorbeurs de vibration non linéaires. L’objectif de cette thèse est d’analyser le comportement dynamique non linéaire des systèmes présentant des absorbeurs de vibration non linéaires. Pour cela, un modèle dynamique d’un système à deux degrés de liberté est développé mettant en équations le comportement non linéaire. La résolution des équations de mouvement est faite par la Méthode Asymptotique Numérique (MAN). La performance de cette méthode est montrée via une comparaison avec la méthode de Newton-Raphson. L’analyse des modes non linéaires du système ayant une non-linéarité cubique est faite par une formulation explicite des Fonctions de Réponse en Fréquence non linéaires (FRFs) et les Modes Normaux Non linéaires (MNNs). Un démonstrateur sur la base d’un système simple à deux degré de liberté est mis en place afin de recaler les modèles envisagés sur la base des résultats expérimentaux trouvés. / Vibrations are usually undesired phenomena as they may cause discomfort, disturbance, damage, and sometimes destruction of machines and structures. It must be reduced or controlled or eliminated. For this reason, the vibrations attenuation became a major issue for scientists and researchers in order to develop effective solutions for these problems. Many technologies have been developed. Among these technologies, the nonlinear vibration absorbers have significant performance in the vibration attenuation over a wide frequency band. In this context, this thesis focuses on the analysis and optimization of nonlinear vibration absorbers. The objective of the thesis is to analyze the nonlinear dynamic behavior of systems with nonlinear vibration absorbers. For this, a dynamic model of a two degrees of freedom system is developed. The Asymptotic Numerical Method (ANM) is used to solve the nonlinear equations of motion. The performance of this method is shown via a comparison with the Newton-Raphson method. The nonlinear modal analysis system with cubic nonlinearity is made by an explicit formulation of the nonlinear Frequency Response Functions (FRFs) and Nonlinear Normal Modes (MNNs). An experimental study is performed to validate the numerical results.

La maîtrise du bruit et des vibrations est un objectif fréquemment rencontré dans le domaine industriel. Qu’il s’agisse de questions de confort ou de sécurité, les domaines d’applications sont nombreux et variés : transport, BTP, ingénierie civile et militaire… Dans cette thèse, un problème de vibroacoustique interne avec couplage fluide-structure est étudié. Il s’agit d’une cavité remplie de fluide dont les parois sont constituées d’une structure sandwich viscoélastique. Les difficultés numériques associées à ce modèle portent sur la non linéarité du matériau et sur les propriétés des opérateurs matriciels manipulés (conditionnement, non symétrie). Le calcul des vibrations du système dissipatif couplé nécessite une valeur initiale, choisie comme la solution du problème conservatif. Cette solution n’étant pas aisée à déterminer, deux solveurs aux valeurs propres basés sur la Méthode Asymptotique Numérique (MAN) sont proposés pour résoudre le problème des vibrations libres du système conservatif. Associant des techniques de perturbation d'ordre élevé et de continuation, la MAN permet de transformer le problème non linéaire de départ en une suite de problèmes linéaires, plus simples à résoudre. Les solutions obtenues sont ensuite utilisées comme point initial pour déterminer la réponse libre du système dissipatif. Un solveur de Newton d’ordre élevé, basé sur les techniques d’homotopie et de perturbation est développé pour résoudre ce problème. Enfin, le régime forcé est étudié. Pour toutes les configurations envisagées, les résultats obtenus mettent en évidence des performances numériques améliorées par rapport aux méthodes classiquement utilisées (Arpack, Newton…). / Noises and vibrations control is a common objective in the industrial field. Whether it is a question of comfort or safety, the fields of application are numerous and varied: transport, building, civil and military engineering… In this thesis, a vibroacoustics interior problem with fluid-structure coupling is studied. A cavity filled of fluid whose walls are made of a sandwich viscoelastic structure is considered. The numerical difficulties associated with this model relate to the non-linearity of the viscoelastic material and the properties of the matrix operators used (conditioning, non-symmetry). The calculation of the vibrations of the coupled dissipative system requires an initial value, chosen as the solution to the conservative problem. Since this solution is difficult to determine, two eigenvalue algorithms based on the Asymptotic Numerical Method (ANM) are proposed to solve the problem of free vibrations of the conservative system. Combining high order perturbation and continuation techniques, ANM transforms the initial non-linear problem into a set of linear problems that are easier to solve. The solutions obtained are then used as the initial point to determine the free vibrations of the dissipative problem. A high order Newton solver, based on homotopy and perturbation techniques, is developed to solve this problem. Finally, the forced harmonic response of the damped system is computed. For all the configurations tested, the results obtained show improved numerical performance compared to the methods conventionally used (Arpack solver, Newton algorithm…).

L'objectif ultime de ce travail est de fournir des méthodologies robustes et efficaces sur le plan des calculs pour la modélisation et la résolution des problèmes de contact adhésifs basés sur le potentiel de Lennard-Jones (LJ). Pour pallier les pièges théoriques et numériques du modèle LJ liés à ses caractéristiques nondéfinies et non-bornées, une méthode d'adaptativité en modèle est proposée pour résoudre le problème purement-LJ comme limite d'une séquence de problèmes multiniveaux construits de manière adaptative. Chaque membre de la séquence consiste en une partition modèle entre le modèle microscopique LJ et le modèle macroscopique de Signorini. La convergence de la méthode d'adaptativité est prouvée mathématiquement sous certaines hypothèses physiques et réalistes. D'un autre côté, la méthode asymptotique numérique (MAN) est adaptée et utilisée pour suivre avec précision les instabilités des problèmes de contact à grande échelle et souples. Les deux méthodes sont incorporées dans le cadre multiéchelle Arlequin pour obtenir une résolution précise, tout en réduisant les coûts de calcul. Dans la méthode d'adaptativité en modèle, pour capturer avec précision la localisation des zones d'intérêt (ZDI), une stratégie en deux résolutions est suggérée : une résolution macroscopique est utilisée comme une première estimation de la localisation de la ZDI. La méthode Arlequin est alors utilisée pour obtenir une résolution microscopique en superposant des modèles locaux aux modèles globaux. En outre, dans la stratégie MAN, la méthode Arlequin est utilisée pour supprimer les oscillations numériques, améliorer la précision et réduire le coût de calcul. / The ultimate goal of this work is to provide computationally efficient and robust methodologies for the modelling and solution of a class of Lennard-Jones (LJ) potential-based adhesive contact problems. To alleviate theoretical and numerical pitfalls of the LJ model related to its non-defined and nonbounded characteristics, a model-adaptivity method is proposed to solve the pure-LJ problem as the limit of a sequence of adaptively constructed multilevel problems. Each member of the sequence consists of a model partition between the microscopic LJ model and the macroscopic Signorini model. The convergence of the model-adaptivity method is proved mathematically under some physical and realistic assumptions. On the other hand, the asymptotic numerical method (ANM) is adapted to track accurately instabilities for soft contact problems. Both methods are incorporated in the Arlequin multiscale framework to achieve an accurate resolution at a reasonable computational cost. In the model-adaptivity method, to capture accurately the localization of the zones of interest (ZOI), a two-step strategy is suggested: a macroscopic resolution is used as the first guess of the ZOI localization, then the Arlequin method is used there to achieve a fine scale resolution. In the ANM strategy, the Arlequin method is also used to suppress numerical oscillations and improve accuracy.

L'étude par continuation des solutions périodiques et quasi-périodiques est appliquée à plusieurs modèles issus du violon. La continuation pour un modèle à un degré de liberté avec friction régularisée permet de montrer la préservation, par rapport à la friction de Coulomb, des bifurcations de cycle limite (une vitesse maximale et une force minimale permettant le mouvement de Helmholtz) et de propriétés globales de la branche de solution (croissance de l'amplitude avec la vitesse, décroissance de la fréquence avec la force normale). L'équilibrage harmonique est évalué sur la friction régularisée et a des propriétés de convergence intéressantes (erreur faible, monotone, à décroissance rapide). La continuation sur un modèle à deux modes donne accès aux solutions de registres supérieurs, dont la stabilité coïncide avec l'expérience. La valeur retenue pour l'inharmonicité peut modifier fortement le diagramme de bifurcation. Une nouvelle méthode de continuation des solutions quasi-périodiques est proposée. Elle associe l'EH étendu à deux pulsations avec la Méthode Asymptotique Numérique. Une attention particulière est portée à la rapidité des calculs, face à la croissance rapide de la taille des systèmes à inverser. Un modèle de friction prenant en compte la température au point de contact est reformulé à l'aide d'une dérivée fractionnaire. Nous proposons une méthode de continuation de solutions périodiques de systèmes contenant des dérivées ou intégrales fractionnaires. Nous établissons une condition suffisante pour que les cycles asymptotiques du cadre causal (Caputo) soient solutions du cadre que nous avons choisi. / The continuation of periodic and quasi-periodic solutions is performed on several models derived from the violin. The continuation for a one degree-of-freedom model with a regularized friction shows, compared with Coulomb friction, the persistence of limit cycle bifurcations (a maximum bow speed and a minimum normal force allowing Helmholtz motion) and of global properties of the solution branch (increase of amplitude with respect to the bow speed, decrease of frequency with respect to the normal force). The Harmonic Balance Method is assessed on this regularized friction system and shows interesting convergence properties (the error is low, monotone and rapidly decreasing). For two modes the continuation shows higher register solutions with a plausible stability. A stronger inharmonicity can greatly modify the bifurcation diagram. A new method is proposed for the continuation of quasi-periodic solutions. It couples a two-pulsations HBM with the Asymptotic Numerical Method. We have taken great care to deal efficiently with large systems of unknowns. A model of friction that takes into account temperature of the contact zone is reformulated with a fractional derivative. We then propose a method of continuation of periodic solutions for differential systems that contain fractional operators. Their definition is usually restricted to causal solutions, which prevents the existence of periodic solutions. Having chosen a specific definition of fractional operators to avoid this issue we establish a sufficient condition on asymptotically attractive cycles in the causal framework to be solutions of our framework.

L'objectif de cette thèse est de développer des techniques de modélisation et de simulation multi-échelle avancées et efficaces pour étudier les phénomènes d'instabilité dans trois structures d'ingénierie courantes: membrane, film/substrat et structures sandwich, en combinant la technique des coefficients de Fourier lentement variables (TSVFC) et la méthode numérique asymptotique (ANM). À cette fin, basée sur les équations de la plaque de Von Karman, la TSVFC été utilisée pour développer un modèle de Fourier à bidimensionnel (2D) qui a également été implémenté dans ABAQUS via sa sous-routine UEL. Ensuite, un 2D modèle de Fourier est construit pour le film/substrat. En outre, en utilisant leurs caractéristiques de déformation, un 1D modèle de Fourier est développé en utilisant à la fois le TSVFC et le CUF. Par la suite, sur la base d'une cinématique Zig-Zag d'ordre supérieur, un 2D modèle de Fourier est déduit pour une plaque sandwich. Les équations directrices pour les modèles ci-dessus sont discrétisées par la méthode des éléments finis, et les systèmes non linéaires résultants sont résolus par le solveur non linéaire efficace et robuste ANM. Ces modèles sont ensuite adoptés pour étudier les instabilités dans ces structures. Les résultats montrent que les modèles établis peuvent simuler avec précision et efficacité divers phénomènes d'instabilité. En outre, on constate que l'instabilité membranaire est sensible aux conditions aux limites et qu'il existe un paramètre sans dimension presque constant près du point de bifurcation pour différents cas de charge et paramètres géométriques, ce qui peut être utile pour prédire rapidement l'apparition des rides / The main aim of this thesis is to develop advanced and efficient multi-scale modeling and simulation techniques to study instability phenomena in three common engineering structures, i.e., membrane, film/substrate and sandwich structures, by combining the Technique of Slowly Variable Fourier Coefficients (TSVFC) and the Asymptotic Numerical Method (ANM). Towards this end, based on the Von Karman plate equations, the TSVFC has been firstly used to develop a two-dimensional (2D) Fourier double-scale model for membrane, which has also been implemented into ABAQUS via its subroutine UEL. Then a 2D Fourier model is constructed for film/substrate. Further, making use of deformation features of the film/substrate, a 1D Fourier model is developed by using both the TSVFC and the Carrera’s Unified Formulation (CUF). Subsequently, based on high-order kinematics belonging to Zig-Zag theory, a 2D Fourier model is deduced for sandwich plate. The governing equations for the above models are discretized by the Finite Element Method, and the resulting nonlinear systems are solved by the efficient and robust nonlinear solver ANM. These models are then adopted to study instabilities in these structures. Results show that the established models could accurately and efficiently simulate various instability phenomena. Besides, it’s found that the membrane instability is very sensitive to boundary conditions, and there exists a dimensionless parameter that is almost constant near bifurcation point for various loading cases and geometric parameters, which may be helpful for fast predicting the occurrence of wrinkles

Le procédé de soudage par friction et malaxage est un procédé récent qui a été développé au sein de l'institut de soudure britannique "The Welding Institute" au début des années 90. Ce procédé, utilisé généralement en aéronautique, est sans apport de matière et permet de souder principalement des alliages d'aluminium difficilement soudables par les procédés classiques de soudage. Il consiste à malaxer le matériau de base à l'aide d'un outil constitué d'un pion et d'un épaulement frottant sur les faces supérieures des tôles à souder. La modélisation de ce procédé est très complexe puisque ce dernier implique des couplages entre des phénomènes mécaniques, thermiques et métallurgiques. Le malaxage dans le procédé de soudage FSW est difficile à simuler à l'aide de la méthode des éléments finis (en lagrangien) puisque la zone proche de l'outil de soudage est le siège de grandes déformations. Donc le remaillage est nécessaire. Cependant, le remaillage est cher et très difficile pour les problèmes tridimensionnels. Par ailleurs, après un remaillage, il est nécessaire d'interpoler les champs (vitesses, contraintes,...) correspondant à la solution courante, ce qui peut introduire des erreurs supplémentaires dans le calcul (on parle de diffusion numérique). Nous proposons dans ce travail des modèles basés sur la méthode sans maillage dite "Smoothed Particle Hydrodynamics SPH" et la méthode des moindres carrés mobiles (Moving Least Square MLS) pour la simulation de ce procédé. Ces modèles sont formulés dans le cadre lagrangien et utilisent la forme forte des équations aux dérivées partielles. Le premier modèle basé sur SPH considère la zone de soudure comme un fluide non newtonien faiblement compressible et dont la viscosité dépend de la température. Ce modèle est proposé pour la simulation numérique du comportement thermomécanique d'un matériau soudé par le procédé FSW. Dans le deuxième modèle, un algorithme itératif implicite de premier ordre a été proposé, pour simuler le malaxage de la matière dans le cas d'un matériau viscoplastique, en utilisant la méthode MLS et la technique de collocation. Le troisième modèle est un algorithme implicite d'ordre élevé basée sur le couplage de la méthode MLS et la Méthode Asymptotique Numérique MAN. Cet algorithme permet de réduire le temps de calcul par rapport à l'algorithme itératif implicite de premier ordre. La validation de ces trois modèles proposés a été faite par le code industriel Fluent / Friction stir welding is a recent process that has been developed by the British Welding Institute TWI "The Welding Institute" since 1990s. This process, generally used in aerospace, does not need additional material and allows mainly joining plates of aluminum alloys which are difficult to weld by the classical welding processes. It consists in mixing the base material using a tool comprising a pin and a shoulder which heats the plates to be welded by friction. The modeling of this process is very complex since it involves the coupling between mechanical, thermal and metallurgical phenomena. The mixing in welding process FSW is difficult to simulate using finite element method in lagrangian framework since the area near the welding tool is submitted to large deformations. So remeshing procedure is often required. However, remeshing can be very expensive and difficult to perform for three-dimensional problems. Moreover, after remeshing step, it is necessary to interpolate the fields (velocities, constraints ...) corresponding to the current solution, which may lead to additional errors in the calculation (called numerical diffusion). We propose in this work models based on meshless methods called "Smoothed Particle Hydrodynamics SPH" and Moving Least Square method for the simulation of this welding process. These models are formulated in lagrangian framework and use the strong form of partial differential equations. The first model based on SPH considers the welding zone as a weakly compressible non-newtonian fluid and whose viscosity depends on the temperature. This model is proposed for the numerical simulation of thermo-mechanical behavior of a welded material by FSW process. The second model is a first order implicit iterative algorithm proposed to simulate material mixing in the case a visco-plastic behavior using the MLS method and the collocation technique. The third model is a high order implicit algorithm based on the coupling of MLS method and Asymptotic Numerical Method (ANM). This algorithm allows reducing the computation time by comparison with the first order implicit iterative algorithm. The validation of these three proposed models was done by the industrial code Fluent

Avec l’essor de l’électronique embarquée, les câbles électriques constituentune part importante des pièces automobiles tandis que l’espace à bord n’a cessé de diminuer. Leur flexibilité requiert la prédiction de leur déformation durant leur montage afin d’éviter le contact avec d’autres pièces du véhicule et leur endommagement. Les outils actuels ne permettent pas une prédiction assez réaliste et précise de leur comportement, nécessaire dans un volume de travail très restreint. Les étapes de montage sont donc validées via la réalisation de maquettes réelles coûteuses. Cette thèsea pour but d’améliorer la simulation numérique de ces pièces souples. Nous proposonsici un code de simulation 3D basé sur un modèle de poutre géométriquement exact résolu par la méthode des éléments finis. Son originalité tient dans le couplage des quaternions pour modéliser les rotations 3D et de la méthode asymptotique numérique pour la continuation du système non linéaire qui lui confère une grande robustesse. Un banc d’essai permettant l’identification des paramètres homogénéisés nécessaires au modèle numérique et sa validation par comparaison de la géométrie finale et du chemin d’équilibre est présenté. Combinés à des développements analytiques sur les modèles de poutres avec cisaillement, les essais mènent à une évaluation critique du modèle deTimoshenko 3D pour la représentation des torons de câbles. / With on-board electronics expansion, electrical cables are an essential partof automotive pieces and the space on board has plummeted. Their flexibility requires to predict their deformation during vehicle assembly in order to avoid the contact with other pieces and damaging. Current numerical tools do not allow a realistic and accurate prediction, which is necessary in the obstructed car space. Assembly steps thus are validated on costly physical mock-ups. This thesis aims at improving numerical simulation of these flexible pieces. We herein propose a 3D algorithm based on a geometrically exact beam model solved by the finite element method. This work’s originality stands in coupling quaternions as rotational parameters and the asymptotic numerical method as nonlinear solver which results in a very robust algorithm. A test bench designed to identify the homogenized beam parameters of the numerical model and to validate it by offering a comparison on the final geometry and the equilibrium path is presented. Analytical developments on shear beams and the results of these experimental tests lead to a critical evaluation of the 3D Timoshenko model for representing stranded cables.

Le plissement est généralement observé dans les structures minces ayant un comportement de type membrane. Ces structures minces ne supportent pas d'effort de flexion et sont donc sollicitées en traction. Dans cette thèse, nous avons développé une technique de réduction de modèle pour la modélisation du plissement des membranes minces. Cette technique, basée sur les séries de Fourier à double échelle, permet de déduire d'un modèle complet de membrane, un modèle réduit capable de prendre en compte les instabilités globales et locales. Les valeurs critiques de charge et longueur d'onde sont déterminées analytiquement puis numériquement. Des exemples numériques nous ont permis de valider le modèle numérique par rapport au modèle analytique. Les modèles numériques étudiés prennent en compte le modèle complet et le modèle réduit de la membrane. Le modèle complet est simulé dans Abaqus et résolu numériquement à l'aide de la méthode de la longueur d'arc et le modèle réduit est implémenté dans Matlab et résolu numériquement à l'aide de la méthode asymptotique numérique. Nous avons étudié le comportement de la membrane sous sollicitation mécanique, thermique et thermo-mécanique. Les résultats obtenus montrent que le modèle réduit est capable de se substituer au modèle complet dans la détermination des contraintes critiques et longueurs d'onde correspondantes. Le gain en temps de calcul obtenu est important, ceci grâce à la très faible densité de maillage requis par le modèle réduit. Le modèle réduit est très sensible aux conditions aux bords de la membrane et requiert d'avoir une longueur d'onde des plis quasiment constante dans la largeur de la membrane / Wrinkling is an instability phenomenon generally observed in thin structures with membrane's behavior. Those thin structures have no rigidity to flexion and are therefore used in traction. In this thesis, we developed a reduction model's technique for the modeling of wrinkling phenomenon in thin membranes. This technique, based on the double scale Fourier series, allow us to deduce from a full membrane model, a reduced membrane model that is able to take into account the global and local instability of the structure. The critical load and critical wavelength are determined analytically on one side, then numerically on the other side. Numerical exemples are conducted to validate the numerical model towards the analytical one. Numerical models studied take into account both full and reduce membrane models. The full model is simulated in Abaqus and solved numerically using the arc length method and the reduced model is implemented in Matlab and solved numerically using the asymptotic numerical method. We studied the membrane behavior under mechanical, thermal and thermo-mechanical loading. The results obtained show that the full membrane model can be replaced by the reduced one in determining critical loads and corresponding wavelengths. The gain in computation time obtained is important, due to the coarse mesh required by the reduced model. The reduced model is very sensitive to membrane's boundaries conditions and requires to have a quasi constant wavelength along the membrane width

Le travail de cette thèse a été réalisé dans le cadre d'une collaboration entre EDF R&D et le LMA de Marseille (CNRS). Le but était de développer des outils théoriques et numériques pour le calcul de modes non-linéaires de structures industrielles possédant des non-linéarités localisées à jeu. La méthode de calcul utilisée est une combinaison de la méthode d'équilibrage harmonique (EH) et de la méthode asymptotique numérique (MAN), appelée EHMAN. Elle est réputée pour sa robustesse sur les problèmes réguliers. L'enjeu de ce travail de thèse est de l'appliquer sur des problèmes non-réguliers régularisés de type butée à jeu pour lequel un grand nombre d'harmonique est nécessaire. Des améliorations ont été apportées à la méthode de base pour rendre effectif le traitement de modèles à "grand" nombre de degrés de liberté (DDL). Les développements réalisés pendant la thèse ont été capitalisés par la création de nouveaux opérateurs dans Code_Aster.Une étude approfondie d'un système à 2 degrés de liberté a permis de faire émerger quelques caractéristiques des systèmes non-linéaires à jeu. Celles-ci ont servi entre autre à établir une méthodologie pour l'étude de systèmes à grand nombre de DDL. Pour finir, la potentialité des modes non-linéaires comme outil de diagnostic vibratoire est démontrée avec l'étude d'un tube cintré de générateur de vapeur. Le calcul des modes non-linéaires a monté l'existence d'une interaction entre un mode hors-plan (basse fréquence) et un mode plan (haute fréquence) expliquant des régimes vibratoires non-standards. Ce résultat, impossible à obtenir avec les outils de l'analyse modale linéaire, est confirmé expérimentalement. / This work is a collaboration between EDF R&D and the Laboratory of Mechanics and Acoustics. The objective is to develop theoretical and numerical tools to compute nonlinear normal modes (NNMs) of structures with localized nonlinearities.We use an approach combining the harmonic balance and the asymptotic numerical methods, known for its robustness principally for smooth systems. Regularization techniques are used to apply this approach for the study of nonsmooth problems. Moreover, several aspects of the method are improved to allow the computation of NNMs for systems with a high number of degrees of freedom (DOF). Finally, the method is implemented in Code_Aster, an open-source finite element solver developed by EDF R&D.The nonlinear normal modes of a two degrees-of-freedom system are studied and some original characteristics are observed. These observations are then used to develop a methodology for the study of systems with a high number of DOFs. The developed method is finally used to compute the NNMs for a model U-tube of a nuclear plant steam generator. The analysis of the NNMs reveals the presence of an interaction between an out-of-plane (low frequency) and an in-plane (high frequency) modes, a result also confirmed by the experiment. This modal interaction is not possible using linear modal analysis and confirms the interest of NNMs as a diagnostic tool in structural dynamics.