Cette recherche porte sur le développement d’un modèle de structure de géométrie complexe équipée de composants à comportement non linéaire viscoélastique dans le but de simuler sa réponse à des excitations définies par des densités spectrales de puissance (DSP). L’application industrielle concerne l’isolation vibratoire d’un module de refroidissement automobile monté sur plot de suspension en élastomère. Une revue du comportement des élastomères en fonction de leurs conditions d’environnement et de sollicitations identifie les paramètres des différents modèles analysés. Des essais préliminaires ont été menés pour définir les intervalles des niveaux de sollicitations et quantifier l’échauffement des plots. La caractérisation expérimentale de la suspension porte sur des plots en élastomère munis de leurs interfaces afin d’agréger dans un seul modèle les non linéarités du comportement viscoélastique, des glissements et des frottements. Les boucles effort-déflexion axiales et radiales mesurées sont traitées avec un système expert développé spécialement pour caler les paramètres du modèle retenu. Il s’agit du modèle de Dahl généralisé qu’il a fallu étendre aux aspects viscoélastiques. Ce processus de calage automatique a été codé avec un logiciel développé sous Octave/Matlab. Des méthodes d’interpolations et extrapolations rendent opérationnel ce modèle sur toute la gamme fonctionnement de l’application définie dans l’espace fréquence-déflexion. Ce processus a été codé dans le module UserSubroutine pour Abaqus. Soumettre le système mécanique non linéaire à des excitations définies par une DSP nécessite de s’appesantir sur le traitement des vibrations aléatoires. En effet il faut, pour calculer les réponses, considérer le passage fréquence-temps et inversement pour les comparer éventuellement aux exigences des normes. De plus, la taille et la complexité du modèle EF de la structure industrielle rendent impossible une résolution temporelle sur l’ensemble de ses degrés de liberté. Il s’agit alors de faire appel à des techniques d’homogénéisation et de condensation dynamique afin de prévoir la réponse aux excitations à large bande fréquentielle dans le but d’analyser les performances de l’isolation vibratoire. / This research work regards the development of a complex structure model with non-linear viscoelastic components. The purpose of this study is to simulate the response of this structure submitted to a random vibration excitation based on a power spectral density definition (PSD). The industrial applicative case is the vibratory insulation of a automotive engine cooling module supported by elastomer mounts. A brief review of elastomers behavior depending on solicitations types enables to identify the parameters of the different investigated models. Preliminary tests have been conducted to define the range of amplitudes of excitations and evaluate the internal warming of rubbers during the full structure validation test. The experimental characterization of the suspension is based on rubbers mounts and their interfaces with the cooling module, in order to take into account in a unique model all nonlinearities due to the viscoelastic behavior, the slidings, and the friction. Measured force-deflection hysteretic cycles in axial and radial direction are post-processed with an expert system developed to obtain the parameters of the retained model: the modified Dahl’s model, generalized to viscoleastic aspect. This process has been developed with Octave/Matlab code. Interpolation and extrapolation methods enable to obtain a good model response on the global operating range. These methods have been coded in an Abaqus UserSubroutine. Imposing random vibration excitation of a non linear mechanical system based on PSD imposes to take into account signal processing aspects. To evaluate response levels versus norms requirements, it’s mandatory to consider the time-frequency transfer. In addition, the size and the complexity of the total finite element model of the industrial structure don’t allow a global resolution in the time domain for all the degrees of freedom. Homogenization and dynamic reduction techniques are used to evaluate the response of the system submitted to large frequency range excitations, and to analyse the behavior of the suspension.