La sécurité routière est aujourd’hui un des objectifs principaux des constructeurs automobiles. Pour garantir un haut niveau de sécurité aux conducteurs et passagers mais aussi aux autres usages de la route, comme les piétons et les deux roues, les véhicules sont testés sous tous les angles avant leur commercialisation. Qu’ils soient virtuels ou réels, les essais jalonnent la conception et le développement d’une voiture. Aujourd’hui, la simulation numérique joue un rôle primordial dans le développement des nouveaux dispositifs et équipements. Pour que les modèles numériques offrent toujours plus de prédictivité, le domaine de la simulation numérique est en constante évolution. L’objectif de cette thèse est de contribuer à l’amélioration de la simulation des crash de véhicules. Le sujet de l’étude présentée dans ce document est composé de deux problématiques. La première consiste en la prise en compte des non-linéarités matérielles de type rupture dans le processus de dimensionnement au choc. Dans le domaine du crash, deux types de pièces sont distingués. Celles qui doivent casser en choc comme les suspensions moteur et celles dont rupture est interdite comme le carter moteur ou bien le carter d’embrayage. La rupture est un phénomène délicat à modéliser et également coûteux en temps de calcul. Par conséquent, les premiers modèles pour la simulation de la rupture ont été seulement fonctionnels (désactivation des ressorts de rupture). Avec l’arrivée de méthodes numériques telles que la méthode de suppression d’éléments, de nombreux modèles de rupture ont été développés. Comme ces modèles sont relativement récents, les bases de données matériaux sont souvent incomplètes ce qui rend difficile leur utilisation. Un des objectifs de la thèse est donc l’identification des paramètres matériaux pour le modèle de Mohr-Coulomb et le modèle de Johnson-Cook à l’aide d’essais sur éprouvettes. Une campagne d’essais sur un sous-système automobile a été également effectuée. La corrélation numérique des essais a permis de valider la modélisation« volumique » de la rupture et d’identifier des voies de progrès de la méthodologie actuelle. Le deuxième sujet présenté dans ce document est la modélisation du comportement élastique linéaire des composants qui ne doivent pas casser en choc et la prise en compte des non-linéarités géométriques de type grandes rotations et de contact. Une modélisation par des super-éléments élastiques linéaires adaptés pour les calculs explicites avec grandes rotations a été proposée pour les pièces massives de fonderie. La modélisation par super-élément est réalisée dans Radioss via un fichier.fxb. Une nouvelle méthode de réduction dérivée à partir de la méthode de MacNeal a été proposée. La base de transformation de celle-ci est constituée par des modes libres complétés par des modes de souplesse résiduelle (MLSR). Le modèle réduit à l’aide de la méthode MLSR est généré avec des routines Matlab développées dans le cadre de la thèse. Afin d’intégrer la modélisation par super-élément dans le processus de dimensionnement au choc PSA, une filière de calcul Abaqus–Matlab–Radioss a été proposée. Un processus de génération du fichier.fxb à l’aide des outils CFAO utilisés chez PSA a été établi. Afin d’optimiser la performance des calculs avec super-élément, trois nouvelles modélisations pour la gestion de contact avec super-élément ont été proposées. Une application de la modélisation par super-élément a été faite pour un modèle de test et également pour un modèle de choc de véhicule. / Nowadays, road safety is one of the main aims of car makers. To guarantee a high level of security to car drivers and passengers and also to other road users, such as for example pedestrians and two wheelers, the cars are submitted to a lot of tests before their commercialisation. Whether they are virtual or real, the experiments set out the conception and the development of a car. Nowadays, the numerical simulation plays a prime role in the development of new systems and equipments. In order to offer more and more predictivity by numerical models, the domain of the numerical simulation is inconstant evolution. The aim of this thesis is to contribute to the improvement of a car crash simulation. The subject of the study presented in this document is composed by two themes. The first one consists in inclusion of the material nonlinearities as fracture in the car crash dimensioning process. In car crash, two types of components are distinguished. Those which have to brake during the shock such as engine suspensions and those which have not to brake such as crankcase or clutch housing. The fracture is a phenomenon difficult to model, moreover, its computational time is important. Thus, the first fracture models were only functional (deactivation of fracture springs). With the arrival of numerical methods such as element deletion method, a lot of fracture models were developed. Because these models are relatively recent, the material databases are often incomplete which makes their use more difficult. Thus, one of the aims of this study is the material parameters identification of Mohr-Coulomb and Johnson-Cook fracture model with the aid of experiments on specimens. Another experiment campaign on an automotive subsystem was also realized. The numerical correlation using the material fracture parameters identified in the first step enabled to validate the "volumic"fracture modelisation and to identify the ways of progress of the current methodology. The second theme presented in this document is the modeling of linear elastic behavior of components which have not to brake in car crash and inclusion of geometrical nonlinearities as large rotations and contact. A modeling by linear elastic super elements adapted for explicit simulations with large rotations was proposed for massive cast iron components. The modeling by super elements was realized in Radioss by means of the file.fxb. A new reduction method derived from that of MacNeal was proposed. Its transformation basis is constituted by a free eigen-modes completed by residual flexibility modes. The model is reduced by means of Matlab routines which were developed in the framework of the thesis. In order to integrate the modeling by superelements in the dimensioning car crash process used at PSA, a procedure Abaqus–Matlab–Radioss was proposed. A process for the generation of the file.fxb by means of the CAD/CAMtools used at PSA was established. In order to optimize the simulation performance with superelements, three new models for contact management with a superelement were proposed. The modeling by superelement was applied in a test model and also in a frontal car crash simulation.

Ce travail de thèse présente la démarche utilisée pour construire un modèle éléments finis en statique ou en dynamique d'un pneumatique d'avion gros porteur prenant en compte la géométrie, la structure matérielle complexe, les différents matériaux et leurs propriétés ainsi que les interactions entre le pneumatique et le sol (contact, frottement et couplage thermomécanique). Des essais ont été effectués afin d'identifier les paramètres géométriques et matériaux.Ces simulations utilisant une approche Lagrangienne et une approche mixte Eulérienne/Lagrangienne ont été réalisées sur le modèle proposé. La seconde approche qui réduit considérablement le temps de calcul a été validée pour simuler le roulement en dérapage du pneumatique. Le modèle permet ainsi d'estimer le torseur des efforts dans le contact pneumatique /sol. Nous montrons l'influence des paramètres de chargement (charge verticale, pression de gonflage et vitesse de roulage) et de l'angle de dérapage sur le moment d'autoalignement (MZ) et sur le potentiel d'adhérence (µY) correspondant au rapport entre l'effort latéral et l'effort vertical dû au poids de l'avion. Nous présentons aussi une étude de sensibilité aux paramètres géométriques et matériels.Des essais de Coulomb et de diffusion thermique ont permis d'identifier la loi d'évolution du coefficient de frottement en fonction de la température (béton, asphalte) et l'évolution de la température dans l'épaisseur du pneumatique. Ceci a permis de prendre en compte les effets thermiques dans le modèle et de proposer un modèle de couplage thermomécanique qui met en évidence la décroissance de µY et la chute rapide de MZ vers des valeurs négatives au-delà d'un angle de dérapage critique βmax variant avec les conditions de chargement du pneumatique. Ces variations sont observées expérimentalement. / This PhD Thesis presents the approach adopted for the setting of numerical model based on Finite Element Method for jumbo-jet tire. The model takes into account the real geometry, the complex material structure, the various materials and their properties as well as the interactions between the tire and the ground (contact, friction and thermal-mechanical coupling due to friction). Tests are performed in order to identify geometrical and material parameters.Static and dynamic simulations using a Lagragian approach and an Eulerian/ Lagrangian mixed approach were performed on this proposed model. This second approach which significantly reduces the computational cost time was validated for cornering tire simulation. The model allows thereby to estimate the forces in the tire/ ground contact patch. We show the influency of loading parameters (vertical load, inflating pressure and rolling velocity) and of the slip angle on the self aligning torque (MZ) and on the lateral friction coefficient (µY) corresponding on the ratio between lateral force and vertical load due to the aircraft weight. We also present a sensitivity study on geometrical and material parameters.Coulomb's and thermal diffusion tests were performed in order to identify the friction coefficient law as function of temperature (on concrete and asphalte surfaces) and the temperature evolution in the aircraft tire thickness. These tests allowed to take into account thermal effects in the model and to propose a thermal-mechanical coupling model which emphasized the decreasing of µY and the rapid vanishing of MZ towards zero beyond a critical slip angle βmax varying with the tire loading conditions. These variations were observed experimentally