Nous présentons deux approches différentes pour décrire le couplage hydromécanique des géomatériaux. Dans une approche de type phénoménologique nous traitons le milieu poreux comme un milieu continu équivalent dont les interactions entre la phase fluide et le squelette solide constituent le couplage du mélange à l'échelle macroscopique. En caractérisant le comportement de chaque phase nous arrivons à décrire le comportement couplé du milieu couplé saturé.Nous utilisons cette approche pour modéliser des essais expérimentaux faits sur un cylindre creux pour une roche argileuse (argile de Boom). Les résultats expérimentaux montrent de façon claire que le comportement de cette roche est fortement anisotrope. Nous avons choisi de modéliser ces essais en utilisant une lois de comportement élasto-plastique pour laquelle la partie élastique est transversalement isotrope.Le problème aux conditions aux limites étudié met en évidence des déformations localisées autour du forage intérieur. Afin de décrire de façon objective le développement de ces bandes de cisaillement nous avons considéré un milieu continu local de type second gradient qui permet d'introduire une longueur interne. De ce fait nous avons pu étudier le problème d'unicité en montrant qu'un changement de la discrétisation temporelle du problème aux limites peut conduire à des solutions différentes.Dans la deuxième approche étudiée nous caractérisons la microstructure du matériau avec des grains et un réseau de canaux pour la phase fluide. À l'aide d'un processus numérique d'homogénéisation nous arrivons à calculer numériquement la contrainte du mélange et le flux massique. Cette méthode d'homogénéisation numérique a été implémentée dans un code aux éléments finis afin d'obtenir des résultats macro. Une validation de l'implentation est proposée pour des calculs en mecanique pure et en hydromécanique.

Les micronageurs et parmi eux les microangeurs artificiels sont en général, limités à exister dans des écoulements dominés par des forces visqueses. Ces écoulements sont caractérisés par un bas nombre de Reynolds (Re). Cela impacte la stratégie de nage et plus particulièrement les séquences de forme possibles, qui doivent nécessairement être non-réciproques dans l'espace de déformation pour espérer induire un déplacement net non-nul. De plus, due aux forts effets de traînée, les vitesses de nage sont limités à des valeurs faibles.Dans cette thèse, on examine la possibilité d'utiliser un mécanisme de nage basé sur l'instabilité de flambage d'une sphère creuse. Cette instabilité est provoquée en soumettant la sphère à une onde de pression. La particularité de ce mécanisme est qu'il satisfait par construction la condition nécessaire de nage à bas Reynolds exposée précédemment. De plus, la rapidité de la déformation lors de l'instabilité pousse à prévoir l'apparition d'effets inertiels, et ce même à l'échelle microscopique.Une étude expérimentale a été conduite à l'échelle macroscopique dans le but de comprendre la dynamique de l'instabilité et son impact sur le fluide qui entoure la coque creuse. Ces expériences nous permettent de montrer qu'un déplacement net non-nul est produit pour tous les régimes d'écoulements.On met en évidence le rôle de paramètres géométriques, des propriétés du matériau composant la coque creuse et de la rhéologie du fluide sur l'efficacité de la nage.On montre l'existence d'un optimum de déplacement net pour des valeurs intermédiaires du nombre de Reynolds. Pour expliquer cela, on se sert de mesures de PIV résolues temporellement pour mettre en évidence la présence d'effets d'histoire non-triviaux qui augmentent le déplacement net.On dérive un simple modèle en se basant sur les observations expérimentales pour montrer que ce régime optimal de nage est atteignable pour des sphères microscopiques, ceci est possible grâce l'activation rapide de l'instabilité. Cette propriété permet aussi une excitation à haute fréquence en utilisant des ultrasons. Une étude d'échelle nous permet de prédire une vitesse de nage de 1 cm/s pour un micro-robot contrôlé à distance. Cet ordre de grandeur de vitesse est idéal pour des applications biologiques comme la distribution ciblée de médicaments. / Microswimmers, and among them aspirant microrobots, are generally bound to cope with flows where viscous forces are dominant, characterized by a low Reynolds number (Re). This implies constraints on the possible sequences of body motion, which have to be nonreciprocal. Furthermore, the presence of a strong drag limits the range of resulting velocities.Here, we propose a swimming mechanism which uses the buckling instability triggered by pressure waves to propel a spherical hollow shell. The particularity of this mechanism is that it fulfills naturally the necessary condition of swimming at low Re. In addition, the swiftness of the instability might produce inertial effects even at the microscopic scale.With a macroscopic experimental model we show that a net displacement is produced at all Re regimes. We put in evidence the role of geometrical parameters, shell material properties and rheology of the surrounding fluid on the swimming efficiency.An optimal displacement is reached at intermediate Re. Using time-resolved PIV measurements, we explain that non-trivial history effects take place during the instability and enhance net displacement.Using a simple model, derived from the study of shell dynamics, we show that due to the fast activation induced by the instability, this regime is reachable by microscopic shells. The rapid dynamics would also allow high frequency excitation with standard traveling ultrasonic waves. Scale considerations predict a swimming velocity of order 1 cm/s for a remote controlled microrobot, a suitable value for biological applications such as drug delivery.

Dans le cadre d’une Cifre avec Airbus Helicopters, le projet a pour but le développement d’un modèle pour le pré-dimensionnement de pièces lamifiées élastomère-métal dont le rôle est critique en termes de conception et de sécurité pour les architectures de rotors d’hélicoptères. Pour cela, un premier volet de la thèse a consisté à caractériser le comportement élasto-dissipatif du matériau d’étude (un butadiène chargé de noir de carbone) via une campagne d’essais statiques et dynamiques, sous différents modes de sollicitations (uniaxiales et biaxiales) et sur un spectre assez large de fréquences, d’amplitudes et de températures. A partir de ces observations, un modèle phénoménologique de comportement hyper-viscoélastique est proposé. Sur la base d’un modèle de Maxwell généralisé, il permet de traduire les phénomènes observés sur la gamme de sollicitations visées. Un accent particulier a été porté sur la prise en compte de l’effet Payne en adoptant une approche originale. Les paramètres du modèle sont identifiés par une méthode robuste et rapide. Le modèle est ensuite développé à la fois dans un code commercial de calcul par éléments finis et dans un outil de calcul basé sur une méthode de réduction de modèles. Enfin, une étude du comportement en fatigue est réalisée à travers une campagne d’essais originaux servant de point de départ à la proposition d’une loi d’endommagement continu. / In association with Airbus Helicopters, the aim of the project is to develop a model to pre-size elastomer-metal laminated devices whose role is critical in terms of design and safety for helicopters rotor architectures. To do so, the first part of this thesis consisted in characterizing the elasto-dissipatice behavior of the studied material (a carbon black filled butadiene rubber) thanks to static and dynamic tests, with several kind of loading (uni-axial and bi-axial) and over a wide range of frequences, amplitudes and temperatures. From these observations, a phenomenological hyper-viscoelastic model is proposed. Based on a generalized Maxwell model, it is able to describe the phenomena over the loading range of concern. A particular focus is made to take in account the Payne effect thanks to an original approach. The model parameters are identified with a fast and robust method. The model is then implemented in a commercial finite element code and in a tool based on a model reduction method. Last, a study of the behaviour in fatigue is performed with an original characterization campaign from which a continuous damage law is proposed.

L'objectif de ce travail de thèse est d'étudier analytiquement et numériquement le comportementdynamique non-linéaire des rotors, en prenant en compte des effets significatifs comme les grandesdéformations en flexion, les non-linéarités géométriques et le cisaillement. Le manuscrit est diviséen trois parties principales. Dans la première partie, le principe de Hamilton est utilisé pour formulerles équations du mouvement qui prennent en compte un ensemble d’effets non-linéaires comme desdéformations d'ordre supérieur en flexion et le cisaillement. De plus, si les supports du rotor nepermettent pas à l'arbre de se déplacer dans la direction axiale, il y a alors une force dynamiqueharmonique agissant axialement sur le rotor en fonctionnement. Ces modèles se composentd’équations différentielles non-linéaires du deuxième et du quatrième ordre.Les deux parties suivantes sont consacrées à la résolution des différents modèles non-linéairesdéveloppés dans la première partie. Des méthodes analytiques et numériques sont appliquées afin detraiter les équations non-linéaires du mouvement. Une méthode basée sur des développementsasymptotiques, la méthode des échelles multiples (MEM) est utilisée. Les courbes de réponse sonttracées pour différentes résonances possibles et l'effet de la non-linéarité est discuté par rapport àl'analyse linéaire. La réponse forcée du système provoquée par un balourd est également présentéepour plusieurs configurations du rotor. Lorsque les déformations de cisaillement sont prises encompte, l'analyse est effectuée pour différents élancements afin de mettre en évidence cet effet sur ladynamique d’un système arbre-disque / The objective of the present work is to investigate the nonlinear dynamic behavior of the rotorsystems analytically and numerically, taking into account the significant effects, for example, higherorder large deformations in bending, geometric nonlinearity and shear effects.This thesis is dividedinto two major parts. In the first part, Hamilton’s principle is used to derive the equations of motionwhich take into account various effects, for example, nonlinearity due to higher order largedeformations in bending and shear effects. In addition, if the supports of the rotor do not allow theshaft to move in the axial direction, then there will be a dynamical force acting axially on the rotoras it operates. The mathematical models are composed of coupled nonlinear differential equations ofthe 2nd and the 4th order.In the second part, the resolution of various nonlinear models developed in the first part isaddressed. Analytical and numerical methods are applied for treating the nonlinear equations ofmotion. A method based on asymptotic developments, the method of multiple scales (MMS) is used.The response curves are plotted for different possible resonance conditions and the effect ofnonlinearity is discussed with respect to the linear analysis. The forced response of the system due toa mass unbalance is also presented for various configurations of the rotor. When shear deformationsare taken into account, the analysis is performed for various slenderness ratios to highlight sheareffects on the dynamics of the shaft-disk rotor systems

Nous analysons les propriétés dynamiques de milieux architecturés périodiques et de réseaux fibreux aléatoires en petites et grandes déformations, à partie de méthodes d’homogénéisation afin de calculer leurs propriétés statiques et dynamiques. Des modèles effectifs de type micropolaire et du second gradient sont élaborés afin de prendre en compte l’impact de la microstructure sur le comportement effectif. L’influence des degrés de liberté en rotation additionnels et des gradients d’ordre supérieur du déplacement sur les relations de dispersion sont analysés pour des comportements élastique et viscoélastique du matériau constitutif. Les milieux continus généralisés ainsi construits conduisent à des effets dispersifs, en accord avec les observations. Dans la seconde partie du travail, nous analysons l’influence des grandes déformations sur la propagation des ondes élastiques dans des milieux architecturés périodiques. Des méthodes théoriques assortis de schémas numériques sont développés afin de prédire l’influence des déformations finies générées au sein des structures sur l’évolution de leur diagramme de bande. Un schéma incrémental d’évolution de la fréquence et de la vitesse de phase du milieu continu homogénéisé est établi, à partir d’une méthode de perturbation établie pour des structures 1D, 2D et 3D, en considérant plus particulièrement des structures auxétiques. Ce schéma montre un effet important de l’état de déformation appliquée et de la densité effective sur l’évolution de la fréquence et de la vitesse de phase des ondes. Une méthode de perturbation spécifique aux structures périodiques nonlinéaires est développée afin de généraliser le théorème de Bloch pour couvrir les non linéarités tant géométriques que matérielles. Des modèles hyperélastiques du premier et du second gradient de différentes structures sont identifiés par des tests virtuels reposant sur une méthode d’homogénéisation dédiée, qui permettent de formuler des équations d’onde spécifiques – équations de Burgers et de Boussinesq – dont les propriétés dispersives sont analysées / Micropolar and second gradient effective continua are constructed as two different strategies to account for microstructural effects. The influence of additional degrees of freedom or higher order displacement gradients on the dispersion relations is analyzed in both situations of elastic and viscoelastic behaviors of the material. Generalized effective continua lead to dispersive waves, as observed in experiments. In the second part of the thesis, we analyze the influence of large deformations on the propagation of acoustic waves in repetitive network materials. Both theoretical and numerical methods are developed in order to assess the influence of finite strains developing within the networks on the evolution of their band diagrams. An incremental scheme for the update of frequency and phase velocity of the computed homogenized medium is developed based on a perturbation method for 1D, 2D and 3D structures, considering with a special emphasis auxetic networks. This scheme shows an important effect of the applied finite deformation on the frequency and phase velocity of the propagating waves. A perturbation method for nonlinear periodic structures is developed to extend Bloch’s theorem to cover both geometrical and material nonlinearities. Hyperelastic first and second order gradient constitutive models of different network materials are identified based on dedicated homogenization methods, from which specific wave equations are formulated - Burgers and Boussinesq equations - the dispersion properties of which are analyzed

La striction diffuse et localisée sont des phénomènes précurseurs à la rupture ductile et représentent l'une des principales causes de rebut des pièces métalliques au cours de leur mise en forme. La mise en œuvre d'outils théoriques et numériques capables de prédire l'apparition de ces défauts s'avère nécessaire pour des raisons économiques et environnementales. Ces outils nécessitent en partie l'intégration d'un modèle de comportement adéquat permettant de reproduire les phénomènes physiques mis en jeu. Un tel modèle de comportement est ensuite couplé à un indicateur d'instabilité plastique offrant la possibilité de prédire de manière fiable les phénomènes de striction diffuse et localisée. Dans ce travail de thèse, nous avons considéré certains modèles d'endommagement micromécaniques basés sur l'approche de Gurson, et qui ont été couplés à différents critères d'instabilités plastiques reposant sur l'Analyse de Bifurcation. L'implantation numérique des modèles retenus a été réalisée dans le code de calcul par éléments finis Abaqus/Standard. En ce qui concerne les modèles d'endommagement, plusieurs schémas d'intégrations ont été testés afin d'analyser leur performance et leur robustesse lorsque le comportement présente un effet adoucissant. L'approche combinant le modèle de comportement et les critères de striction a été utilisée pour prédire les limites de formabilité en striction diffuse et localisée de plusieurs matériaux métalliques. Les résultats obtenus ont permis de dresser une hiérarchisation théorique et numérique des critères de striction utilisés. / Diffuse and localized necking are precursors to ductile fracture, and represent one of the main causes of metal parts rejection during forming operations. The implementation of theoretical and computational tools to predict the occurrence of these defects turns out to be necessary for economic and environmental reasons. These tools require in part the introduction of an appropriate behavior model in order to reproduce the physical phenomena involved during forming operations. Such a behavior model is then coupled to a plastic instability indicator providing the ability to reliably predict diffuse and localized necking. In the present work, we considered a micromechanical damage models based on Gurson's approach, which were coupled to different plastic instabilities criteria, based on Bifurcation Analysis. The numerical implementation of these models was carried within the implicit finite element code Abaqus/Standard. With regard to the damage models, several integration schemes were tested to analyze their performance and robustness when the behavior exhibits softening effect. The approach combining the Gurson's damage model and necking criteria has been applied for the prediction of formability limits of several metallic materials. The obtained results allowed establishing a theoretical and numerical classification between the necking criteria used in this work.

La dégénérescence des disques intervertébraux (DD) est un processus naturel qui touche une grande partie de la population. Toutefois, elle peut devenir pathologique et ainsi être accélérée par de multiples facteurs et conduire à un dysfonctionnement précoce des disques lombaires. La DD est un phénomène asymptomatique, ce qui le rend délicat à diagnostiquer précocement. Dans ce cadre, le présent travail vise à développer une méthodologie de diagnostic de la DD utilisant imagerie médicale et modélisation biomécanique.Les disques intervertébraux (DIV) sont des articulations fibro-cartilagineuses qui lient les vertèbres du rachis entre elles. Ils ont pour rôle de supporter et redistribuer les chargements appliqués au rachis tout en assurant sa mobilité. Le DIV peut être assimilé d’un point de vue mécanique à un milieu poreux biologique dont le squelette est formé par la matrice extracellulaire (MEC) et la phase fluide par son contenu hydrique. C’est un milieu hétérogène et anisotrope. Possédant une pression osmotique interne élevée, les DIV ne sont pas vascularisés et possèdent un contenu hydrique important dans lesquels les nutriments cellulaires diffusent à partir des vertèbres adjacentes. Cette particularité rend le fonctionnement des disques très dépendant de la régulation en eau en particulier lorsque ces derniers sont soumis à des sollicitations mécaniques non physiologiques.L’une des clés dans l’élaboration d’un diagnostic fiable de la DD peut être de déterminer de manière objective certaines des caractéristiques internes au disque relatives à son hydratation. Dans le présent travail, une méthodologie de diagnostic quantitatif de la viabilité discale et plus largement de la DD est ainsi proposée. Cette méthodologie se base sur un couplage entre l’imagerie médicale et l’analyse par modélisation du comportement biomécanique du DIV. Un modèle physique biphasique du DIV est développé à l’aide de données déduites d’images IRM. Ce modèle tient compte de l’anisotropie du DIV et des effets induits par les déformations mécaniques sur le processus de transport nutritif.La partie mécanique du modèle est tout d’abord validée en utilisant des résultats d’essais de relaxation réalisés précédemment par l’équipe de biomécanique d’IRPHE. Le modèle est à même de reproduire la distribution de la porosité à l’intérieur du DIV et la perte d’eau suite à un effort mécanique. Cela permet de confirmer la fiabilité de la méthodologie de diagnostic développée pour les cas étudiés. En intégrant le couplage mécano-métabolique, la réponse biomécanique du DIV à des chargements quasi statiques et cycliques est ensuite analysée. / Intervertebral disc degeneration (DD) is a natural process affecting a large part of the population. However, it can be accelerated by several factors and then become associated to lumbar disc pathologies. Given that the DD is symptomless, it is complicated to diagnostic it early. In this context, the present work aims to develop a new methodology of DD diagnostic based on medical imaging and biomechanical modeling.The intervertebral discs (IVD) are fibro-cartilaginous joints that connect vertebrae together ensuring their relative motion. They support and distribute loads applied on the spine. Mechanically, the IVD can be assimilated to a biological porous media in which the solid phase is composed of the extracellular matrix (MEC) and the fluid phase is formed by water content. It is a non-homogeneous and anisotropic component. The IVD is nonvascularized and characterized by a high water content in which cell nutrients diffuses from the adjacent vertebrae. Given this latter particularity, the porosity of the IVD presents a key factor in in its functioning especially when exposed to mechanical non-physiological loads. One of the main keys in the elaboration of a reliable DD diagnostic is to determine objectively some of the IVD interne properties related to its water content. In the present work, a new quantitative methodology of DD diagnostic is proposed. This method is based on a coupling between the medical imaging and the IVD biomechanical behavior modeling. A biphasic IVD model is developed using data from MRI imaging. This model takes into account the anisotropy of the IVD and the effect of the mechanical deformation on the nutrient transport process.The mechanical part of the developed model was firstly validated using results of relaxation tests previously performed by the biomechanics staff of IRPHE. The model is able to predict water distribution within the IVD and the water loss as response to a mechanical load. These results confirm the reliability of the developed diagnostic methodology for the studied cases. The biomechanical response to quasi-static and cyclic loads is then analyzed. The present work proposes and evaluates a new methodology of DD quantitative diagnostic. It also analyzes and compares mechanical and metabolic responses to quasi-static and cyclic loads.

Ce travail de recherche propose un modèle numérique pour simuler les procédés par infusion de résine en utilisant la méthode des éléments finis. Ce modèle permet de représenter l'écoulement d'une résine liquide dans des préformes poreuses subissant de grandes déformations. Dans cette étude, une modélisation macroscopique est utilisée. Au niveau du procédé, une zone de résine liquide est déposée sur les préformes. Ces dernières étant considérées comme un milieu poreux. Les équations de Stokes et de Darcy sont utilisées pour modéliser l'écoulement de la résine respectivement dans le drainant et dans les préformes. L'originalité du modèle réside dans le fait qu'un seul maillage est utilisé pour les deux milieux. La discrétisation est réalisée avec des éléments mixtes : dans Stokes, des éléments P1+/P1 sont utilisés et dans Darcy, des éléments P1/P1 stabilisés avec une formulation multi-échelle sont employés. Des fonctions distances signées sont utilisées pour représenter l'interface entre Stokes-Darcy et pour représenter le front de résine. Concernant la déformation des préformes, une formulation Lagrangienne réactualisée est utilisée. Dans cette formulation Lagrangienne, le comportement des préformes humides est représenté à l'aide du modèle de Terzaghi dans lequel les préformes sèches ont un comportement élastique non-linéaire. La perméabilité est reliée à la porosité via la relation de Carman-Kozeny. Celle-ci est déterminée à partir de l'équation de conservation de la masse. Ce modèle a été implémenté dans ZéBuLoN. Plusieurs simulations numériques d'infusion de résine sont présentées à la fin de ce manuscrit. / This work proposes a numerical model to simulate the manufacturing processes by resin infusion using the finite element method. This model allows to represent the resin flow into porous preforms, which are themselves subject to large deformations. In this study, a macroscopic description is used. The preforms are considered as a porous medium. The Stokes and the Darcy equations are used respectively to describe the resin flow into the liquid zone and into the preforms.The originality of the model consists in using one single unstructured mesh. The discretization is ensured by using a mixed velocity-pressure formulation. Indeed, a P1/P1 formulation is employed throughout the entire discretized domain, stabilized in the Darcy region with a multi-scale formulation and in the Stokes subdomain with a hierarchical-based bubble, i.e. a P1+/P1 finite element. Signed distance functions are used both to represent the Stokes-Darcy interface and to capture the moving flow front. Concerning the deformations of the preforms, an updated Lagrangian scheme is used. In the Lagrangian formulation, the behavior of the wet preforms is represented by using the Terzaghi model in which the dry preforms have a non-linear elastic behavior. The permeability depends on the porosity through the Carman-Kozeny relationship. This model has been implemented in Zset. Several numerical simulations of manufacturing processes by resin infusion are presented at the end of this manuscript.

Le forgeage à chaud est un procédé de formage des métaux utilisé pour former des matériaux qui sont difficiles à former à froid ainsi que pour réaliser des géométries complexes. La réduction de la limite d’élasticité à haute température et une augmentation subséquente de l’aptitude à la mise en forme constituent le principal mécanisme à l’origine du procédé. Les méthodes numériques constituent un moyen efficace de prédire les états de contrainte / déformation du produit à différentes étapes de la mise en forme. Bien que les méthodes classiques soient suffisamment précises pour fournir une représentation appropriée du procédé, elles ont tendance à être coûteuses en ressources informatiques. Cela limite leur utilisation dans des cas concret, en particulier pour des études d’optimisation du procédé. L’approche pseudo inverse (API), développée dans le contexte du forgeage à froid 2D axisymétrique, fournit une estimation rapide des champs de contrainte et de contrainte dans le produit final pour une forme initiale donnée. Dans ce travail, l’API est étendue pour inclure les effets thermiques et visco-plastiques dans le procédé de forgeage ainsi que dans le cas général 3D. Les résultats sont comparés aux codes commerciaux disponibles basés sur les approches classiques pour montrer l’efficacité et les limites de l’API. Les résultats obtenus indiquent que l’API est un outil assez efficace pouvant être utilisé à la fois pour des simulations 2D et 3D du forgeage à chaud. / Hot forging is a metal forming process used to form difficult-to-form materials as well as to achieve complex geometries. The reduction of yield stress at high temperatures and a subsequent increase in formability is the primary mechanism that drives the process. Numerical methods provide an efficient means to predict the material yield and the stress/strain states of the product at different stages of forming. Although classical methods are accurate enough to provide a suitable representation of the process, they tend to be computationally expensive. This limits its use in practical cases especially for process optimization. Pseudo Inverse Approach (PIA) developed in the context of 2D axisymmetric cold forming, provides a quick estimate of the stress and strain fields in the final product for a given initial shape. In this work, the PIA is extended to include the thermal and viscoplastic effects on the forging process as well as to the general 3D case. The results are compared with commercially available software, based on the classical approaches, to show the efficiency and the limitations of PIA. The results obtained indicate that PIA is a quite effective tool that can be used for both 2D and 3D simulations of hot forging.

Une nouvelle approche appelée “Approche Pseudo Inverse” (API) est développée pour la modélisation du procédé de forgeage à froid des pièces axisymétriques. L'API est basée sur la connaissance de la forme de la pièce finale. Certaines configurations intermédiaires « réalistes » ont été introduites dans l'API pour considérer le chemin de déformations. Elles sont créées géométriquement sans traitement de contact et ensuite corrigées par la méthode de surface libre afin de respecter l'équilibre, les conditions aux limites et la condition d'incompressibilité. Un nouvel algorithme direct de plasticité est développé, conduisant à une méthode d'intégration plastique très rapide, précise et robuste même dans le cas de très grands incréments de déformations. Un modèle d'endommagement en déformation, est couplé à la plasticité et implémenté dans l'API. Les validations numériques montrent que l'API donne des résultats très proches des résultats de l'approche incrémentale mais en utilisant beaucoup moins de temps de calcul.L'API est adoptée comme solveur du forgeage pour la conception et l'optimisation des préformes du forgeage multi-passes. La rapidité et la robustesse de l'API rendent la procédure d'optimisation très performante. Une nouvelle technique est développée pour générer automatiquement le contour initial d'un outil de préforme pour la procédure d'optimisation. Les variables de conception sont les positions verticales des points de contrôle des courbes B-spline définissant les formes des outils de préforme. Notre optimisation multi-objectif consiste à minimiser la variation de la déformation plastique équivalente dans la pièce finale et la force du poinçon au cours du forgeage. Un algorithme génétique et un algorithme de recuit simulé sont utilisés pour trouver les points d'optimum de Pareto. Pour réduire le nombre de simulations de forgeage, un méta-modèle de substitution basé sur la méthode de krigeage est adopté pour établir une surface de réponse approximative. Les résultats obtenus par l'API en utilisant les outils de préforme optimaux issues de l'optimisation sont comparés à ceux obtenus par les approches incrémentales classiques pour montrer l'efficacité et les limites de l'API. La procédure d'optimisation combinée avec l'API peut être un outil numérique rapide et performant pour la conception et l'optimisation des outillages de préforme. / A new method called “Pseudo Inverse Approach” (PIA) is developed for the axi-symmetrical cold forging modelling. The PIA is based on the knowledge of the final part shape. Some « realistic » intermediate configurations are introduced in the PIA to consider the deformation path. They are created geometrically without contact treatment, and then corrected by using a free surface method in order to satisfy the equilibrium, the boundary conditions and the metal incompressibility. A new direct algorithm of plasticity is proposed, leading to a very fast, accurate and robust plastic integration method even in the case of very large strain increments. An isotropic damage model in deformation is coupled with the plasticity and implemented in the PIA. Numerical tests have shown that the Pseudo Inverse Approach gives very close results to those obtained by the incremental approach, but using much less calculation time.The PIA is adopted as forging solver for the design and optimization of preform tools in the multi-stage forging process. The rapidity and robustness of the PIA make the optimization procedure very powerful. A new method is developed to automatically generate the initial preform tool shape for the optimization procedure. The design variables are the vertical positions of the control points of B-spline curves describing the preform tool shape. Our multi-objective optimization is to minimize the equivalent plastic strain in the final part and the punch force during the forging process. The Genetic algorithm and Simulated Annealing algorithm are used to find optimal Pareto points. To reduce the number of forging simulations, a surrogate meta-model based on the kriging method is adopted to build an approximate response surface. The results obtained by the PIA using the optimal preform tools issued from the optimization procedure are compared to those obtained by using the classical incremental approaches to show the effectiveness and limitations of the PIA. The optimization procedure combined with the PIA can be a rapid and powerful tool for the design and optimization of the preform tools.