Ce travail de thèse, mené dans le cadre d’une convention CIFRE entre Snecma, le CERFACS et le LMFA, s’inscrit dans un contexte d’amélioration des performances et d’extension de la plage de fonctionnement des compresseurs de type axial équipant les turboréacteurs. L’une des principales difficultés rencontrée dans cette démarche concerne la maîtrise des écoulements dans la zone de jeu en tête des aubes rotors et qui peuvent entraîner une perte de stabilité du système (pompage et décollement tournant).Une solution technologique prometteuse pour améliorer la stabilité est le traitement de carter qui consiste en un dispositif passif complexe de fentes implantées au carter au droit des rotors. En vue d’en améliorer sa conception, les travaux de thèse visent plus particulièrement à approfondir la compréhension des mécanismes d’action grâce à une approche numérique CFD avec le code elsA développé par l’ONERA et le CERFACS, en modélisation stationnaire et instationnaire. Ces travaux s’articulent autour de trois axes principaux. Le premier a eu pour objectif de développer un outil numérique d’aide à la compréhension des mécanismes d’action des traitements de carter et de diagnostic de leur efficacité. Le principe de l’outil, qui est une extension du modèle initialement proposé par Shabbir et Adamczyk, repose sur une évaluation des contributions des termes des équations de Navier-Stokes stationnaires et instationnaires sur un volume de contrôle pris dans l’écoulement. Dans le cas pratique, cela revient à quantifier les efforts appliqués sur le fluide. Le second axe traite de l’analyse des mécanismes d’action des traitements de carter axisymétriques dans deux compresseurs axiaux : l’un subsonique à carter cylindrique (CREATE) et l’autre transsonique à carter conique (NASA Rotor 37). Les enseignements de cette étude indiquent que ce type de géométrie est marqué par son effet d’aspiration de fluide dans la veine. Ce mécanisme est d’autant plus amplifié par un phénomène d’interaction complexe des fentes avec l’écoulement de jeu et la proximité de l’intrados de l’aube adjacente. Cette partie s’est également attardé à la réponse des rainures à un phénomène instationnaire de type sillage de roue amont. Les résultats ont montré que les fentes amortissent les fluctuations de gradient de pression adverse. Le troisième axe porte sur l’analyse des mécanismes des traitements de carter non-axisymétriques à travers l’étude numérique d’un cas test transsonique à carter cylindrique (CBUUA). Le mécanisme d’action améliorant la stabilité de la machine tient en la capacité des fentes à limiter la migration dans la direction circonférentielle du vortex de jeu. Les résultats montrent que ce type de géométrie est caractérisé par son effet de réinjection d’air qui vient ré-énergétiser l’écoulement proche carter. / This thesis work, conducted as part of a CIFRE agreement between Snecma, CERFACS and LMFA, deals with the context of improving performance and extending the operating range of axial compressors fitted turbojets. One of the main difficulties in this approach is the flow control in the rotor tip region, which can cause the loss of the system stability (surge and rotating stall). A promising technology known to bring substantial stability is the casing treatment. This passive control device consists of slots of complex geometry within the rotor casing. In order to improve its design, the thesis aimed specifically at improving the understanding of their mechanisms through a numerical approach using the CFD code elsA developed by ONERA and CERFACS, with steady and unsteady approaches. This work focused on three main axes. The first concerns the development of a numerical tool to support the understanding of casing treatment mechanisms and the diagnosis of their efficiency. The principle of the tool, which is an extension of the model originally proposed by Shabbir and Adamczyk, is based on an assessment of the contributions of the terms of the steady and unsteady Navier-Stokes equations on a control volume taken in the flow. In practice, this permits to quantify the forces applied to the fluid. The second axis deals with the analysis of the flow mechanisms induced by axisymetric casing treatments in two axial compressors : one subsonic with a cylindrical casing (CREATE) and the other transonic with a conical casing (NASA Rotor 37). The findings of this study indicate that this type of geometry is characterized by its bleeding effect. This mechanism is further amplified by a complex phenomenon of interaction between grooves, tip leakage vortex and the proximity to the pressure side of the adjacent blade. This part has also dwelt on the groove’s response to unsteady upstream stator wakes. The results showed that the slots are able to damp fluctuations of adverse pressure gradient. The third area concerns the analysis of the flow mechanisms induced by non-axisymmetric casing treatment through the numerical study of a transonic compressor with cylindrical casing (CBUUA). The mechanism leading to an enhancement of the stability results in slots ability to limit the migration in the circumferential direction of the tip leakage vortex. The results show that this type of geometry is characterized by its effect of re-injection of fluid that comes re-energize the near casing flow.

La rectification est un procédé couramment utilisé dans l’industrie pour la finition des surfaces. L’optimisation du procédé consiste à trouver un compromis entre la qualité des pièces, minimiser les temps d'usinage et augmenter l’efficacité énergétique grâce au choix judicieux des paramètres de rectification. Par ailleurs le taux de production des pièces rectifiées est souvent limité par des contraintes sur la topographie de la surface et des problèmes liés à l’apparition de brûlures de rectification ou de micro-fissures à la surface des pièces. Ces défauts d’aspect engendrent généralement, lorsqu’ils concernent une surface fonctionnelle, une réduction de la durée de vie du composant ainsi rectifié.’effet des conditions de rectification et des propriétés des matériaux sur la nature des contraintes résiduelles a été analysé par modélisation numérique. Le modèle élément finis permet la prédiction non seulement des contraintes résiduelles, mais aussi des phases en présence et des déformations associées. L'objectif de cette étude est de construire un modèle numérique fiable en se basant sur la méthode des éléments finis pour analyser les contraintes résiduelles induites par la rectification et d'explorer, par conséquent, les mécanismes en termes de conditions de rectification. La variation des contraintes résiduelles et des déformations aux points d'intégration a été analysée. Les effets du coefficient de frottement (µ), du nombre de Peclet (Pe), de la conductance de paroi (H) et du flux de chaleur (Q) sur la microstructure et l’état de contraintes résiduelles ont été analysés. Enfin, sur la base des nouveaux résultats de ce travail de recherche, une méthodologie plus complète est proposée pour la suite. / Grinding is a commonly used finishing process to produce components of desired shape, size and dimensional accuracy. The ultimate goal is to have the maximum workpiece quality, minimum machining time and high economic efficiency by making a selective adaptation of the possible process strategy and chosen parameter selection. The focus of this study arose from a limitation that challenges the grinding industry. The production rate of the ground parts is generally constrained by surface topography and subsurface damage appearing as residual tensile stress, localized burns, and phase transformation induced micro and macro-cracking. This motivates the need for a reliable numerical modelling to simulate the grinding process. The numerical model sought should be able to predict not only the required grinding residual stresses but also the deformation history. The objective of this thesis is to build up a reliable finite element model for grinding-induced residual stress analysis and thus to explore thoroughly the mechanisms in terms of grinding conditions. The variations of the residual stresses and strains at integration points have been examined, and the effects of the friction coefficient (µ), Peclet number (Pe), non dimensional heat transfer coefficient (H) and different magnitudes of input heat flux (Q) on both the microstructure and the residual stress state are analyzed. Finally, based on the new findings in this research, a more comprehensive methodology is suggested for further study.

Afin d’augmenter le rendement et la puissance spécifique des turbines à gaz, la température en entrée de la turbine haute pression a été continuellement augmentée, au point d’atteindre une valeur de l’ordre de la température de fusion des matériaux. La capacité à prédire l’écoulement (température en particulier) dans la turbine est donc un élément essentiel pour la conception des turbines à gaz, particulièrement pour celle qui ne sont pas refroidies. Toutefois, cette tâche est rendue difficile par l’extrême complexité de l’environnement dans lequel la turbine évolue (interaction avec la chambre de combustion, effets technologiques, couplage thermique fluide / solide, etc.).L’objectif de ce travail de thèse est de proposer, d’une part une stratégie de simulation numérique afin de prédire au mieux les champs aérothermiques dans une turbine haute-pression, d’autre part une méthode d’analyse permettant de quantifier l’impact des différents éléments de l’environnement sur les performances. Pour répondre à ces objectifs, des simulations instationnaires d’un étage de turbine ont été réalisées avec le code elsA, prenant en compte les effets technologiques (baignoire, refroidissement externe grâce aux évents du distributeur et aux cavités moyeu et carter du rotor qui sont alimentés par le système d’air secondaire, congés de raccordement) et les transferts thermiques conjugués. Une analyse de la production locale d’entropie a aussi été menée afin de comparer les performances aérodynamiques pour plusieurs niveaux d’approximations dans la définition de la turbine (prise en compte ou non de tel ou tel effet technologique) et de localiser l’origine de ces différences. L’analyse des résultats a montré que la hauteur de jeu et le système de refroidissement externe ont l’impact le plus significatif sur les performances aérodynamiques de la turbine. La température de paroi de la pale de rotor est de son côté fortement influencée par l’écoulement issu de la chambre de combustion, le refroidissement externe et le couplage thermique fluide / solide. Ce travail est un premier pas vers la réalisation de simulations totalement intégrées qui doivent permettre d’améliorer la précision des conceptions. / In order to increase the thermodynamic efficiency of gas turbine engines, the high-pressure turbine inlet temperature has been continually increased up to reach levels of the order of magnitude of the vanes and blades melting temperatures. The ability of predicting the flow through the turbine(especially the temperature) is a key point for the design of gas turbines, especially for the uncooled ones. However, this is challenging because of the complex environment that interacts with the turbine(hot-streak migration, technological details, fluid/solid thermal coupling . . .).The aim of this work is to develop a strategy based on CFD in order to predict aerothermal fields in a high-pressure turbine as well as an analysis allowing to quantify the impact of the environment on the turbine performances. To achieve these goals, the elsA code has been used to perform unsteady simulation of a turbine stage, taking into account technological details (squealer tip, external cooling thanks to vane trailing edge cooling and rotor hub and shroud cooling cavities, fillets) and conjugate heat transfers. An analysis of the local entropy production rate was also used to compare the performances resulting from different modeling of the turbine (taking into account or not one particular technological detail) and to localize the origin of these discrepancies. The results show that the tip height and the external cooling have the greatest impact on the turbine performances. The rotor blade temperature is mainly affected by the flow coming from the combustion chamber, the external cooling and the fluid / solid thermal coupling. This work is the first step towards the realization of integrated simulations allowing to improve the accuracy of design.

Les procédés additifs, auxquels appartient la fusion laser de poudres métalliques, ont la capacité de créer des structures à géométries complexes, avec la possibilité d'intégrer des formes creuses, par exemple des canaux de refroidissement assurant un contrôle thermique optimum. Ce procédé permet de fabriquer des pièces réelles à partir de poudres métalliques, par fusion du matériau, couche par couche, en accord avec le modèle CAO. Au cours du procédé, de nombreux cycles thermiques et d'importants gradients thermiques se produisent dans la pièce au cours de sa fabrication. Ces gradients de température induisent des déformations plastiques hétérogènes et de ce fait des contraintes résiduelles. Ces contraintes peuvent nuire à la qualité de la pièce obtenue, par exemple sa résistance mécanique. Ces travaux ont pour objectifs de proposer un modèle numérique, s’appuyant sur la méthode des éléments finis afin d'étudier l'apparition des contraintes résiduelles lors du procédé de fusion laser de poudres métalliques. Le logiciel multiphysique ABAQUS® a été utilisé pour effectuer les analyses thermiques et mécaniques. La technique « d'ajout et de suppression des éléments » a été utilisée afin de simuler la fusion et la solidification de la matière au cours du procédé. Les propriétés mécaniques dépendantes de la température de l'acier maraging, utilisé dans notre cas, ont été obtenues à l’aide d’essais expérimentaux de caractérisations et intégrées dans le modèle. Les calculs sont réalisés de manière découplée, dans un premier temps le calcul thermique est effectué, puis les résultats sont utilisés pour réaliser le calcul mécanique et finalement prédire les champs de contraintes. Dans le cadre de ce travail, une méthode originale s'appuyant sur la technique de mesure des contraintes résiduelles par enlèvement de couches successives a été mise au point pour mesurer ces contraintes en direct au cours du procédé. Les résultats renseignent sur le niveau et la distribution des contraintes dans la pièce créée et le support. Deux paramètres ont été testés afin d'étudier leur influence sur le niveau des contraintes résiduelles : le temps d’étalement de la poudre entre deux couches successives et la hauteur des couches. Le modèle numérique paramétrable permet d'analyser les effets de paramètres liés au procédé sur la répartition des contraintes résiduelles dans les pièces fabriquées. Les résultats montrent que la variation de l'épaisseur du support n'affecte pas la répartition des contraintes dans la pièce créée. Le préchauffage du support à une température de 800°C réduit les contraintes résiduelles L'étude de quelques trajectoires laser montre leurs influences sur la répartition des déformations plastiques cumulées ainsi que la hauteur des couches de poudres ou de la forme du support (embase, colonnes). / The Selective Laser Melting process, belonging to Additive processes , have the ability to create structures with complex geometries , with the possibility of including cavities, such as cooling channels providing optimum temperature control. This process enables the manufacture of three-dimensional parts from metal powders by melting the material , layer by layer, in agreement with the CAD model. In the process , high temperatures and thermal gradients cycles occur in the part during the process. These temperature gradients induce heterogeneous plastic strain and residual stresses. These residual stresses may affect the quality of the part obtained, for example the fatigue life. This work aims to propose a numerical model , based on the finite element method to study the appearance of residual stresses during laser melting process of metallic powders . The ABAQUS® Multiphysics software was used to perform the thermal and mechanical analyzes. The movement of the laser beam and the resolution of the thermal problem can predict the evolution of the temperature in the part and support. The "birth and death elements" technique was used to simulate the melting and solidification of the material during the process. Dependent mechanical properties of the temperature of the maraging steel used in this case were obtained using experimental testing and characterization and were established in the model. The calculations are decoupled : initially thermal calculation is performed and the results are used to perform mechanical calculations and finally predict the residual stress fields. In this work, a novel method based on the technique of measuring residual stresses by removing layers was developed to measure these stresses directly in the process. The results provide information on the level and distribution of stresses in the created part and support. Two parameters were tested to study their influence on the level of residual stress : time to spread the powder between two successive layers and layer height. The model is used to analyze the effects of process parameters related to the distribution of residual stresses in the manufactured parts. The results show that the variation of the thickness of the support does not affect the distribution of stresses in the part created. Preheating the substrate to a temperature of 800 °C reduces the residual stresses. The study of some laser strategies shows their influence on the distribution of plastic strain thus the height of the layers of powder or in the form of support (base, columns).

L’amélioration des performances des turboréacteurs actuels est un enjeu crucial dans un contexte de contraintes économiques et environnementales fortes. Au sein du turboréacteur, le canal inter-turbines, localisé à l’interface entre la turbine Haute Pression (HP) et le premier distributeur Basse Pression (BP), est le siège d’écoulements très complexes. Ainsi, les structures aérodynamiques issues de la turbine HP (sillages, tourbillons et ondes de choc) interagissent fortement entre elles et impactent l’écoulement du distributeur BP, engendrant ainsi des pertes de rendement de l’ensemble de la configuration. Ce travail de thèse s’attache à étudier les phénomènes d’interactions aérodynamiques entre une turbine HP et le premier distributeur BP et à analyser les mécanismes à l’origine des pertes aérodynamiques dans le distributeur BP. Une campagne expérimentale antérieure, réalisée sur un banc d’essai comprenant une turbine HP couplée à un distributeur BP, avait permis de recueillir des mesures de l’écoulement dans des plans situés dans le canal inter-turbines et à l’aval du distributeur BP. En lien avec ces résultats expérimentaux, les simulations numériques menées dans cette étude avec le logiciel elsA s’attachent à restituer précisément la nature tridimensionnelle, instationnaire et turbulente de l’écoulement au sein de cette même configuration. Ces travaux se développent alors en trois étapes principales. Dans un premier temps, une étude stationnaire avec traitement plan de mélange permet de comprendre et quantifier les aspects généraux de l’écoulement. Une évaluation de l’effet de la modélisation turbulente RANS (Reynolds-Averaged Navier-Stokes) et du schéma numérique spatial sur les structures aérodynamiques présentes dans la configuration est réalisée. Dans un deuxième temps, une modélisation turbulente avancée de type ZDES (Zonal Detached-Eddy Simulation) est employée pour la résolution de l’écoulement dans le distributeur BP. Les structures aérodynamiques instationnaires issues de la roue HP amont sont modélisées par une condition limite à l’entrée du domaine de calcul. L’approche ZDES est comparée à une approche Unsteady RANS (URANS) sur la même configuration. La formation et la dissipation des sillages et des tourbillons est significativement différente entre les deux modélisations, ce qui impacte de manière importante la génération des pertes aérodynamiques. Enfin, des simulations URANS de plusieurs configurations permettent de mieux comprendre les effets d’interaction entre les différentes rangées d’aubes. Ainsi, les approches instationnaires chorochroniques prenant en compte un seul rotor et un seul stator évaluent des effets instationnaires importants dans le canal inter-turbines. Ces approches conduisent à la mise en oeuvre d’un calcul sur une configuration multipassages-chorochronique prenant en compte les deux stators et le rotor afin de modéliser complètement les interactions déterministes existantes. Afin de quantifier celles-ci avec précision, une décomposition modale du champ instationnaire est mise en place. Les niveaux d’interactions liées aux différentes roues sont alors quantifiés et l’impact sur les pertes aérodynamiques est évalué. / Improving the performance of current aeronautical turbines is an important issue in a context of severe economical and environmental constraints. In a turbofan, the inter-turbine channel which is located between the High-Pressure (HP) turbine and the first Low Pressure (LP) vane is characterized by a complex flow. Therefore aerodynamic structures coming from the HP turbine (wakes, vortices and showkwaves) strongly interact between each other and affect the LP vane flow field. This generates efficiency losses of the overall configuration. This PhD thesis aims at studying the aerodynamic phenomena between a HP turbine and the first LP vane and at analyzing the mechanisms creating aerodynamic losses. A previous experimental campaign, which was carried out on a facility including a HP turbine coupled to a LP vane, enabled to gather flow field measurements in planes located in the inter-turbine channel and downstream of the LP vane. In comparison with these experimental data, the numerical simulations done with elsA software intend to reproduce accurately the 3D, unsteady and turbulent nature of the flow within this configuration. The work can be divided into three mains steps. As a first step, steady simulations with a sliding mesh treatment enable to understand the general aspects of the flow. An assessment of the effects of RANS (Reynolds-Averaged Navier-Stokes) turbulent predictions and of spatial numerical schemes on the aerodynamic structures present in the configuration is carried out. As a second step, the advanced turbulence approach ZDES (Zonal Detached-Eddy Simulation) is considered for the LP vane flow prediction. The unsteady aerodynamic structures coming from the upstream HP rotor are set as an inlet boundary condition of the computational domain. The ZDES approach is compared to a URANS (Unsteady RANS) approach on the same computational domain. The generation and dissipation of the wakes and vortices are significantly different on the two simulations, and thus impact the creation of aerodynamic losses. Finally, URANS simulations enable to better understand the interaction effects between the different blade rows. First, the unsteady phase-lagged approaches that take into account a single rotor and stator assess the important unsteady effects in the inter-turbine channel. They finally lead to the implementation of a multipassages phase-lagged computation that takes into account the two stators and the rotor in order to model all the existing determinist interactions. In order to quantify them accurately, a modal decomposition of the unsteady flow field is set up. The interaction levels linked to the different blade rows are therefore quantified and the impact of the aerodynamic losses is evaluated.

Des simulations numériques rapides et précises du contact élastohydrodynamique (EHD) sont recherchées pour aider au développement de produits. L'objectif de cette thèse est de proposer un modèle compact pour le problème du contact EHD en appliquant des méthodes de réduction de modèle. Dans ce but l'équation de Reynolds (non-linéaire), l'équation d'élasticité (linéaire) et l'équilibre de la charge, sont résolus dans un système d'équations unique par la méthode de Newton. La réduction s'effectue par projection sur un sous espace de faible dimension, qui repose sur des solutions du système complet. De plus, une approximation du système est effectuée, dans laquelle les matrices du système réduit sont approximées. Pour le problème du contact EHD stationnaire, un algorithme de génération automatique des modèles compacts est présenté. L'algorithme fournit des modèles réduits stables et rapides sur une région de paramètres définies. La méthode de Newton réduite est également étendue aux fluides non-newtoniens. Les résultats du modèle réduit sont en très bon accord avec ceux du système complet, malgré un temps de calcul clairement plus petit. Par ailleurs, une nouvelle formulation pour le problème de contact EHD transitoire est introduite, dans laquelle la région de calcul est adaptée à la taille du contact. Ceci permet d'obtenir des modèles réduits efficaces, en particulier pour des excitations à grandes amplitudes. Alternativement, la méthode "Trajectory-Piecewise-Linear" (TPWL) est appliquée au problème du contact EHD transitoire. Cette méthode permet une accélération du calcul conséquente. / In today's product development process, fast and exact simulational models of complex physical problems gain in significance. The same holds for the elastohydrodynamic (EHD) contact problem. Thus, the objective of this work is to generate a compact model for the EHD contact problem by the application of model order reduction. Thereto, the EHD contact problem, consisting of the nonlinear Reynolds equation, the linear elasticity equation and the load balance, is solved as a monolithic system of equations using Newton's method. The reduction takes place by projection onto a low-dimensional subspace, which is based on solutions of the full system. Moreover, a so-called system approximation is executed at which the reduced system matrices are substituted by less complex surrogates. For the stationary EHD contact problem, an algorithm for the automated generation of the compact model is presented. This algorithm provides fast and numerically stable reduced systems on a given parameter range. Additionally, the reduced Newton method is extended to the consideration of Non-Newtonian fluids whereat highly accurate results are obtained requiring a very low computational time. Furthermore, a new formulation for the transient EHD contact problem is introduced, at which the computational area is adapted to the current contact size. This kind of morphing enables efficient reduced models in particular for excitations of large amplitude. Beside of the reduced Newton-method with system approximation, the method Trajectory Piecewise Linear (TPWL) is applied to the transient EHD contact problem. Here, further speed-up potential arises. Despite a distinctly lower computational time, the reduced model is in very good accordance with the full system.

Le but de cette thèse est d'étudier la modélisation détaillée des écoulements qui ont lieu lors des inondations urbaines. Dans une première partie, des écoulements en bifurcation incluant des petits obstacles génériques ou des profils de canaux avec trottoirs sont étudiés sur une maquette expérimentale, puis simulés numériquement avec le modèle bidimensionnel Rubar20. Les résultats expérimentaux et numériques montrent l'avantage d'inclure des obstacles de petite taille dans un modèle d'inondation urbaine, alors qu'il n'y a qu'un intérêt limité à utiliser une topographie détaillée des rues. Dans une deuxième partie, les interactions entre écoulements de surface et écoulements en conduites souterraines sont étudiées. Un modèle physique de système de drainage urbain permet de valider un modèle analytique prédisant les débits d'échange entre les deux couches d'écoulement. Une modélisation 1D/2D (conduite/rue) est mise en place avec les modèles Rubar3/Rubar20 et validée sur des écoulements expérimentaux observés sur le modèle physique. Dans une troisième partie, les inondations dans la ville d'Oullins (près de Lyon, France) sont étudiées. La modélisation des écoulements de surface est validée avec des données de terrain, et nous discutons l'intérêt de plusieurs représentations du milieu urbain. L'intégration du réseau d'assainissement dans un modèle 1D/2D reste affectée par plusieurs incertitudes, mais cette étape montre l'intérêt de la modélisation couplée pour décrire les interactions complexes des écoulements lors d'inondations urbaines, ainsi que les limites de l'approche développée pour les écoulements à faible profondeur / Aim of this thesis is to study the detailed modelling of flows that occur during urban floods. In a first part, bifurcation flows including small obstacles or channel profiles with sidewalks are studied on an experimental facility, and then numerically simulated with the two dimensional model Rubar20. Experimental and numerical results show the benefits of including small obstacles in an urban flood model, whereas there is only little benefit of using a detailed representation of the streets topography. In a second part, interactions between surface and underground pipe flows are studied. A physical model of an urban drainage system allows the validation of an analytical model predicting exchange discharges between both flow layers. A 1D/2D modelling (pipe/street) is set up with the models Rubar3/Rubar20 and validated on experimental flows observed on the physical model. In a third part, floods in the city of Oullins (near Lyon, France) are studied. Surface flows modelling is validated with field data, and we discuss the interest of several representations of the urban area. Integration of the sewer system in a 1D/2D model remains impacted by several uncertainties, yet this step shows the interest of the coupled modelling to describe complex flows interactions during urban floods, as well as limitations of the developed approach for shallow flows

En France, en 2015 les motocyclistes représentent 43% des blessés graves victimes d'accidents de la circulation routière. Au même moment, même si les accidents de chocs de motocyclistes contre les dispositifs de retenue routiers ne représentent en moyenne que 10% des accidents de la circulation routière, les motocyclistes présentent un risque de décès 30 à 80 fois supérieur au risque de décès des usagers de véhicules impliqués dans le même type d'accidents.Pour renforcer la sécurité des motocyclistes, les dispositifs de retenue pour automobiles (dispositifs traditionnels) sont équipés à leur partie inférieure d'écran de protection motards. Les écrans motards ont pour but d'amortir de manière progressive les chocs des motocyclistes en les empêchant dans le même temps de passer à travers la partie inférieure des dispositifs de retenue pour automobiles, de percuter un poteau et d'être sévèrement blessés. Les écrans motards doivent subir des essais de chocs normatifs et répondre à certaines exigences afin de pouvoir être installés sur le bord des routes européennes et de garantir surtout la sécurité des motards. Cependant, il n'existe que des procédures normatives pour tester ces dispositifs, la procédure LIER dans le contexte Français et la procédure normative EN 1317-8 dans le contexte Européen. Afin d'évaluer ces deux procédures normatives et de proposer de nouvelles dispositions pouvant servir de procédure normative plus avancée sur le protocole d'essais d'évaluation de performances d'écrans motards, des simulations numériques ont été réalisées. Un modèle numérique complet de dispositif de retenue routier équipé d'écran motard en acier a été développé et corrélé sur un essai de choc réel sur un dispositif de même type. Une fois corrélé,le modèle a été utilisé afin d'évaluer l'influence des conditions d'essais et des caractéristiques mécaniques matériaux de l'écran sur les performances de l'écran motard. Enfin, une nouvelle démarche d'évaluation de performances d'écran de protection motard a été mise en place / In France, in 2015,motorcyclists represent 43% of seriously injured accident victims in road traffic. At the same time, even though motorcyclists impact crashes against the road restraint systems represent on average 10% of accidents in road traffic, motorcyclists present a risk of death from30 to 80 times upper than the risk of death for users of vehicles involved in the same type of accidents. To enhance the safety of motorcyclists, Vehicle Restraints Systems (VRS) arefitted to the lower part of restraint systems for motorcyclists. Motorcyclist’s Restraint Systems (MRS) are designed to absorb shocks gradually and to preventmotorcyclists to pass through the lower part of VRS, to hit a post and be severely injured. Before being installed on the roadsides, MRS have to be crash-tested according to standards in order to evaluate their safety and ensure motorcyclists safety. However, there are only normative procedures to testMRS, the LIER procedurein the French context and the procedure EN1317-8 in the European context. To assess both normative procedures and propose new provisions which can be used as normative procedure more advanced on the protocol of MRS performances evaluation tests, numerical simulations were performed. A numerical model of MRS was developed and correlated on one real impact test on other device of same type. Once correlated, the model was used to conduct a parametric study on the test conditions and thematerialsmechanical properties of the MRS. Finally, a new demarche of MRS performances evaluation was developed

La construction des biosystèmes d’oxydation passive du méthane (BOPM) est une option économique et durable pour réduire les émissions de méthane des sites d’enfouissement de déchets et des effets subséquents du réchauffement climatique. Les BOPM sont constitués de deux couches principales: la couche d'oxydation du méthane (MOL) et la couche de distribution du gaz (GDL). L'oxydation du méthane se produit dans la MOL par les réactions biochimiques des bactéries méthanotrophes, et la GDL est construite sous la MOL pour intercepter et distribuer les émissions fugitives de biogaz à la base de la MOL. Fondamentalement, l'efficacité d'un BOPM est définie en fonction de l'efficacité d'oxydation du méthane dans la MOL. Par conséquent, il est indispensable de fournir des conditions adéquates pour les activités bactériennes des méthanotrophes. En plus des paramètres environnementaux, l'intensité et la distribution du biogaz influencent l'efficacité des BOPM, et ils peuvent rendre le matériau de la MOL - avec une grande capacité d'accueillir les activités bactériennes - inutilisables en termes d'oxydation du méthane sur place. L'effet de barrière capillaire le long de l'interface entre la GDL et la MOL peut provoquer des émissions localisées de méthane, due à la restriction ou la distribution non uniforme de l’écoulement ascendant du biogaz à la base de la MOL. L'objectif principal de cette étude est d'incorporer le comportement hydraulique non saturé des BOPM dans la conception des BOPM, afin d’assurer la facilité et la distribution adéquates de l'écoulement du biogaz à la base de la MOL. Les fonctions de perméabilité à l'air des matériaux utilisés pour construire la MOL des BOPM expérimentaux au site d’enfouissement des déchets de St Nicéphore (Québec, Canada), ainsi que celles d'autres de la littérature technique, ont été étudiés pour évaluer le comportement d'écoulement non saturé du gaz dans les matériaux et pour identifier le seuil de migration sans restriction du gaz. Ce dernier seuil a été introduit en tant que un paramètre de conception avec lequel le critère de conception recommandé ici, c’est à dire la longueur de la migration sans restriction de gaz (LMSG), a été défini. La LMSG est considérée comme la longueur le long de l'interface entre la GDL et la MOL où le biogaz peut migrer à travers la MOL sans restriction. En réalisant des simulations numériques avec SEEP/W, les effets de la pente de l'interface, des paramètres définissant la courbe de rétention d'eau, de la fonction de la conductivité hydraulique du matériau de la MOL sur la valeur de la LMSG (représentant la facilité d'écoulement du biogaz à l'interface) et de la distribution de l'humidité (et par conséquent celle du biogaz) ont été évalués. Selon les résultats des simulations, la conductivité hydraulique saturée et la distribution des tailles de pores du matériau de la MOL sont les paramètres les plus importants sur la distribution de l'humidité le long de l'interface. Ce dernier paramètre influe également sur la valeur du degré de saturation et donc la facilité du biogaz à la base de la MOL. La densité sèche du matériau de MOL est un autre paramètre qui contrôle la facilité d'écoulement ascendant du biogaz. Les limitations principales de la présente étude sont associées au nombre de matériaux de MOL testés et à l'incapacité de SEEP/W de considérer l'évapotranspiration. Toutefois, compte tenu des hypothèses raisonnables dans les simulations et en utilisant les données de la littérature, on a essayé de réduire ces limitations. En utilisant les résultats des expériences et des simulations numériques, des étapes et des considérations de conception pour la sélection du matériau de MOL et de la pente d'interface ont été proposées. En effet,le comportement hydraulique non saturé des matériaux serait intégré dans les nécessités de conception pour un BOPM efficace, de sorte que la capacité maximale possible d'oxydation du méthane du matériau de la MOL soit exploitée. / Abstract : Implementation of passive methane oxidation biosystems (PMOB) is a cost effective and sustainable solution to eliminate the methane emissions of landfills to the atmosphere and ensuing global warming effects. PMOBs consist of two main layers: methane oxidation layer (MOL) and gas distribution layer (GDL). The oxidation of methane occurs in MOL through the biochemical reactions of methanotrophic bacteria, and GDL is constructed beneath the MOL to intercept and distribute the fugitive biogas emissions at the base of MOL. Basically, the efficiency of a PMOB is defined based on the methane oxidation efficiency in MOL. Therefore, it is indispensable to provide adequate conditions for the bacterial activities of methanotrophs. In addition to the environmental parameters, the intensity and the distribution of the biogas reaching the MOL material influence the efficiency of PMOBs, and they may cause the MOL material possessing great capacity to host the bacterial activities to be unserviceable in terms of in field methane oxidation. The capillary barrier effect along the GDL MOL interface may provoke localized surface methane emissions, resulted from the restricted and/or non uniform distribution of upward flow of biogas at the base of MOL. The main focus of present study is to incorporate the unsaturated hydraulic behavior of PMOBs into the design of PMOBs, providing adequate ease and distribution of upward flow of biogas at the base of MOL. The air permeability functions of the materials used to construct the MOL of experimental PMOBs at the St Nicephore landfill (Quebec, Canada), along with other materials from the technical literature, were studied to evaluate the unsaturated gas flow behavior of the materials and to identify the threshold of unrestricted gas migration. This latter threshold was introduced as a design parameter based on which the recommended design criterion herein, i.e. the length of unrestricted gas migration (LUGM), was defined. LUGM is considered as the length along the GDL MOL interface along which biogas can migrate upwards without restriction. Performing sets of numerical simulations in SEEP/W, the effect of slope of interface and the parameters defining the water retention curve and hydraulic conductivity function of MOL material on value of LUGM (representing the ease of upward flow of biogas at the interface) and distribution of moisture (and therefore biogas) along the GDL MOL interface were assessed. The saturated hydraulic conductivity and the pore size distribution of the MOL material were the most influencing parameters in distribution of moisture along the interface. The latter parameter influences also the value of degree of saturation and therefore, the ease of biogas at the base of MOL. Dry density of MOL material is another parameter that controls the ease of upward flow of biogas. The main limitations of the present study are associated with the number of tested MOL materials and the inability of SEEP/W in considering the evapotranspiration. However, considering reasonable assumptions in simulations and using the data from the literature, it was attempted to reduce the limitations. Based on the results of experiments and numerical simulations, some design steps and considerations for selection of the MOL material and the slope of interface were suggested that incorporate the unsaturated hydraulic behavior into the design necessities for an efficient PMOB so that the maximum possible methane oxidation capacity of MOL material is exploited.

La tenue des surfaces des contacts rugueux roulants est un problème crucial l’évaluation de la durée de vie des mécanismes. Cette durée de vie est conditionnée dès les premiers cycles par le rodage puis par les mécanismes de fatigue des surfaces. Le rodage est défini par le temps nécessaire à l’accommodation géométrique des surfaces rugueuses entre elles, à l’interface du contact. La charge transmise sur une faible aire de contact par rapport à l’aire apparente, crée des pressions importantes qui induisent de fortes contraintes en couche superficielle et des déformations plastiques de la microgéométrie. Cette plastification a lieu dans les tous premiers cycles puis la surface se stabilise, c’est le rodage. La répétition cyclique des sollicitations au cours du fonctionnement conduit enfin à l’endommagement du matériau et des avaries en surface telles que des micro-écailles. Après une étude bibliographique sur le contact roulant rugueux et les dispositifs expérimentaux existants, la difficulté de ce type d’analyse est mise en évidence. Elle consiste à effectuer un suivi en continu de l’évolution de l’état de surface du contact à une échelle suffisamment fine et précise. Une micromachine bi-disque a été développée afin de réaliser ce suivi quasi "in-situ" à l’échelle des rugosités permettant d’identifier les mécanismes de rodage et de dégradation. Un protocole expérimental précis permet de mesurer les surfaces antagonistes dans les premiers cycles correspondant au rodage. Les surfaces vierges mesurées sont utilisées comme paramètre d’entrée d’une simulation numérique du contact rugueux d’une sphère sur un plan. La déformée de surface numériquement obtenue à l’état stabilisé est comparée à celle mesurée expérimentalement à la fin du rodage. La très bonne superposition de ces résultats permet de valider cette méthode et les résultats numériques tels que les contraintes résiduelles et déformations plastiques. Les surfaces à l’état stabilisé obtenues, sont exploitées à travers différents critères de fatigue multiaxiaux. Les résultats numériques sont également comparés aux observations expérimentales pour déterminer le critère le plus adapté à cette analyse et permettant d’expliquer la formation de fissures et d’avaries de surfaces. / The surface life of rolling rough contacts is an important problem in the evaluation of the life expectancy of a machine. This life span is conditioned by the first cycles of the running-in process and then by the surface fatigue. The running-in period is defined by the time necessary for the rough surfaces to accommodate. The real area of contact is small compared with the apparent area, hence the load creates important pressures which lead to important stresses in the superficial layer and to plastic deformation of the microgeometry. The plastic deformation takes place over the first cycles then the surface stabilizes, this is the end of running-in process. The repeated cyclic loading finally leads to material damage below the surface and to surface micropitting. After a bibliographical study on the rough rolling contact and the existing experimental test machines, the difficulty of analysing the roughness evolution is pointed out. It requires a precise, continuous monitoring of the contact surface evolution on a small enough scale. A two-disk micro-machine was developed to perform this almost "in situ" monitoring at the roughnesses scale, allowing one to identify the mechanisms of running-in and surface degradation. An accurate experimental protocol allows one to measure the opposing surfaces in the first cycles corresponding to the running-in period. The initial surfaces are used as entrance parameters for a numerical simulation of the rough contact of a sphere on a plane. The deformed surface numerically obtained in the stabilized state is compared with the measured one at the end of the running-in period. The very good agreement between these results allows one to validate this method and the numerical results such as the residual stresses and the plastic deformation. Different multiaxial fatigue criteria are applied to the numerical results obtained in the stabilized state. The results are compared to the experimental observations to determine the criterion that is the most suited for this analysis and allows one to explain the crack formation and surfaces damage.