On étudie l'écoulement turbulent d'un fluide sur une canopée, que l'on modélise comme un milieu poreux déformable. Ce milieu poreux est en fait composé d'un tapis de fibres susceptibles de se courber sous la charge hydrodynamique du fluide, et ainsi de créer un couplage fluide-structure à l'échelle d'une hauteur de fibre (honami). L'objectif de la thèse est de développer un modèle macroscopique de cette interaction fluide-structure, afin d'en réaliser des simulations numériques. Une approche numérique de simulation aux grandes échelles est donc mise en place pour capturer les grandes structures de l'écoulement et leur couplage avec les déformations du milieu poreux. Pour cela nous dérivons les équations régissant la grande échelle, au point de vue du fluide ainsi que de la phase solide. À cause du caractère non-local de la phase solide, une approche hybride est proposée. La phase fluide est décrite d'un point de vue Eulerien, tandis que la description de la dynamique de la phase solide nécessite une représentation Lagrangienne. L'interface entre le fluide et le milieu poreux est traitée de manière continue. Cette approche de l'interface fluide/poreux est justifiée par un développement théorique sous forme de bilan de masse et de quantité de mouvement à l'interface. Ce modèle hybride est implémenté dans un solveur écrit en C$++$, à partir d'un solveur fluide disponible dans la librairie CFD \openfoam. Un préalable nécessaire à la réalisation d'un tel modèle macroscopique est la connaissance des phénomènes de la petite échelle en vue de les modéliser. Deux axes sont explorés concernant cet aspect. Le premier consiste à étudier les effets de l'inertie sur la perte de charge en milieu poreux. Un paramètre géométrique est proposé pour caractériser la sensibilité d'une microstructure poreuse à l'inertie de l'écoulement du fluide dans ses pores. L'efficacité de ce paramètre géométrique est validée sur une diversité de microstructures et le caractère général du paramètre est démontré. Une loi asymptotique est ensuite proposée pour modéliser les effets de l'inertie sur la perte de charge, et comprendre comment celle-ci évolue en fonction de la nature de la microstructure du milieu poreux. Le deuxième axe d'étude de la petite échelle consiste à étudier l'effet de l’interaction fluide-structure à l'échelle du pore sur la perte de charge au niveau macroscopique. Comme les cas présentent de grands déplacements de la phase solide, une approche par frontières immergées est proposée. Ainsi deux méthodes numériques sont employées pour appliquer la condition de non-glissement à l'interface fluid/solide: l'une par interface diffuse, l'autre par reconstitution de l'interface. Cela permet une validation croisée des résultats et d'atteindre des temps de calcul acceptables tout en maîtrisant la précision des résultats numériques. Cette étude permet de montrer que l'interaction fluide-structure à l'échelle du pore a un effet considérable sur la perte de charge effective au niveau macroscopique. Des questions fondamentales sont ensuite abordées, telles que la taille d'un élément représentatif ou la forme des équations de transport dans un milieu poreux souple. / We study the turbulent flow of a fluid over a canopy, that we model as a deformable porous medium. This porous medium is more precisely a carpet of fibres that bend under the hydrodynamic load, hence initiating a fluid-structure coupling at the scale of a fibre's height (honami). The objective of the thesis is to develop a macroscopic model of this fluid-structure interaction in order to perform numerical simulations of this process. The volume averaging method is implemented to describe the large scales of the flow and their interaction with the deformable porous medium. An hybrid approach is followed due to the non-local nature of the solid phase; While the large scales of the flow are described within an Eulerian frame by applying the method of volume averaging, a Lagrangian approach is proposed to describe the ensemble of fibres. The interface between the free-flow and the porous medium is handle with a One-Domain- Approach, which we justify with the theoretical development of a mass- and momentum- balance at the fluid/porous interface. This hybrid model is then implemented in a parallel code written in C$++$, based on a fluid- solver available from the \openfoam CFD toolbox. Some preliminary results show the ability of this approach to simulate a honami within a reasonable computational cost. Prior to implementing a macroscopic model, insight into the small-scale is required. Two specific aspects of the small-scale are therefore studied in details; The first development deals with the inertial deviation from Darcy's law. A geometrical parameter is proposed to describe the effect of inertia on Darcy's law, depending on the shape of the microstructure of the porous medium. This topological parameter is shown to efficiently characterize inertia effects on a diversity of tested microstructures. An asymptotic filtration law is then derived from the closure problem arising from the volume averaging method, proposing a new framework to understand the relationship between the effect of inertia on the macroscopic fluid-solid force and the topology of the microstructure of the porous medium. A second research axis is then investigated. As we deal with a deformable porous medium, we study the effect of the pore-scale fluid-structure interaction on the filtration law as the flow within the pores is unsteady, inducing time-dependent fluidstresses on the solid- phase. For that purpose, we implement pore-scale numerical simulations of unsteady flows within deformable pores, focusing for this preliminary study on a model porous medium. Owing to the large displacements of the solid phase, an immersed boundary approach is implemented. Two different numerical methods are compared to apply the no-slip condition at the fluid-solid interface: a diffuse interface approach and a sharp interface approach. The objective is to find the proper method to afford acceptable computational time and a good reliability of the results. The comparison allows a cross-validation of the numerical results, as the two methods compare well for our cases. This numerical campaign shows that the pore-scale deformation has a significant impact on the pressure drop at the macroscopic scale. Some fundamental issues are then discussed, such as the size of a representative computational domain or the form of macroscopic equations to describe the momentum transport within a soft deformable porous medium.

La simulation des procédés de mise en forme des composites à renforts tissés de type RTM est un enjeu majeur pour les industries de pointe mettant en œuvre ce type de matériaux. Au cours de ces procédés, la préforme tissée est souvent soumise à des déformations importantes. La connaissance et la simulation du comportement mécanique de la préforme à l’échelle macroscopique et à l’échelle mésoscopique s’avère souvent nécessaire pour optimiser la phase de conception de pièces composites formées par de tels procédés. Une analyse du comportement mésoscopique des préformes tissées de composites est d’abord proposée. Une loi de comportement hyperélastique isotrope transverse est développée, permettant de décrire le comportement mécanique de chacun des modes de déformation de la mèche : élongation dans la direction des fibres, compaction et distorsion dans le plan d’isotropie de la mèche, cisaillement le long des fibres. Une méthodologie est proposée pour identifier les paramètres de cette loi de comportement à l’aide d’essais sur la mèche et sur le tissu, et une validation par comparaison avec des essais expérimentaux est présentée. Une analyse du comportement macroscopique des renforts interlocks est ensuite proposée : une loi de comportement hyperélastique orthotrope est développée et implémentée. Cette loi, extension de la loi de comportement pour la mèche, est également basée sur une description phénoménologique des modes de déformation de la préforme. Une méthode d’identification des paramètres de cette loi de comportement est mise en œuvre, utilisant des essais expérimentaux classiques dans le contexte des renforts tissés (tension uniaxiale, compression, bias extension test, flexion). Cette seconde loi de comportement est validée par comparaison avec des essais de flexion et d’emboutissage hémisphérique. / Simulating the preforming stage of RTM-like fabric-reinforced composites manufactoring processes is a major stake for industries which use such materials. During such processes, the woven preform often undergoes finite deformations. Simulation methods are then required to optimize the conception of composites parts formed by RTM. An analysis of the mechanical behaviour of woven preforms at mesoscale is first presented. A transversely isotropic hyperelastic behaviour law is developped in order to describe the mechanical behaviour of each deformation mode of the yarn : elongation in the direction of fibres, compaction and distorsion in the transverse plane and along-fibres shear. An identification method is set up for this behaviour law which allows to compute its parameters by use of simple experimental tests on the yarn and on the fabric. The behaviour law is then validated by comparizon between simulations ans experimental tests. An analysis of the mechanical behaviour of interlock woven preforms at macroscale is the presented. An orthotropic hyperelastic behaviour law is developped and implemented as an extension of the behaviour law for the yarn. A phenomenological approach is also used to describe the mechanical behaviour of each deformation mode of the preform. An identification method is set up and put into place, based on tests well known in the field of fabric reinforcements : tensile test, crushing test, bias extension test, flexure test. A hemispherical stamping simulation is set up and compared to experiment for validation purpose.

Afin de pouvoir prédire le comportement des renforts de composites 3D interlock au cours d'un procédé de mise en forme, il est nécessaire de connaitre la position des mèches dans le renfort durant la phase de préformage du procédé. Les travaux présentés ici traitent de la simulation du préformage de renforts 3D épais à l'aide d'un élément fini hexaédrique semi-discret spécifique. En utilisant le principe des travaux virtuels, on distingue le travail interne virtuel dû à la tension des mèches des autres travaux virtuels. La raideur due aux tensions de mèches, qui constitue la contribution principale de la rigidité du matériau, est prise en compte à l'aide de barres incluses dans les éléments. Les rigidités dues aux autres sollicitations, comme la compression transverse, les cisaillements ou les frottements inter-mèches, sont décrites par un matériau continu additionnel. La combinaison de ce modèle discret du premier ordre et d'un matériau continu hyperélastique anisotrope dit du second ordre, pour formuler le comportement du matériau va permettre la simulation du préformage des renforts tissés épais. Conjointement aux travaux sur la simulation, des travaux expérimentaux pour l'identification des paramètres matériau de la loi de comportement ont été définis et réalisés. Ces paramètres concernent les deux parties de la formulation du comportement. / In order to simulate 3D interlock composite reinforcement behavior during forming process, it is necessary to predict yarns positions in the fabric during the preforming stage of the process. The present work deals with thick 3D interlock fabric forming simulation using specific hexahedral semi-discrete finite elements. Using the virtual work principle, we distinguish the virtual internal work due to tensions in yarns from other internal virtual works. The stiffness relative to yarns tension which is the main part of the rigidity is described by bars within the elements. The other rigidities - like transverse compression, shears or friction between yarns - are depicted by a continuous additional material. A combination of this "first order" discrete model and a continuous orthotropic hyperelastic "second order" material formulation will enable us to simulate the interlock preforming process. Jointly to the simulation work, we also had to specify and perform experimental testing identification of material parameters. These parameters concern both parts of the model.

L'industrie aéronautique doit faire face aux nouvelles exigences environnementales, tout particulièrement concernant la réduction de la consommation des énergies fossiles. L'utilisation de matériaux composites plus léger permet de répondre en partie à cette attente. Pour limiter les coûts lors de la fabrication et du développement des composites à renforts tissés 3D, il est nécessaire d'utiliser des outils de simulation performants. En particulier, les outils existants, qui discrétisent à une échelle mésoscopique l'architecture des tissus 3D, ne tiennent pas compte de l'influence du procédé de fabrication sur la constitution de la structure textile. Si des outils numériques dédiés à la modélisation du procédé de tressage et de tricotage sont disponibles, il n'en est rien concernant le tissage. Cette étude avait donc pour but de s'intéresser plus particulièrement à la simulation du prodécé de tissage pour pouvoir obtenir une structure de tissu sèche déformée numériquement. La production de différentes architectures de tissu en verre E dans notre laboratoire nous a permis d'observer les différents éléments en contact avec le fil ou le tissu sur la machine à tisser, par le biais de l'utilisation d'une caméra rapide par exemple. Le développement d'un modèle numérique par éléments finis reproduisant le procédé de tissage a été réalisé. Une loi de comportement isotrope transverse fut utilisée pour modéliser les fils de verre. Des premières simulations numériques encourageantes pour la fabrication d'un tissu d'armure toile et d'un tissu d'armure croisé 2-2 sont présentées et comparées avec les tissus réels produits correspondants. / The aeronautical industry faces new challenges regarding the reduction of fossil fuel consumption. One way to address this issue is to use lighter composite materials. The ability to predict the geometry and the mechanical properties of the unit cell is necessary in order to develop 3D reinforcements in composite materials for these aeronautical applications. There is a difficulty to get realistic geometries for these unit cells due to the complexity of their architecture. Currently, existing tools which model 3D fabrics at a meso scale don't take into account manufacturing process influence on the shape modification of the textile structure. There is already some numerical tools that can model the braiding or knitting process, but none have been developed for weaving so far. Consequently, this study deals with the numerical simulation of the weaving process to obtain a deformed dry fabric structure. During the weaving process of E-glass fabrics, achieved in our laboratory, it has been observed that large deformations led to the modification of transverse section of meshes, or local density changes, that can modify the fabrics mechanical resistance. For this reason, a numerical tool of the weaving process, based on finite element modelling, has been developped to predict these major deformations and their influences on the final textile structure. The correlation between numerical results and fabrics produced with glass fibres has been achieved for plain weave and 2-2 twill.

Dans le cadre du transport de marchandises par voie ferroviaire, un certain nombre de verrous limitent les efforts consentis pour un fret plus "propre". Des améliorations doivent être faites afin de limiter le nombre de locotracteurs assurant les charges/décharges de marchandises en bout de ligne. Dans cette perspective, FEMTO-ST et Alstom Transport ont pour objectif de concevoir et développer un système de gestion d'énergie pour locomotive électrique hybride. Composée d'un groupe électrogène couplé à des accumulateurs électrochimiques et des super-condensateurs, la locomotive électrique hybride permet d'accroître la souplesse et l'efficacité du transport ferroviaire électrique, tout en réduisant encore son impact environnemental. Cette étude a consisté dans un premier temps à développer une modélisation macroscopique de la chaîne de traction électrique hybride et à proposer une structuration de la commande avec l'identification des capteurs matériels et logiciels nécessaires au contrôle optimal de la chaîne de traction. La caractérisation expérimentale des moyens de stockage a permis l'amélioration comportementale et dynamique des modèles correspondant. Dans un second temps, le contrôle de la tension de sortie d'un hacheur dévolteur a permis d'étudier l'application de la logique floue de type-2 (intervalle et générale) dans le cas d'applications industrielles (relativement simples). Enfin, une gestion innovante des flux d'énergie au sein d'un système plus complexe: locomotive électrique hybride, a été développée. Mettant en œuvre la logique floue de type-2 et des algorithmes s'inspirant de la théorie de l'évolution, cette gestion d'énergie optimale opère en temps réel et sans connaissance à priori du cycle de conduite. La gestion fréquentielle ainsi que le contrôle des états de charge des sources secondaires du véhicule contribuent également à leur bon fonctionnement. / To achieve a "greener" freight transport, efforts are still needed to overcome some technological barriers. New improvements must be carried at to limit the shunting locomotives' use intended for the goods' load/unload. Considering this aim, FEMTO-ST and Alstom Transport decided to conceive and develop an energy management strategy system for hybrid electric locomotives. Such locomotives include a diesel driven generator set which is coupled with batteries and ultra-capacitors. The architecture aims at improving the flexibility, the effectiveness of the electrical railway transport and at reducing the environmental impact of these activities again. Firstly, the study consists in implementing a macroscopic modelling of the hybrid electric powertrain. Then, the control is optimally designed by identifying the hardware and software sensors of the powertrain. The dynamics and the behavior of the secondary sources' models are improved thanks to their experimental characterizations. Secondly, the use of Type-2 Fuzzy Logic (interval and general) controllers permits to study their efficiency in the control of a very simple industrial system: the output voltage of a buck converter. Lastly, based on the obtained results, an innovative management of the system's energy flows is developed in the case of the hybrid electrical locomotive. The use of the type-2 fuzzy logic and evolutionary algorithms permit to optimally perform a real time energy management strategy without a priori knowledge of the duty cycle. The frequency approach and the secondary sources' state of charge control contribute to their efficient use.