L'objectif de ce travail de thèse est de proposer un modèle numérique permettant de simuler la mise en forme de viroles métalliques par un procédé de formage incrémental. Ce procédé de mise en forme est particulièrement intéressant pour les industriels car il permet d'une part de réduire les pertes de matière première et d'autre part d'obtenir des pièces avec des propriétés mécaniques améliorées. La mise au point du procédé reste néanmoins une étape fastidieuse et le recours à l'outil numérique devient indispensable si l'on souhaite atteindre le meilleur rapport qualité/coût-délai. À cette fin, le comportement mécanique d'un acier faiblement allié a été étudié à travers divers essais de traction et de cisaillement. Les résultats de ces essais mettent en évidence que le matériau présente une faible anisotropie de comportement, une faible sensibilité à la vitesse de déformation et un écrouissage cinématique (effet Bauschinger). Cette base de données expérimentales a ensuite servi à l'identification de plusieurs modèles de comportement élastoplastique phénoménologique. Deux modèles numériques différents ont alors été développés au sein du code de calculs par éléments finis Abaqus afin de simuler le formage d'une pièce de référence et d'une étude de cas industriel. Les résultats de ces simulations montrent que le modèle tridimensionnel fournit une bonne représentation du procédé de formage a contrario du modèle bidimensionnel axisymétrique construit notamment afin de réduire les temps de calculs. Par ailleurs, une étude de la sensibilité du modèle éléments finis aux différents modèles de comportement identifiés a été menée. Les résultats de cette étude ont permis de mettre en évidence que la description du comportement plastique n'a qu'une faible influence sur la géométrie de la pièce déformée alors même que le comportement élastique joue un rôle bien plus important notamment au regard du calcul du retour élastique des pièces. Il a également été montré que le choix d'une forme particulière de géométrie de pièce permettait de réduire considérablement ce phénomène. Enfin, une étude sur la modélisation du comportement du matériau avec un modèle moins phénoménologique a été entreprise. Elle permet d'avoir un caractère plus générique avec un nombre de paramètres inférieurs aux modèles phénoménologiques en considérant les hétérogénéités du matériau. Elle ouvre la voie à l'utilisation de modèle de comportement prenant en compte des phénomènes physiques au sein d'un modèle éléments finis sur une opération de formage.

Afin de participer à l’émergence de technologies innovantes et contribuer aux objectifs de développement durable, un consortium composé de 13 partenaires (3 laboratoires académiques et 10 entreprises industrielles) a été formé. Le projet Thermofluid-RT consiste à développer un système de refroidissement à boucle de fluide diphasique à pompage mécanique. Le LMT-Cachan contribue à la conception de l'élément essentiel de la pompe, à savoir, le soufflet. Celui-ci est obtenu en soudant des feuilles d’acier inoxydable à durcissement structural très mince (70 µm) et doit fonctionner sans défaillance pendant 20 ans. Un dimensionnement fiabiliste du soufflet basé sur la théorie du maillon le plus faible est réalisé. Une méthode d'optimisation basée sur la technique d’identification intégrée à la corrélation d'images numériques aboutit à une géométrie d’éprouvette qui minimise l'incertitude des paramètres recherchés. Toutes les données brutes sont combinées à leur juste valeur grâce à une formulation Bayésienne basée sur l’hypothèse de bruits blanc gaussiens. La géométrie optimisée est testée sur Mini-Astrée, la nouvelle machine biaxiale du LMT. Plusieurs lois de comportements sont étudiées et testées sur le matériau étudié. De très fines feuilles du même acier, sont de même testées à l’aide d’essais uni-axiaux et multiaxiaux. Une analyse microscopique est conduite au travers d’un micro-essai de traction sur 2 grains et l’identification de paramètres d’une loi de plasticité cristalline. Enfin, le dimensionnement probabiliste du soufflet est validé à partir de plusieurs mesures expérimentales infra-rouges sur un nouveau banc d’essai. Un modèle probabiliste à deux échelles permet la caractérisation des phénomènes observés relatifs au matériau d’étude et au composant lui-même. / In the context of leveraging and accelerating innovative technological solutions that contribute to meeting sustainable goals, a consortium composed of 13 partners (3 academic laboratories and 10 industrial companies) has been formed. The THERMOFLUID-RT project consists in developing a two-phase fluid loop driven by a mechanical pump. LMT-Cachan contributes to the design of the critical component of the pump, i.e., bellows. This component obtained by welding very thin (70 µm) pre-deformed sheets made of precipitation-hardened stainless steel is expected to operate without failure for 20 years. First, a new probabilistic design methodology allows for the fatigue design of the component. Second, an optimization methodology based on full field measurements and Integrated Digital Image Correlation (IDIC) allows the sample geometry to be designed with the least uncertainty of the sought parameters. All data account for the latter thanks to a Bayesian foundation that equitably weights all measurements. Third, the optimized cruciform geometry has been tested in a new compact biaxial machine, mini-Astrée that allows for a fast, yet robust identification. Several elasto-plastic models with increasing complexity are investigated to probe the material behavior. Ultra-thin sheets are also tested in uniaxial and biaxial experiments thanks to an anti-wrinkling setup. Fourth, a microscopic analysis is also performed via quasi-3D IDIC and a uniaxial micro-specimen allow a crystal plasticity law to be characterized. Last, the probabilistic fatigue design of the bellows based on the weakest link theory and a two-scale probabilistic model is validated with infrared measurements in the high cycle fatigue regime.

Le travail de thèse porte sur la modélisation de la mise en forme et de la tenue mécanique des assemblages par déformation plastique avec comme principale application le rivetage auto-poinçonneur. Dans un premier temps, le développement d'un modèle de comportement, capable de prendre en compte l'ensemble des phénomènes physiques apparaissant lors de la pose du rivet ou de la tenue mécanique du point d'assemblage, est abordé. Le comportement retenu est un comportement élasto-plastique endommageable. Le modèle d'endommagement, basé sur le modèle de Lemaitre et le principe d'équivalence en déformation, a été amélioré afin de prendre en compte l'effet de fermeture des fissures ainsi qu'une limite de triaxialité en dessous de laquelle l'endommagement ne peut plus évoluer. Ce modèle de comportement a alors été implémenté dans la suite logicielle Forge2005®. Afin de s'affranchir du phénomène de localisation, des méthodes d'endommagement non locales ont été introduites. Dans un deuxième temps, la caractérisation des paramètres des lois d'écrouissages et d'endommagement est réalisée en utilisant une méthode d'indentification par analyse inverse basée sur un algorithme évolutionnaire couplé à un méta-modèle. L'ensemble des essais mis en place pour alimenter cet algorithme est abordé de façon précise. Finalement, le modèle numérique développé est utilisé pour simuler le rivetage auto-poinçonneur et la tenue mécanique du point d'assemblage. Pour ce dernier point, un montage ARCAN a été conçu pour solliciter les éprouvettes assemblées en cisaillement, traction ou encore pour des sollicitations mixtes. Les résultats obtenus sont alors comparés avec l'expérience et montrent une bonne corrélation. Le modèle numérique ainsi développé peut également être utilisé pour la majorité des technologies d'assemblage par déformation plastique, mais il constitue également un outil d'aide à l'optimisation de la tenue mécanique, comme cela est illustré en dernière partie de ce travail.

Le déminage mécanique consiste à scarifier le sol avec une " charrue ", outil formé d'une lame en V munie de plusieurs dents, poussée par un engin : les dents déstructurent le sol et le font remonter devant la lame qui l'évacue sur le côté avec les mines qu'il contient. L'objectif de nos travaux est de mettre en œuvre la simulation numérique par éléments finis du problème fortement non linéaire issu de la modélisation de la scarification du sol. Le code d'éléments finis implicite Forge3®, dédié à la mise en forme des métaux, a été choisi comme support numérique. Il permet de prendre en compte les grandes déformations, en particulier à l'aide de son remaillage automatique. Nous avons dans un premier temps implanté dans ce code deux modèles élastoplastiques de comportement, l'un incompressible réservé aux sols fins saturés, purement cohésifs, l'autre compressible, fondé sur la notion d'état critique, pour les matériaux purement frottants ou cohésifs et frottants. Ces modèles adoucissants sont intégrés par un schéma de retour radial généralisé, au sein d'une formulation implicite du problème aux limites. Nous montrons que la matrice de raideur est non symétrique dans le cas compressible et que la symétrisation du système ne conduit pas à une approche robuste. Nous avons donc importé et testé un solveur itératif non-symétrique : Bi-CGSTAB. Nous avons validé la programmation de ces modèles sur des essais triaxiaux. Pour les comportements adoucissants, on constate des oscillations dans la relation contrainte/déformation passé le pic de contrainte. Ces difficultés numériques sont traitées par linéarisation et régularisation. Dans un deuxième temps, nous avons mis en œuvre des simulations de scarification pour différents niveaux de complexité : une dent seule, une dent + une tranche de lame, plusieurs dents, en faisant suffisamment avancer l'outil pour atteindre le régime stationnaire ; cela se révèle d'autant plus long que l'outil est large, ce qui nous amène au calcul très intensif. L'influence de paramètres géométriques comme l'inclinaison de la dent ou l'angle d'étrave du système a été mise en évidence, elle est qualitativement conforme aux observations expérimentales. L'étude de l'influence des paramètres des modèles de comportement montre la prépondérance de ceux liés à l'état critique, i.e. aux propriétés mécaniques après de grandes déformations. Enfin nous avons validé qualitativement le code en termes de modes d'écoulements et de répartition des efforts pour des outils multi-dents. La comparaison quantitative des efforts reste à affiner en revenant sur le comportement choisi et son implémentation.

Grâce à leurs bonnes propriétés mécaniques spécifiques, les matériaux cellulaires architecturés présentent un fort intérêt pour répondre aux problématiques du secteur aéronautique. Cependant, la modélisation d'une structure macroscopique incluant un matériau cellulaire nécessite, soit de modéliser complètement l'architecture à l'échelle mésoscopique - ce qui est coûteux en temps de calcul - soit d'utiliser un Milieu Homogène Equivalent (MHE). Ainsi, cette thèse propose de caractériser un matériau cellulaire modèle constitué d'un empilement de tubes, selon un motif carré ou hexagonal, puis d'identifier un modèle phénoménologique rendant compte du comportement mécanique inélastique du matériau. Dans un premier temps, le matériau est caractérisé sous chargements multi-axiaux à l'aide de simulations éléments finis périodiques en petites déformations. Le comportement homogénéisé en petites déformations est ensuite utilisé pour l'identification d'une Loi Homogène Equivalente (LHE) compressible et anisotrope, qui permet la modélisation de structures sandwichs en remplaçant le coeur cellulaire par son MHE. Une comparaison est réalisée entre les réponses mécaniques des simulations de référence complètement maillées et celles utilisant l'approche par MHE, validant ainsi la pertinence de la méthode multi-échelle de modélisation proposée. La caractérisation en grandes déformations des deux types d'empilement est ensuite menée. D'abord, les effets de bords et les instabilités qui gouvernent le comportement macroscopique sont étudiés. Puis, après une étude du volume élémentaire représentatif des empilements, la caractérisation du comportement inélastique par la technique de l'homogénéisation périodique est réalisée. Le comportement adoucissant en compression de l'empilement hexagonal est ainsi étudié. Finalement, une extension des LHE identifiées en petites déformations est proposée pour rendre compte du comportement en compression du matériau observé en grandes déformations. / Cellular materials have excellent specific properties, which make them attractive for aeronautical applications. However, modelling macroscopic structures including a cellular material is either very costly in terms of computational time if the whole mesoscopic structure is considered or a Homogeneous Equivalent Medium (HEM) has to be used. This Ph.D. dissertation presents, the characterisation of a cellular material built from a stacking of tubes with a square or hexagonal based pattern and the identification of a phenomenological model of their inelastic mechanical behaviour. First, the material is characterised for multi-axial loadings through a periodic finite element model in small deformations for each tube stacking pattern. The macroscopic behaviour is then used to identify a compressible anisotropic Homogeneous Equivalent Law (HEL). Within the infinitesimal strain hypothesis, a comparison is carried out between reference full scale models and HEM based ones of sandwich structures with a cellular core, confirming the relevance of the proposed multi-scale method. Then, the mechanical behaviour of each tube stacking is characterised for large deformations in order to study the influence of the boundary size effects and the instabilities in the core on the macroscopic behaviour of sandwich structures. After a study on the representative volume element, the macroscopic inelastic behaviour is characterised through the periodic homogenisation technique, especially the softening observed in compression for the hexagonal pattern. Finally, an extension of the HELs identified in small deformations is proposed to model the behaviour observed in large deformations.

Les paramètres de loi de comportement s'identifient habituellement à partir d'éprouvettes standardisées présentant des champs de déformations homogènes. L'essor actuel des méthodes inverses incite à l'utilisation d'essais plus hétérogènes. L'objectif de cette thèse est de développer une procédure d'identification inverse en utilisant une méthode de recalage de modèles éléments finis basée sur des expériences à hétérogénéités croissantes des champs de déformation. À cette fin, des essais plans avec différentes géométries d'éprouvettes, en titane commercialement pur, ont été effectués et un essai original hors plan fortement hétérogène a été mis au point. Dans ce dernier cas, un essai unique réalisé à l'aide d'une machine de traction uni-axiale permet d'obtenir des sollicitations d'expansion, mais aussi de traction et de cisaillement selon deux directions orthogonales. Des mesures de champs cinématiques planes et non planes (stéréo-corrélation d'images) sont utilisées afin d'évaluer les déformations jusqu'à striction. À partir de ces mesures en grandes déformations, on propose de mener une comparaison entre les informations issues de l'expérience et un modèle dupliquant les conditions expérimentales (conditions aux limites issues de la corrélation d'images). Cette comparaison conduit, via un processus d'optimisation, à l'actualisation itérative du jeu de paramètres de la loi de comportement utilisée dans le modèle éléments finis. Le nombre de paramètres ainsi identifiés varie selon le modèle de comportement : élastique anisotrope puis élasto-plastique anisotrope. Différentes validations ont mis en évidence l'intérêt de considérer des champs de déformation fortement hétérogènes. Dans ce cadre, des prédictions numériques d'opérations d'emboutissage, réalisées à l'aide des différents jeux de paramètres identifiés sont comparées à l'expérience. On montre que la prise en compte d'un modèle élastique anisotrope plutôt qu'isotrope fournit un état de contrainte plus réaliste en début de plasticité, améliorant considérablement l'identification du modèle plastique. Enfin, une approche couplée entre mesures de champs cinématiques et thermiques (thermographie infrarouge) a été développée. Les paramètres d'un modèle mécanique sont identifiés puis utilisés pour l'évaluation d'un modèle thermique calculant les sources de chaleur. Par cette méthode, plusieurs modèles d'évolution de la fraction d'énergie dissipée sous forme de chaleur ont été identifiés (coefficient de Taylor et Quinney).

Les catastrophes naturelles sont de plus en plus fréquentes et sont à l'origine de conséquences socio-économiques souvent graves et coûteuses. La prévision de ces risques, notamment en ce qui concerne les glissements de terrain, est dès lors indispensable. Dans cette problématique, la modélisation numérique doit être pertinente et prédictive. Actuellement, les glissements sous une faible pente (inférieure à 14°) voire très faible (inférieure à 8°), où la rupture est souvent de type diffuse, ne peuvent être expliqués ni par les méthodes empiriques classiques ni par l'analyse limite traditionnelle (théorie de plasticité). De nouvelles approches sont donc nécessaires pour rendre compte de ces mécanismes de glissement. La condition suffisante de stabilité de Hill (1958) basée sur le signe du travail du second ordre pourrait être un critère pertinent pour analyser ces modes de rupture diffus.<br /><br />Ce travail de thèse porte d'une part sur l'analyse et l'implémentation du critère de stabilité de Hill dans deux codes de calculs par éléments finis (LAGAMINE et PLAXIS) et d'autre part sur leur application, dans le cadre des modélisations numériques de problèmes aux limites non linéaires, aux glissements de terrain.<br /><br />L'étude du travail du second ordre, menée sur des modèles de comportement incrémentaux non linéaires, montre que certains sols exhibent de larges domaines potentiellement instables à l'intérieur strict du critère limite de plasticité. Ces analyses sont corroborées par les résultats expérimentaux.<br /><br />Nous montrons, dans ces travaux, par des approches analytique et numérique dans le cadre de la théorie de l'élasto-plasticité, que le critère de Hill est toujours vérifié avant la condition limite de plasticité (critère de Mohr-Coulomb) et le critère de localisation de déformation (critère de Rice).<br /><br />L'utilisation du critère de Hill nous permet d'analyser le glissement de la côtière de Trévoux après une période de fortes précipitations en 1983. Ces travaux de modélisations ont été menés dans le cadre d'un couplage hydromécanique en milieu non saturé.<br /><br />La même démarche a été menée pour l'analyse du glissement catastrophique du versant de Las Colinas (El Salvador) provoqué par un séisme en 2001. La définition du chargement considéré est basée sur une méthode de type pseudo-statique.

La plasticite non-associee est un facteur important si ce n'est primordial dans le processus de localisation de la deformation cassante. Un code numerique particulierement adapte a reproduire les comportements mecaniques tres varies des roches nous permet de modeliser efficacement l'initiation et la propagation des deformations localisees a travers la lithosphere. La region de l'Ocean Indien est un cas spectaculaire de flambage lithospherique , accompagne d'une activite sismique intense. Plusieurs modeles analytiques, analogiques et numeriques rendent compte du processus de ambage. On developpe un modele numerique qui pour la premiere fois reproduit le developpement de ambage simultan'avec celui de failles. En utilisant une rheologie elasto-visco-plastique (critere de Mohr-Coulomb et loi de puissance dependant du geotherme), la fracturation di use a ecte d'abord l'ensemble de la lithosphere competente, est suivie de ambage et simultanement de failles qui se localisent progressivement aux points d'inflexion des plis. Les modeles numeriques testes dans ce memoire renseignent sur le comportement rheologique de la lithosphere, en particulier sur l'importance de l'angle de friction interne, qui permet de supporter les de ections liees au ambage. Les contraintes reelles sur les amplitudes montrent une dependance combinee du taux de raccourcissement, de la friction, et de l'epaisseur de la couche inferieure ductile. Des domaines intracontinentaux tels que l'Ouest du Gobi et le centre de l'Australie presentent les caracteristiques d'un ambage d'echelle lithospherique, malgre les heterogeneites pre-existantes. Une modelisation numerique confirme le role de l'age thermotectonique et de la rheologie crustale dans le developpement d'un mode de flambage. La dimension caracteristique des couches qui ambent determine l'evolution du modele vers un epaississement localise et une orogenese. Dans ce cas la lithosphere s'epaissit sous de larges bassins intracontinentaux de tailles comparables a la realite. L'insertion d'un modele d'erosion confirme son role stabilisant dans l'evolution de la topographie, alors qu'en profondeur le ambage est ampli e, et la localisation facilitee. L'elevation du relief obtenu varie de 0.1 a 1mm/an.

Cette thèse concerne la modélisation du comportement mécanique de poudres ductiles au cours de leur mise en forme par le procédé de compression à froid. L’approche utilisée consiste à modéliser la poudre par un échantillon numérique de 50 particules sphériques distribuées aléatoirement. Chaque particule est maillée en éléments finis avec un comportement élasto-plastique et interagit par contact avec ses voisines. Sous l’effet du chargement, les particules se déforment et subissent de grandes transformations.La méthode développée consiste à imposer aux frontières de l'échantillon numérique des conditions aux limites en déplacement ou en force afin de simuler une réponse continue équivalente en termes de contrainte-déformation pour plusieurs chemins de chargement depuis les faibles valeurs de densité relative jusqu'aux fortes valeurs. Outre le fait que la méthode permet d’appliquer par la simulation des chemins de chargement inaccessibles aux techniques expérimentales, elle fournit de nombreuses indications sur la phénoménologie du comportement aux frontières de l'échantillon en relation avec la morphologie des particules déformées plastiquement et l’évolution des surfaces de contact inter-particules.Cette réponse obtenue sera étudiée au moyen des notions de l'élasto-plasticité classique, à savoir la surface de charge et la loi d'écoulement. L'accent est mis sur l’analyse de l’écoulement plastique et l'existence du potentiel plastique. Les résultats montrent que la direction d’écoulement est approximativement unique pour les états de contrainte éloignés du point de chargement. En revanche, la direction d’écoulement cesse d’être unique à proximité de point de chargement, ce qui révèle le caractère anisotrope de l’écrouissage. Ce comportement, propre aux poudres ductiles, est intimement lié à la déformation des grains. La dégradation anisotrope des propriétés mécaniques résultant du mécanisme de dilatance est également étudiée. / This thesis is concerned with the modelling of the mechanical behaviour of ductile powders during cold compaction process. A numerical method was implemented, in which a powder sample was assimilated to an assembly of 50 meshed particles in a cubic box which was submitted to compressive loadings simulated by means of the finite element method. Each particle was meshed and assigned an elastic-plastic constitutive behaviour. Particles interacted through mechanical, frictional contact and underwent large transformation during loading.The method developped in the thesis consisted in imposing boundary conditions as displacements or forces on the bounding walls of the numerical sample. The mechanical response of the sample was averaged to obtain the equivalent response of a continuum in terms of stress and strain. Various loading paths, including loading paths that are technically not attainable by experimental means, were applied up to different values of relative densities, from low (60%) to high values (98%). Much information on the phenomenological behaviour of the numerical sample could be obtained, in relation with the morphology of deformed particles and the evolution of contact surfaces.The results were interpreted based on the concepts of classical elasto-plasticity, i.e., yield surface and flow rule. The main focus of the study was the analysis of incremental plastic flow and the existence of a plastic potential. Results showed that the flow rule postulate, based in the plastic potential, could be considered as valid for stress states relatively far from the loading point. In the vicinity of the loading point, the direction of the plastic strain increment vector ceases to be unique. This behaviour, which is an original feature of ductile powders mechanical behaviour, was attributed to the anisotropic strain-hardening processes at stake, related to the anisotropic formation of contact surfaces between particles. The drop in mechanical properties related to dilatancy was also studied.

L'objectif de ce travail de thèse est de proposer un modèle numérique permettant de simuler la mise en forme de viroles métalliques par un procédé de formage incrémental. Ce procédé de mise en forme est particulièrement intéressant pour les industriels car il permet d'une part de réduire les pertes de matière première et d'autre part d'obtenir des pièces avec des propriétés mécaniques améliorées. La mise au point du procédé reste néanmoins une étape fastidieuse et le recours à l'outil numérique devient indispensable si l'on souhaite atteindre le meilleur rapport qualité/coût-délai. À cette fin, le comportement mécanique d'un acier faiblement allié a été étudié à travers divers essais de traction et de cisaillement. Les résultats de ces essais mettent en évidence que le matériau présente une faible anisotropie de comportement, une faible sensibilité à la vitesse de déformation et un écrouissage cinématique (effet Bauschinger). Cette base de données expérimentales a ensuite servi à l'identification de plusieurs modèles de comportement élastoplastique phénoménologique. Deux modèles numériques différents ont alors été développés au sein du code de calculs par éléments finis Abaqus afin de simuler le formage d'une pièce de référence et d'une étude de cas industriel. Les résultats de ces simulations montrent que le modèle tridimensionnel fournit une bonne représentation du procédé de formage a contrario du modèle bidimensionnel axisymétrique construit notamment afin de réduire les temps de calculs. Par ailleurs, une étude de la sensibilité du modèle éléments finis aux différents modèles de comportement identifiés a été menée. Les résultats de cette étude ont permis de mettre en évidence que la description du comportement plastique n'a qu'une faible influence sur la géométrie de la pièce déformée alors même que le comportement élastique joue un rôle bien plus important notamment au regard du calcul du retour élastique des pièces. Il a également été montré que le choix d'une forme particulière de géométrie de pièce permettait de réduire considérablement ce phénomène. Enfin, une étude sur la modélisation du comportement du matériau avec un modèle moins phénoménologique a été entreprise. Elle permet d'avoir un caractère plus générique avec un nombre de paramètres inférieurs aux modèles phénoménologiques en considérant les hétérogénéités du matériau. Elle ouvre la voie à l'utilisation de modèle de comportement prenant en compte des phénomènes physiques au sein d'un modèle éléments finis sur une opération de formage. / This Ph.D. thesis aims to propose a numerical model for modelling the incremental forming process of a metallic thin wall short tube. This process is especially interresting for manufacturer because it allows to reduce the material raw losses and to increase the mechanical properties of the part. The technical development of the process is a tedious step and the numerical simulation becomes necessary if the best ratio quality/cost-delay must be reach. To achieve this goal, the mechanical behaviour of a low-allow steel was studied with uniaxial and shear tests. The results of this tests highlight that the material has a low anisotropic behaviour, a low strain rate sensitivity and kinematic hardening (Bauschinger effect). Then, this experimental database was used to identify several elastic-plastic phenomenological behaviour laws. Two different numerical models were developed in the finite element code Abaqus to simulate the forming of a reference part and an industrial study. Results of simulation show that tridimensional model is a good representation of the forming process contrary to the bidimensional model built to reduce the CPU time. Furthermore, a finite element model sensibility study of the identifying behaviour models has been carried out. The results of this study highlight that plastic behaviour description has a low influence on the formed part geometry. However, the elastic behaviour has an important influence in particular on the calculated springback. Moreover, it was shown that the choice of a particular part geometry reduce in an important way this phenomenon. Finaly, a study about the material behaviour modeling with a less phenomenological model was done. This latter is a more generic model with a lower parameters number compared with phenomenological models taking into account material heterogeneities. This study innovates for the behaviour model taking account of physical phenomenon in finite elements model on a forming process.