La grande densité de couple qui fait la réputation des machines à flux transverse (MFT) est attribuable à l’agencement particulier de ses circuits magnétiques et électriques. Cependant, cet agencement complique aussi l’étape de modélisation nécessaire à la conception de la machine. Des facteurs tels que le flux de fuite, la saturation et l’aspect tridimensionnel du flux expliquent en grande partie la difficulté à élaborer un modèle magnétique analytique précis pour la machine. Par ailleurs, l’emploi d’un modèle imprécis pour l’optimisation des caractéristiques de la machine peut conduire à des résultats erronés. Pour cette raison, la plupart des concepteurs fait appel à une modélisation numérique avec des outils de calcul des champs 3-D, coûteux en ressources informatiques, en temps de calcul et malheureusement peu adaptée au processus de conception itératif d’une machine électrique. La thèse défendue dans ce document propose une alternative : elle présente une approche de modélisation hybride, faisant appel à un modèle magnétique analytique assisté par le calcul des champs. L’approche est appliquée pour l’optimisation des performances d’une structure de MFT particulière : la MFT à griffes (ou "claw-pole") à stator hybride (MFTCPSH). Dans un premier temps, on montre comment des facteurs de corrections, dont les valeurs sont déterminées à partir de simulations par éléments finis 3-D, peuvent permettre de palier aux problèmes de précision d’un modèle analytique au cours de l’optimisation de la machine. Un exemple de maximisation du flux à vide d’une MFTCPSH démontre la rapidité, l’efficacité et la précision de l’approche proposée. Dans un second temps, le modèle hybride et la méthode d’optimisation sont développés pour maximiser le couple d’une MFTCPSH. En se basant sur le cahier des charges d’une application de traction électrique de type moteur-roue, il est montré comment le modèle et la méthode proposés permettent, en quelques heures seulement, de déterminer les dimensions géométriques optimales qui maximisent la densité de couple de la machine. / The high power and torque capabilities of Transverse Flux Machines (TFM) are well known. The particular arrangement of the TFM’s electric and magnetic circuits, which mainly explains its exceptional features, also greatly complicates the TFM design process. Building an accurate analytical magnetic model for a TFM is a hard task: factors as the leakage flux, magnetic saturation or 3-D flux paths greatly affect the model accuracy. As the use of an inaccurate model can lead to erroneous results, TFM designers extensively use numerical models, as 3-D finite element analysis (FEA). Despite their accuracy, such numerical tools are still time and resource consuming but also tricky to use for optimization purpose. This thesis proposes an alternative approach for the design of TFM: a hybrid model, based on analytical calculations, corrected with finite element results inside the optimization process. Such an approach is used in order to optimize the performances of a claw-pole TFM with hybrid stator (CTFMHS). This work first presents how FEA-derived correction factors can be inserted in the analytical model to improve the model accuracy, in the course of the design process. The principles and the effectiveness of the proposed approach are shown through an example of CTFMHS no-load flux maximization. The approach is then derived in order to maximize the torque density of the CTFMHS. Considering the specifications of an in-wheel motor application, it is shown how the developed approach can determine, in a few hours only, the optimal machine geometry maximizing its torque density.

Les guides d'ondes optiques sont utilisés dans plusieurs dispositifs de télécommunications pour le transport d'informations numériques. Nous les retrouvons aussi dans les équipements médicaux ou d'ingénierie, comme par exemple les coupleurs, les interféromètres optiques, les convertisseurs de polarisation, etc. Le transport des informations dans ces structures demeure une question capitale et donc les méthodes utilisées doivent être capables de véhiculer ces informations avec un minimum de perte. Dans cette thèse, nous présentons une méthode numérique d'approximation qui permet la simulation du transport des ondes dans une direction choisie avec le maximum de précision et à un coût abordable. Cette méthode est basée sur une approximation au sens de Padé d'opérateurs pseudo-différentiels associés à la fonction algébrique y/1 — x . Essentiellement, nous approchons une équation différentielle ordinaire (EDO) d'ordre 2 par une EDO d'ordre 1 dont les solutions ont la même phase que celle des solutions de l'EDO d'ordre 2 qui se propagent dans la direction choisie. Nous avons évalué les capacités de cette méthode sur de nombreux problèmes bidimensionnels et tridimensionnels. Elle s'est révélée à la fois robuste et précise

Ce projet était composé d’une campagne expérimentale en deux phases. Dans la première phase, des simulations à l’aide de la méthode des éléments finis ont été réalisées afin d'identifier les paramètres les plus sensibles. Dans la seconde phase, une série de 12 poutres en béton armé de 3,4 m de long ont été confectionnées. Ces poutres ont par la suite subi une cure thermique et un séchage à l’air. Ensuite, 10 des 12 poutres confectionnées ont été réparées dans leurs zones centrales en compression à l’aide de deux types de matériaux de réparation. Le premier matériau possédait une plus faible rigidité que le substrat, et le second, une rigidité supérieure. Les conclusions principales de l’étude montrent que le fait d’utiliser un matériau de réparation avec une rigidité supérieure à celle du substrat permet d’obtenir un comportement semblable à celui d’une poutre monolithique. En revanche, l’utilisation d’un matériau à plus faible rigidité tend à affecter négativement le comportement de la poutre réparée. Cette tendance est davantage marquée lorsque l’on considère une charge soutenue, car les matériaux à plus faible rigidité tendent à se déformer davantage que les matériaux plus rigides. Mots clés : Béton, capacité structurale, fluage, compression, MEF, poutre, propriétés mécaniques, réparation / A two-phase research program has been undertaken. In the first phase, computer simulations using a finite-element analysis (FEM) software have been performed to highlight the most sensitive. In the second phase, a series of 12 real-size reinforced concrete beams (3,4 m long) has been cast, air-conditioned and repaired at midspan to different extents, using two formulations of repair self-compacting concrete (high-stiffness and low-stiffness SCC’s), 10 of the 12 beams are repaired beam. The last 2 beams are monolithic beams and used as references. After the installation of the instrumentation, the beams are tested on a 4-point bending setup using two servo-hydraulic actuators, for the measurement of the instantaneous and long-term load deformation. The analysis that follows illustrates the correlation between laboratory results and numerical simulations. The main findings of the study show that using a repair material with a higher stiffness than that of the substrate results in a behavior similar to the monolithic beam. However the use of a material with a lower stiffness tends to negatively affect the behavior of the repaired beam. This trend is more pronounced when considering a sustained load, because the materials with lower stiffness tend to deform more than the rigid materials, it follows that the creep at the young age of the material also affects the behavior of the repaired beams, and that the use of a more rigid material allows the repaired beam to behave like a monolithic beam. Keywords: Beam, concrete, compression, creep, FEM, mechanical repair, structural capacity

Le calcul des surfaces libres a fait l'objet de plusieurs ouvrages. L'objectif de ces travaux est de calculer avec précision les interfaces entre fluides immiscibles en minimisant la diffusion numérique et les pertes de masse. Plusieurs techniques existent pour effectuer de telles simulations et dans ce travail, on présente une variante améliorée de la méthode des surfaces de niveau (« level-set method »). Cette méthode requiert la résolution d'une équation de transport pour l'interface, couplée à la résolution des équations de NavierStokes. Dans la méthode proposée, l'interface est représentée par la surface de niveau 0 de la fonction distance signée. Cette interface présente une zone de transition de largeur 2e des paramètres rhéologiques des deux fluides comme la viscosité et la densité. Ces paramètres sont régularisés à l'interface pour éviter de gérer des discontinuités. La tension superficielle caractérise la résistance à la déformation de l'interface et joue un rôle très important. Il est donc primordial de bien l'imposer comme terme source dans l'équation de Navier-Stokes, ce qui est réalisé par l'introduction d'un tenseur qui agit uniquement à l'interface. Pour remédier au problème de conservation de la masse, nous avons développé une technique qui consiste à modifier la surface de niveau 0 de manière à conserver le volume. Cette correction est calculée avec la méthode de la sécante. Toutes ces stratégies sont facilitées par l'introduction d'une méthode d'adaptation de maillage instationnaire permettant de concentrer les éléments au voisinage de l'interface, là où le besoin s'en fait sentir. Cela assure une imposition très précise des forces surfaciques ainsi que de la transition des caractéristiques rhéologiques des différents fluides. Ce remaillage adaptatif est basé sur un estimateur d'erreur hiérarchique en dimension 2 et sur un estimateur basé sur une métrique en dimension 3. La résolution des équations de Navier-Stokes est faite par une méthode directe en dimension 2 et par une méthode itérative en dimension 3. L'élément de Taylor-Hood et un schéma implicite de différences arrières d'ordre 2 sont utilisés pour les discrétisations en espace et en temps. On utilise une méthode SUPG pour les équations de transport et de régularisation de l'interface. Enfin, plusieurs problèmes sont traités afin de valider la méthode et de comparer nos résultats avec ceux obtenus par d'autres chercheurs et différentes méthodes. On soulève enfin un certain nombre de questions sur les méthodes employées dans la littérature.

La présente thèse présente le développement d'un premier modèle d'électrode poreuse pour une cathode Hall-Héroult ainsi que le développement de nouvelles méthodes numériques. Les nouvelles méthodes numériques sont l'implantation du critère d'électroneutralité par contrainte de type multiplicateur de Lagrange ainsi que l'application d'un équilibre thermodynamique ionique par la méthode de pénalisation. Une revue des divers processus physico-chimiques associés à la cathode est présentée. Le modèle d'électrode poreuse est développé avec une formulation ionique pour le bain comportant cinq espèces. Le modèle est développé avec une méthode d'homogénéisation pour le volume poreux. Les lois de transport sont formulées en fonction des concentrations ioniques. Les concentrations ioniques d'équilibre ont été calculées par un modèle thermodynamique simple formulé selon des données thermodynamiques empiriques. Le modèle prend en compte la migration, la diffusion et la convection des espèces ioniques. Les divers paramètres de transport ont été calés avec des mesures empiriques de conductivité conjointement aux concentrations ioniques d'équilibre. Le modèle prend en compte la réaction électrochimique de production d'aluminium au niveau des pores du carbone. De plus, un équilibre ionique est imposé entre certaines espèces. La méthode par éléments finis a été utilisée pour effectuer les calculs numériques. De nouvelles méthodes numériques d'implantation du critère d'électroneutralité et d'un équilibre ionique sont présentées. Les méthodes ont été utilisées avec succès et permettent d'éviter d'éliminer une équation de conservation. La méthode de pénalisation s'est avérée comparable aux méthodes d'élimination par le biais d'un terme source mais les résultats sont partiellement différents du modèle principal procédant par élimination d'équations de conservation. Certaines méthodes par élimination négligent en partie la contribution au transport provenant des termes associés aux équations de conservation éliminées. Par la suite, des résultats de simulation du début de l'électrolyse sont montrés. Le modèle donne des résultats conformes aux ordres de grandeur des mesures expérimentales.

L’objectif principal de cette thèse est le développement d’un outil numérique, en utilisant l’approche XFEM, permettant la simulation des problèmes transitoires non linéaires thermoélectriques dans un milieu fissuré en deux dimensions, avec prise en compte des échanges thermiques et électriques entre les lèvres de la fissure. La simulation numérique de la propagation de fissures présente un grand intérêt pour de nombreux secteurs industriels (production d’aluminium, aéronautique, nucléaire, etc.). De plus, c’est un problème complexe sur le plan numérique. La méthode d’éléments finis classiques présentent des contraintes importantes de raffinement de maillage en fond de fissure, de remaillage pendant la propagation de la fissure avec la projection des champs, ce qui a pour effet d’augmenter le temps de calcul et de dégrader la précision des résultats. D’autre part, la méthode des éléments finis étendue XFEM, a reçu un succès grandissant pour le traitement de problèmes avec fissures durant la dernière quinzaine d’années. Elle permet d’utiliser un maillage qui ne se conforme pas à la géométrie des fissures, ceci grâce à un enrichissement de l’approximation éléments finis. Dans cette thèse, on s’intéresse à étendre le champ d’application de la méthode XFEM pour les problèmes non linéaires thermoélectriques avec fissures. En effet, le problème thermique transitoire est couplé avec le problème électrique par la génération de la chaleur dans le solide, et la génération de chaleur à la fissure à cause de la résistance de l’interface. Les échanges thermiques et électriques entre les lèvres de la fissure sont aussi considérés, et dépendent, respectivement, du saut de la température et du potentiel électrique à la fissure. En raison de la génération de la chaleur dans le solide et aux lèvres L’objectif principal de cette thèse est le développement d’un outil numérique, en utilisant l’approche XFEM, permettant la simulation des problèmes transitoires non linéaires thermoélectriques dans un milieu fissuré en deux dimensions, avec prise en compte des échanges thermiques et électriques entre les lèvres de la fissure. La simulation numérique de la propagation de fissures présente un grand intérêt pour de nombreux secteurs industriels (production d’aluminium, aéronautique, nucléaire, etc.). De plus, c’est un problème complexe sur le plan numérique. La méthode d’éléments finis classiques présentent des contraintes importantes de raffinement de maillage en fond de fissure, de remaillage pendant la propagation de la fissure avec la projection des champs, ce qui a pour effet d’augmenter le temps de calcul et de dégrader la précision des résultats. D’autre part, la méthode des éléments finis étendue XFEM, a reçu un succès grandissant pour le traitement de problèmes avec fissures durant la dernière quinzaine d’années. Elle permet d’utiliser un maillage qui ne se conforme pas à la géométrie des fissures, ceci grâce à un enrichissement de l’approximation éléments finis. Dans cette thèse, on s’intéresse à étendre le champ d’application de la méthode XFEM pour les problèmes non linéaires thermoélectriques avec fissures. En effet, le problème thermique transitoire est couplé avec le problème électrique par la génération de la chaleur dans le solide, et la génération de chaleur à la fissure à cause de la résistance de l’interface. Les échanges thermiques et électriques entre les lèvres de la fissure sont aussi considérés, et dépendent, respectivement, du saut de la température et du potentiel électrique à la fissure. En raison de la génération de la chaleur dans le solide et aux lèvres / The main objective of this thesis is the development of a numerical tool, using the XFEM approach, for the simulation of transient nonlinear thermoelectrical problems in fractured media in two dimensions, taking into account thermal and electrical exchanges between the crack’s lips. Numerical simulations of crack propagation are of great interest for many industrial sectors (aluminum production, aerospace, nuclear, etc.). In addition, this is a numerically complex problem. The classical finite element method has important constraints of mesh refinement at the crack tip, remeshing during crack propagation and field projections, which has the effect of increasing the computation time and degrading the accuracy. On the other hand, the eXtended Finite Element Method (XFEM), has received a growing success for the treatment of the problems containing cracks in the last fifteen years. It allows using a mesh that does not conform to the geometry of the crack; this is possible by the enrichment of the finite element approximation. In this thesis, we are interested in extending application field of the XFEM method to the nonlinear thermoelectrical problems with cracks. Indeed, the transient thermal problem is coupled to the electrical problem by the heat generation in the solid, and the heat generation at the crack’s lips due to the interface resistance. The heat and electrical exchanges between the crack’s lips are also considered, and depend, respectively, on the temperature and the voltage jump at the crack. Due to the heat generation in the solid and in crack’s lips (Joule effect), and the temperature dependence of the physical parameters of the material, the problem is nonlinear and fully coupled. The discretized nonlinear system by the XFEM method is solved using the Newton-Raphson algorithm. The robustness of the proposed technique is demonstrated through the simulation of different examples, and the results shows an excellent agreement with the analytical solution, or with the finite element solution using a refined mesh.

Les processus de transport et d'atténuation des contaminants dans un sol saturé sont déterminés en fonction de plusieurs facteurs parmi lesquels la consolidation à grande déformation joue un rôle majeur. L'étude de la consolidation d’un sol saturé tient généralement compte des aspects mécaniques et hydrauliques du milieu. Le facteur qui induit une modification structurale importante d'un sol saturé recevant une charge externe est l'expulsion de l'eau dans les pores. L'objectif de cette recherche vise à intégrer les équations de transport des contaminants avec les équations générales de la consolidation à grande déformation des sols (Foriero et Ladanyi, 1995 et 1998). Il est important pour le spécialiste en sols et en environnement de connaître l’état des contraintes, la perméabilité, l’indice des vides, le degré de consolidation, mais aussi les concentrations de contaminants et leurs gradients en vue de faire des prévisions dans le cadre d’un projet affectant un sol contaminé ou à d’autres fins. Pour ce faire, un programme informatique basé sur la méthode des éléments finis a été développé pour résoudre l'équation différentielle résultant du couplage de la consolidation et de la contamination. Les variables utilisées dans le logiciel sont les propriétés du sol, le coefficient de dispersion et la concentration initiale du contaminant. Les résultats de recherche indiquent que le programme informatique peut être utilisé pour prédire les paramètres reliés à la consolidation à grande déformation du sol et à l’évolution de la concentration du contaminant.

Les propriétés physiques, mécaniques et chimiques des matières premières ont un effet majeur sur la qualité des anodes en carbone pour le procédé de production d’aluminium. Ce travail tente d’étudier la faisabilité de l’application de simulation de la Méthode des Élément Discrets (DEM) à la technologie de production d’anodes. L’effet de la forme des particules et de la distribution de leurs tailles sur la densité apparente vibrée (VBD) d’échantillons de coke sec est étudié. Les particules de coke sont numérisées en utilisant des techniques d’imagerie à deux et trois dimensions. Ces images donnent les formes et les aspects réels des particules qui sont utilisées pour les modèles DEM pour les tests VBD pour le codage de flux de particules (PFC). Le coefficient de friction interne des particules de coke est estimé par la méthode de mesure d’angle au repos. Les résultats ont montrés comme attendu, que la VBD des échantillons de coke est affectée par la forme et la distribution de taille des particules. Les simulations à deux dimensions ont confirmé qu’en général, les échantillons formés de particules de tailles poly-dispersées ont une VBD plus haute que ceux dont la taille des particules est mono-dispersée. De plus, la VBD des échantillons augmente lorsque la fraction de grosses particules augmente. Cependant, la présence de 10 % massique de particules fines est nécessaire pour remplir les pores entre les grosses particules. De même pour la simulation 3D, le modèle suit la tendance des données expérimentales montrant que dans une éprouvette de 2,9 cm de diamètre, l’augmentation de la quantité de particules de - 4+6 mesh (de 3,36 à 4,76 mm) engendre une augmentation de la VBD. En conclusion, un modèle DEM approprié est capable de prédire le réarrangement des particules et l’évolution de la densité pendant le processus de vibration. / Physical, mechanical and chemical properties of raw materials have considerable effects on quality of carbon anodes for aluminium smelting process. The present work attempts to investigate the feasibility of application of Discrete Element Method (DEM) simulations in anode production technology. Effects of coke particles shape and size distribution on vibrated bulk density (VBD) of dry coke samples are studied. Coke particles are digitized using two-dimensional and three-dimensional imaging techniques and real-shape particles are used in DEM models of VBD test in Particle Flow Code (PFC). Internal friction coefficient of coke particles were estimated by means of angle of repose tests. Results showed that, as expected, VBD of coke samples is affected by shape and size distribution of the particles. Two-dimensional simulations confirmed that in general, mixed-sized samples have higher VBD than mono-sized cokes and as the fraction of coarse particles increases vibrated bulk density increases. However, existence of 10 wt.% of fine particles to fill the pores between coarse particles is essential. For 3D simulations also, the model follows the trend of experimental data showing that in the container of 2.9 mm diameter, as the content of -4+6 mesh (3.36-4.76 mm) particles increase, VBD increases. It can be concluded that a well-tailored DEM model is capable of predicting the particle rearrangement and density evolution during the vibration process.

L'objectif des travaux présentés dans cette thèse est d'explorer de nouvelles méthodes numériques efficaces et robustes pour la résolution par la méthode des éléments finis des problèmes d'écoulements des fluides viscoélastiques pour des nombres de Weissenberg élevés. La motivation de cette étude provient en grande partie au fait que les fluides viscoélastiques sont à la base de nombreuses applications industrielles, notamment dans l'industrie des polymères, l'industrie des pâtes de papiers, l'industrie alimentaire, etc... Ces écoulements sont régis par un système d'équations aux dérivées partielles fortement non linéaire composé des équations de Navier Stokes (ou Stokes pour les écoulements lents) et d'une équation constitutive modélisant la contribution du polymère. Les rhéologues ont mis à la disposition de la communauté scientifique plusieurs modèles. Dans ce travail, afin de pouvoir effectuer des comparaisons avec les travaux effectués par d'autres chercheurs déjà disponibles dans la vaste littérature du domaine, nos simulations numériques ont été faites pour les modèles de Oldroyd-B, PTT et Giesekus. La méthodologie utilisée est toutefois assez générale. Une des difficultés majeures rencontrées dans la simulation numérique des fluides viscoélastiques est la perte de convergence des algorithmes pour des problèmes à faibles nombres de Weissenberg. Ainsi, dans l'optique d'améliorer la convergence des schémas numériques, nous nous sommes proposés de développer et d'implémenter de nouveaux algorithmes en nous basant sur la nouvelle formulation en logarithmique du tenseur de conformation. Cette formulation a été utilisée pour la première fois dans le contexte des fluides viscoélastiques par Fattal et Kupferman et al. (21). Par la suite, en nous inspirant de Coronado et al. (35), nous avons introduit deux autres formulations logarithmiques. La première est une variante de celle proposée par Coronado et al. (35) avec un traitement spécial pour le terme convectif. La deuxième formulation permet une linéarisation complète des équations pour la résolution par la méthode de Newton. Les deux méthodes que nous proposons ont l'avantage d'être relativement simple à mettre en oeuvre. Pour tester la robustesse de ces nouvelles formulations, quelques calculs ont été effectués pour le problème d'écoulement autour d'un cylindre et l'inévitable problème de la contraction 4:1 sur des maillages non structurés anisotropes obtenus par la stratégie générale de remaillage adaptatif développée au GIREF. Nous avons ensuite jugé utile d'appliquer cette nouvelle formulation logarithmique pour la résolution de quelques problèmes de surfaces libres. Nous avons étudié la présence des écoulements secondaires et leur influence sur la morphologie des écoulements complexes des problèmes multi-couches dans une filière. En effet, nous avons montré qu'il possible de résoudre des problèmes multi-couches dans une filière tridimensionnelle en résolvant un problème bidimensionnel de surface libre instationnaire seulement sur une section. Nous avons également traité un problème de cisaillement d'une goutte dans une matrice fluide. Il s'agit ici d'un écoulement régi par la tension superficielle. Le calcul des sufaces libres est rendu possible grâce à la méthode des surfaces de niveau (level set method). Cette méthode requiert la résolution d'une équation de transport pour trouver la nouvelle position de l'interface. L'interface constitue une zone de transition de largeur 2 ? des paramètres rhéologiques des deux fluides. Ces paramètres sont régularisés à l'interface pour éviter de gérer des discontinuités. Une technique de réinitialisation a été mise en oeuvre pour remédier à d'éventuels problèmes de conservation de masse. Ces techniques sont falicitées par l'introduction de la méthode de remaillage adaptatif. Celle-ci permet de concentrer les éléments dans l'interface permettant ainsi de faire les calculs avec une meilleure précision. Pour la résolution des équations hyperboliques, notamment l'équation constitutive, de transport et de régularisation de l'interface, nous avons utilisé la méthode de stabilisation de type SUPG. Il résulte de ces développements une méthode de résolution nettement plus performante que les méthodes classiques.

Le procédé d'extrusion permet d'obtenir des tubes de section circulaire à partir d'une billette d'aluminium. Cependant, l'extrusion tubulaire de l'alliage 6063 lui confère de faibles propriétés mécaniques ainsi que des épaisseurs de paroi limitées. Elle donne aussi de mauvaises tolérances ainsi qu'un fini de surface qui laisse à désirer. L'opération d'étirage subséquente permet alors de réduire le diamètre ainsi que l'épaisseur de la paroi des tubes afin d'améliorer leurs tolérances, d'augmenter leurs propriétés mécaniques par écrouissage ainsi que d'obtenir un meilleur fini de surface. Pom atteindre la condition T832 prescrite par le client, un certain pourcentage de réduction est nécessaire afin d'écrouir les tubes suffisamment. Le pourcentage de réduction à atteindre est de 50% afin d'obtenir les bonnes propriétés mécaniques. Actuellement, ce haut pourcentage de réduction nécessite deux passes d'étirage afin d'éviter le bris des tubes. Le projet vise à optimiser les paramètres du procédé d'étirage à froid afin de permettre la réduction nécessaire en une seule passe. L'optimisation de la géométrie de l'outillage est effectuée par le biais d'un modèle éléments finis. Une loi de comportement élasto plastique déjà implémentée dans le logiciel commercial LS-DYNA a été utilisée. Le comportement aux interfaces de contact du tube avec l'outillage est modélisé par un contact surface sur surface avec un coefficient de friction de Coulomb. Les deux variables étudiées pour la géométrie des filières sont l'angle d'entrée ainsi que le rayon de courbure entre la région conique et la portée cylindrique. Des expérimentations ont été effectuées afin de valider les résultats numériques. Une bonne corrélation entre les résultats numériques et expérimentaux a été trouvée.