Le traitement thermique des aciers après chauffage par induction localisé est un procédé aujourd'hui répandu, notamment dans l'industrie automobile. Ce type de procédés a fait ses preuves dans le cadre de pièces massives, mais les distorsions générées sont encore parfois insuffisamment maîtrisées lorsqu'il s'agit de tôles minces. PSA-Peugeot-Citroën s'intéresse alors à l'étude et la compréhension des phénomènes en jeu ainsi qu'à la simulation numérique de tels procédés. C'est alors dans cette démarche et dans le cadre d'un partenariat entre le Cemef et PSA, que s'inscrit ce travail appliqué au traitement thermique d'un renfort de pied central. Un procédé aussi complexe nécessite la compréhension de la thermique, de la mécanique, de la métallurgie, de l'électromagnétisme, ainsi que de leurs interactions mutuelles au chauffage comme au refroidissement. Ce travail se trouve alors au carrefour de plusieurs disciplines comme la thermomécanique et l'électromagnétisme, ainsi que les méthodes numériques et l'étude expérimentale. Il fait suite à divers travaux réalisés au laboratoire concernant la thermique ou la mécanique, ainsi qu'une thèse précédente portant sur la modélisation numérique couplée de la trempe. Il constitue également la première approche 3D des procédés de chauffage par induction et des couplages multi physiques qui en découlent. La présentation de ce travail se décompose en 3 grandes parties. Premièrement, on décrit le contexte industriel, la problématique et l'historique de l'étude du procédé, pour en arriver aux études expérimentales réalisées dans ce travail : une porte sur le procédé industriel et l'autre sur un modèle simplifié que nous avons conçu pour une meilleure compréhension des phénomènes physiques qui génèrent les distorsions. Ces études expérimentales, riches en résultats, soulignent la nécessité d'un outil numérique pour aller encore plus loin dans la compréhension physique du procédé. Nous poursuivons alors sur trois chapitres présentant les modèles numériques par éléments finis, développés et utilisés dans le code de calcul pour simuler les problèmes couplés : thermique / électromagnétisme, mécanique / métallurgie et l'ensemble de la structure couplée. La troisième et dernière partie du rapport traite, d'une part, de la mise en donnée d'un cas semi-industriel ainsi que de la caractérisation des paramètres, et d'autre part, des résultats numériques obtenus et de leur comparaison avec les résultats expérimentaux. Le modèle développé permet d'estimer dans une première approche les distorsions subies par une structure mince lors de son traitement thermique localisé par induction. Le code de calcul actuel forme une base solide pour de futures évolutions permettant de simuler des problèmes industriels complexes.

Les travaux présentés dans cette thèse s'inscrivent dans le cadre de l'analyse et la modélisation de la distribution du champ magnétique dans les machines électriques à circuit magnétique feuilleté avec des tôles à grains orientés. L'anisotropie des tôles magnétiques des transformateurs induit des phénomènes 3D complexes au niveau des joints magnétiques où les tôles se chevauchent. Concernant l'augmentation de l'efficacité énergétique des machines tournantes, de nouvelles structures à base de tôles à grains orientés font leur apparition. Néanmoins, simuler finement en 3D un paquet de tôles minces recouvertes d'un isolant de quelques microns d'épaisseur conduit à des temps de calculs très importants.Dans cette optique, nous présentons une méthode d'homogénéisation dont le but est de définir les caractéristiques magnétiques équivalentes d'un empilement quelconque de tôles et d'entrefers. Sa formulation est basée sur la minimisation de l'énergie et la conservation du flux magnétique. Les résultats de cette méthode appliquée à un joint magnétique de transformateur à joints step-lap sont confrontés à des mesures expérimentales et à un modèle élément fini 3D. Ce dernier requiert de connaître les caractéristiques magnétiques des tôles dans les directions de laminage, transverse et normale. La détermination de la perméabilité des tôles dans la direction normale est problématique et constitue un point original de notre étude. Deux méthodes, analytique et numérique, s'appuient sur les mesuresd'un banc de caractérisation en régime statique pour déterminer cette perméabilité normale. / This thesis focuses on the analysis and the modeling the magnetic flux distribution in electrical machines with anisotropic laminated magnetic circuit. The anisotropy of magnetic sheets in transformers induces complex 3D phenomena in step-lap magnetic joints where the sheets are overlapped. Moreover, in order to increase the energy efficiency of rotating machines, new structures based on grain-oriented electrical steel are developed.However, an accurate 3D simulation of a laminated core with thin sheets and insulation of a few microns leads to very large computation time. In this context, we present a homogenization method, which purpose is to define the equivalent magnetic characteristics of any laminated core made of sheets and air gaps. Its formulation is based on the energy minimization and the magnetic flux conservation. The results of this method applied to a step-lap magnetic joint are compared with experimental measurements and a 3D finite element model. The latter requires to know the magnetic characteristics of the sheets in the rolling, transverse and normal directions. The determination of the sheets permeability in the normal direction is problematic and it constitutes an original point of our study.Two methods, analytical and numerical, based on measurements obtained with a static characterization bench makes possible the determination of the normale permeability.

Cette thèse met l'accent sur la recherche scientifique pour développer une classe de procédés à hautesvitesses de déformation des tôles métalliques en alliages d'aluminium et en acier à haute résistance (AHR).Ces technologies emploient une impulsion de pression de courte durée qui propulse la tôle dans une matrice.Ces procédés sont généralement décrits comme procédés de formage par impulsion de pression (PPF). Letravail proposé dans ce mémoire de thèse a permis de surmonter trois obstacles pour l'utilisation desprocédés PPF et la fabrication à moindre coût de structures légères. Le premier obstacle a été le manque decorrélation entre formabilité et vitesses de déformation qui se développent lors d’un procédé PPF. Nous avonsproposé d’analyser la formabilité et la rupture des tôles, et de caractériser les vitesses de déformation et leurl'hétérogénéité pendant le procédé PPF. Le deuxième obstacle a été le manque de lois constitutives validéespour les métaux déformés par le procédé PPF. Nous avons étudié la microstructure et l'évolution despropriétés mécaniques durant le procédé PPF. Le troisième obstacle est le manque de modèles de formabilitéprédictifs validés pour le procédé PPF. Nous avons utilisé la méthode Marciniak-Kuczynski pour la prédictionde la formabilité de l’alliage AA5182 et de l’acier DP600 sous un large éventail de vitesses de déformation etsous différentes directions de ces vitesses. La combinaison de ces résultats de recherche permet une plusgrande capacité prédictive pour concevoir et développer des procédés PPF pour composants d’automobiledésirés à partir d'aluminium et d’acier AHR. / This thesis focuses on scientific investigation to develop and enable a class of high strain ratesheet metal forming of aluminum alloys and advanced high strength steel (AHSS). These technologiesemploy a short duration pressure‐pulse to drive sheet metal into single‐sided dies, and can generally bedescribed as pulse pressure forming (PPF) processes. The work under this thesis has overcome threetechnical barriers to using PPF processing for more cost effective lightweight vehicles. The first technicalbarrier was the lack of understanding of the interrelationship between formability and measured strainrates that develop during PPF processing. The work under this thesis investigated the formability andfracture of sheet metals during PPF, and characterized the strain rate and the strain rate heterogeneity.The second technical barrier was the lack of a validated constitutive model for lightweight materialsduring PPF processing. The work under this thesis investigated the microstructure and mechanicalproperty evolution in metals during PPF. The third technical barrier was the lack of validated andpredictive formability models for PPF processes. The work under this thesis used the Marciniak andKuczynski method of formability prediction to predict the formability of both aluminum alloy AA5182and AHSS alloy DP600 across a wide range of strain rates and strain rate directions. The combination ofthese research results permits a more predictive capacity to design and develop PPF manufacturingprocesses for a desired automotive component made from aluminum alloys and AHSS.

Le travail de ma thèse s'inscrit dans le cadre d'un projet industriel proposé par une entreprise luxembourgeoise (Gottschol-Alcuilux) et en collaboration avec le Centre de Recherche Public Henry Tudor du Luxembourg (Laboratoire des Technologies Industriels (LTI)). L'objectif consiste à mettre au point une méthode numérique de détermination de la forme des outils d'emboutissage et du flan de pièces de précision en tôlerie fine pour que ce dernier, une fois déformé en une ou plusieurs opérations, correspond à la définition tridimensionnelle du cahier des charges. La méthode a pour objectif de remplacer une démarche expérimentale coûteuse par essais–erreur. <br /><br />Deux démarches numériques ont été développées, la première est relative à la détermination de la forme du flan. Elle consiste à estimer la forme du flan par Approche Inverse en partant de la forme 3D demandée. Puis, un logiciel de simulation incrémental par éléments finis en 3D est utilisé dans une procédure d'optimisation heuristique pour déterminer la forme du flan. Dans la deuxième démarche, il s'agit de déterminer la forme des outils d'emboutissage en utilisant le logiciel de simulation incrémental couplé avec une méthode de compensation du retour élastique en 2D. La démarche numérique est validée expérimentalement dans le cas d'un emboutissage réalisé en une ou plusieurs passes, à l'aide d'une presse manuelle, sans serre flan et avec des outils de forme très complexe.

Le travail présenté dans ce mémoire porte sur la conception d'un dépileur électromagnétique de tôles amagnétiques sur les sites d'emboutissage des usines de construction automobile. Après une étude de caractérisation du comportement mécanique du paquet de tôles huilées et soumis à un effort de pelage, deux principes de séparateur électromagnétique ont été exploités. Le premier consiste à faire interagir, dans la tôle supérieure, un courant induit et un champ alternatif produits par une structure adaptée et placée en bordure du paquet de tôles. Plusieurs structures à champ alternatif ont été proposées et modélisées par la méthode des éléments finis et leurs performances, pour l'effort de pelage, analysées et commentées. Le deuxième principe repose sur l'interaction d'un champ magnétique produit par une structure à aimants permanents et un courant électrique circulant en bordure de la tôle supérieure et délivré par des électrodes. Une méthode de calcul tridimensionnelle semianalytique du champ magnétique créé dans le plan de la tôle supérieure est proposée ainsi qu'un modèle analytique de répartition de la densité de courant dans cette même tôle issu de la résolution de l'équation de Laplace par la méthode de séparation des variables. Cette modélisation a permis l'optimisation de l'effort de pelage et un prototype a été construit en laboratoire. Les résultats obtenus ont donné satisfaction.

L'actuel investigation numérique du Gurson-Tvergaard-Needleman (GTN) modèle avancé est une extension du travail de Ben Bettaieb et al. (2011). Le modèle a été implémenté à l'aide d'une sous routine (VUMAT) contenu dans le code commerciale d'éléments finis Abaqus/explicit. Le modèle d'endommagement améliore l'original en intégrant les trois mécanismes d'endommagement, la nucléation, la croissance, et la coalescence des cavités. Le modèle d'endommagement intègre les lois de nucléation et de croissance basés sur les phénomènes purement physiques. Ces nouvelles contributions incluant l'influence de l'écrouissage cinématique, ont été validées par les résultats de mesures expérimentales de tomographie à rayon X à haute résolution. Aussi, l'implémentation numérique de l'écrouissage cinématique dans le modèle modifié a contraint de proposer et de réarranger la définition de la triaxialité que l'on trouve habituellement dans la littérature. A coté de cela, un second critère d'initiation à la rupture basé sur l'ultime distance inter-cavités a été inclue afin de localiser et de quantifier avec plus de précision la distribution des déformations peu avant que le matériau ne casse complètement. L'actuel modèle d'endommagement a été appliqué dans des conditions industrielles pour prédire l'évolution de l'endommagement, l'état de contraintes, et l'initiation à la rupture pour différentes géométries de tôles et sur des essais d'emboutissage de tôles minces. / This numerical investigation of an advanced Gurson–Tvergaard–Needleman (GTN) model is an extension of the original work of Ben Bettaieb et al. (2011). The model has been implemented as a user-defined material model subroutine (VUMAT) in the Abaqus/explicit FE code. The current damage model extends the previous version by integrating the three damage mechanisms: nucleation, growth and coalescence of voids. Physically based void nucleation and growth laws are considered, including an effect of the kinematic hardening. These new contributions are based and validated on experimental results provided by high-resolution X-ray absorption tomography measurements. Also, the numerical implementation of the kinematic hardening in this damage extension has obliged to readapt the classical triaxiality definition. Besides, a secondary fracture initiation criterion based on the ultimate average inter-cavities distance has been integrated to localize and quantify with good accuracy the strain distribution just before the material fails apart. The current damage model is applied in industrial conditions to predict the damage evolution, the stress state and the fracture initiation in various tensile thin flat sheet geometries and the cross-die drawing tests.

L'objectif de ces travaux de thèse est de proposer des approches de modélisation et d'expérimentation de l'impact de structures déformables et indéformables sur différents milieux. Différents modèles analytiques et des simulations numériques sont développés en comparaison aux résultats expérimentaux. Une première partie se consacre à la caractérisation de la similitude entre la réponse à l'impact à l'eau d'un solide et la réponse d'un solide impactant une structure déformable. Des simulations éléments finis (EF) et SPH sont réalisées pour l'impact à l'eau d'un tube cylindrique (sans rupture). Un modèle analytique d'impact à l'eau est proposé pour prédire l'évolution de l'effort (pic, durée). L'analyse des résultats permet de dimensionner un programmateur d'impact solide reproduisant le pic d'effort. Des simulations EF de l'impact sur un tube cylindrique, à géométrie adapté, dans la direction longitudinale, sont réalisées et comparées à quelques expériences tests. Le «flambage dynamique» (dû au comportement inélastique du matériau et aux ondes de déformations) des tubes est alors observé. Une seconde partie traite du cas de la perforation sous impact d'une tôle mince à faibles vitesses d'impact (< 10 m/s, vitesse de déformation < 1000 s-1). Des essais sur puits de chute instrumenté (force, déplacement, déformée de tôle, avancée de fissure) sont analysés. Des simulations EF en éléments coques avec un critère de rupture ductile par endommagement sont réalisées. Les paramètres de rupture dynamique sont identifiés par méthode inverse à l'aide d'essais de résilience Charpy sur l'alliage d'aluminium de désignation 2024 T3. Une analyse des pics de force lors de l'impact permet une meilleure compréhension des mécanismes de perforation. En parallèle, un nouveau modèle analytique, basé sur les énergies impliquées lors de l'impact, est proposé et comparé aux simulations EF. L'étude numérique de la perforation est étendue aux grandes vitesses d'impact et de déformation (100 - 1000 m/s, vitesse de déformation <100 000 s-1) pour identifier les transitions des différents mécanismes de perforation connus (pétalisation, fragmentation des pétales, fragmentation complète). / The objective of this work is to propose approaches to model and to assess experimentally the structural impact on different media. A variety of analytic models and numerical simulations are developed comparing to experimental results. The first part of this work presents a discussion on the similitude between a water impact and an impact on a deformable solid structure. Water impact simulations of a deformable cylinder (without rupture) are performed by finite elements (FE, Coupled Eulerian Lagrangian) and SPH analysis. An analytical model of water impact is proposed for the prediction of peak force evolution. The analysis of results permits to design an impact programmer reproducing this peak force. FE longitudinal impact simulations on cylindrical tubes, with an adapted geometry, are performed and compared with some experiments. The “dynamic buckling” of tubes under impact (due to the material inelastic behavior and to strain waves) is observed. The second part deals with the low velocity perforation (< 10 m/s, strain rate < 1000 s-1) of thin plates. Some experiments on an instrumented drop test (force, displacement, plate shape, crack propagation) are analyzed. Shell FE simulations, with a damage rupture criteria implemented are performed. Parameters are identified by inverse method with the help of Charpy tests made on 2024 T3 aluminum alloy. An analysis of the peak force, during impact, leads to a good understanding of the perforation mechanism. In parallel, a new analytical model, based on an energetic approach of the perforation, is proposed and compared with FE simulations. The numerical perforation study is extended to high velocities and high strain rates (100 - 1000m/s, strain rate < 100 000 s-1) in order to identify different well-known transitions of perforation (Petalisation, petals' fragmentation, total plate's fragmentation).

La simulation numérique par la méthode des éléments finis (MEF) fournit de nos jours une grande aide pour les ingénieurs dans les processus de conception d’optimisation des produits. Malgré le développement croissant des ressources de calcul, la fiabilité et l’efficacité des simulations numériques par la MEF restent à améliorer. Ce travail de thèse consiste à développer une famille d’éléments solide-coque (SHB) pour la modélisation tridimensionnelle des structures minces. Cette famille d’éléments SHB est basée sur une formulation tridimensionnelle en grands déplacements et rotations. La technique dite “d’intégration réduite dans le plan”, en utilisant un nombre arbitraire de points d’intégration dans la direction de l’épaisseur, permet la modélisation des structures minces avec une seule couche d'éléments. Dans ce travail de thèse, deux éléments linéaires SHB prismatique et hexaédrique, ainsi que leurs contreparties quadratiques, ont été implantés dans le code par éléments finis ABAQUS pour l’analyse quasi-statique et dynamique des structures minces. La performance de ces éléments a été validée à travers une série de cas tests académiques, ainsi que sur des problèmes complexes de type impact/crash et des procédés de mise en forme de tôles minces. L'ensemble des résultats numériques obtenus révèle que les éléments SHB représentent une alternative intéressante aux éléments coques et solides traditionnels pour la modélisation tridimensionnelle des structures minces. / Nowadays, the finite element (FE) simulation provides great assistance to engineers in the design of products and optimization of manufacturing processes. Despite the growing development of computational resources, reliability and efficiency of the FE simulations remain the most important features. The current work contributes to the development of a family of assumed strain based solid-shell elements (SHB), for the modeling of 3D thin structures. Based on reduced integration and special treatments to eliminate locking effects and to control spurious zero-energy modes, the SHB solid‒shell elements are capable of modeling most thin 3D structural problems with only a single element layer, while describing accurately the various through-thickness phenomena. In the current contribution, a family of prismatic and hexahedral SHB elements with their linear and quadratic versions have been implemented into ABAQUS using both standard/quasi-static and explicit/dynamic solvers. The performance of the SHB elements is evaluated via a series of popular benchmarks as well as with impact/crash and sheet metal forming processes. All numerical results reveal that the SHB elements represent an interesting alternative to traditional shell and solid elements for the 3D modeling of thin structural problems.

La détermination précise du comportement plastique des tôles métalliques anisotropes est un élément clé d'une simulation numérique fiable des procédés de mise en forme par emboutissage. Dans cette étude, nous avons mis au point une procédure d'identification paramétrique d'un modèle de surface de plasticité anisotrope à 8 paramètres, qui s'appuie sur des essais classiques de traction uniaxiale selon différentes orientations de l'éprouvette, mais aussi sur des essais de traction biaxiale hétérogènes. Ces derniers ont d'abord nécessité la conception de 2 types d'éprouvettes cruciformes qui, sollicitées en traction biaxiale, sont soumises à des champs de déformation couvrant le domaine allant de la traction uniaxiale à la traction équibiaxiale. L'analyse des essais de traction permet tout d'abord de déterminer les paramètres d'écrouissage ainsi que certains des paramètres de la fonction de charge liés à l'anisotropie en déformation et l'anisotropie en contrainte des matériaux. Les autres paramètres de la fonction de charge, qui interviennent dans la définition de la forme de la surface de plasticité dans le domaine de l'expansion, sont ensuite obtenus à l'aide d'une méthode d'identification minimisant l'écart entre les champs de déformation obtenus expérimentalement par une méthode d'analyse d’images, et ceux déterminés par simulation des essais à l'aide d'un code de calcul par éléments finis. La procédure d'identification a été appliquée à deux nuances d'aciers pour emboutissage, un acier inoxydable AISI304, et deux alliages d'aluminium. Les résultats montrent la grande sensibilité de la méthode d'identification proposée. Des essais de validation sont également présentés / The precise determination of the plastic behaviour of anisotropic sheet metals is a key element for obtaining reliable results in numerical simulations of forming processes. In this study, an identification procedure has been developed using an 8- parameter yield surface model. The procedure is based on the analysis of both classical uniaxial tension tests performed along different directions, and heterogeneous biaxial tensile tests. Two types of cruciform specimens have been designed to this end. Under biaxial stretching, the specimens are submitted to strain fields covering the range from uniaxial tension to equibiaxial tension. First, the analysis of uniaxial tensile tests allows us to determine strain-hardening parameters, and yield function parameters related to strain- and stress-anisotropy. Then, the other yield function parameters defining the shape of the yield surface in the biaxial stretching range are obtained using an identification procedure which minimizes the difference between strain fields obtained experimentally by an image correlation method, and strain fields determined by numerical simulations using a finite element code. The identification procedure has been applied to 2 steels of deep-drawing quality, an AISI304 stainless steel and two aluminium alloys. The results display the high sensitivity of the proposed identification method. Validation tests are also presented

Cette étude vise à bien comprendre puis à modéliser le comportement mécanique dans une large plage de vitesse de déformation des tôles d’acier austénitique AISI304 ayant subis un traitement d'attrition mécanique de la surface (SMAT). Ces tôles ainsi traités sont des matériaux multicouches avec un gradient de propriétés. Les principaux résultats obtenus sont résumés comme suit：(1) La sensibilité globale à la vitesse déformation des tôles d’acier austénitique AISI304 traités avec SMAT est caractérisée par des essais de double cisaillement sous chargements quasi-statiques et dynamiques, qui permet d’atteindre une grande déformation sans instabilité géométrique. Des essais de double cisaillement sous impact sont réalisés à l’aide des barres de Hopkinson de grande diamètre et un système d’attache qui a une même impédance acoustique que la barre. Une sensibilité significative a été révélée et on observe ce renforcement n’a pas induit une réduction importante de la ductilité.(2) Dans le but d’un meilleur dépouillement de ces essais de double cisaillement, leur conditions d’essai est analysé dans le détaillé. Le modèle numérique avec le système d’attache a été construit pour étudier l’influence du système d’attache au début de chargement. On trouve un effet limité pour les diverses conditions imparfaites des essais comme la souplesse de système d’attache, des champs mécaniques non-homogènes, l’état de non-équilibre, etc. Par contre, les études numérique et analytique ont démontré que l’hypothèse simple de petites perturbations habituellement utilisé pour le dépouillement de ces essais n’est pas suffisamment précise. La déformation Eulérien cumulée doit être utilisée pour obtenir un résultat numérique correct. A partir de ce résultat, la sensibilité à la vitesse déformation des tôles d’acier austénitique AISI304 traités avec SMAT obtenue expérimentalement a été retouchée.(3) Un modèle multicouche elasto-plastique en dommageable a été proposé pour décrire le comportement des tôles d’acier austénitique AISI304 traités avec SMAT. Les paramètres sont identifiées à partir des essai de traction. La partie elasto-plastique est calée par une loi d’écrouissage de type Ludwig. Par contre, les paramètres d’endommagement sont obtenus avec une méthode d’identification inverse sur la base de simulation numérique de ces essais de traction. Pour valider ce modèle multi-couche elasto-plastique dommageable, un essai d’indentation/perforation est réalisé sur des tôles d’acier austénitique AISI304 traités avec SMAT. Des simulations numériques correspondantes montres que ce modèle multi-couche elasto-plastique en dommageable une prédiction plutôt précise de ces essais de d’indentation/perforation.(4) Pour évaluer la performance anti-perforation des tôles d’acier austénitique AISI304 traités avec SMAT, des essais de perforation sous impact a été réalisés avec des barres de Hopkinson. Des simulations numériques de ces essais de perforation sous impact sont réalisées avec un modèle numérique comparable avec le cas quasi-statique. ́tant donne que la sensibilité globale à la vitesse déformation des tôles d’acier austénitique AISI304 traités avec SMAT est caractérisée par des essais de double cisaillement, la sensibilité à la vitesse a été introduite dans le modèle multi-couche elasto-plastique en dommageable. Le résultat numérique correspond bien à la mesure expérimentale, ce qui indique non seulement l’efficacité du modèle numérique mais aussi celle du modèle multicouche elasto-plastique en dommageable. / This Ph.D dissertation aimed at the comprehensive understanding and the constitutive modeling of the mechanical behaviours of the surface mechanical attrition treatment (SMAT) treated AISI304 stainless steel sheet under a large range of loading rates. SMAT treated AISI304 stainless steel sheets are multi-layered functionally graded materials (FGM). The main research results and conclusions are summarized as followed:(1) The overall rate sensitivity SMAT treated AISI304 stainless steel sheet is characterized by the double shearing test under quasi-static and dynamic loading where a large strain can be achieved without geometry instability. Impact double shear test are performed with a large diameter Hopkinson bar system and an adapted equal-impedance clamping device. Significant rate sensitivity is found. It is also observed that such a rate enhancement does not induce an important reduction of the ductility.(2) In order to extract accurate material information from the double shear tests, their testing conditions are thoroughly analyzed using numerical simulation. Numerical models including clamping devices have been built to investigate the influence of this clamping device at the early stage of loading. A limited effect was found for various imperfect testing conditions such as the clamping device stiffness, non-homogeneous stress and strain fields, non-equilibrium state, etc. On the contrary, numerical and analytical study shows that the simple small strain assumption usually used in double shear tests are not accurate enough. Eulerian cumulated strain definition should be used to get consistent numerical results. From this finding, the experimental rate sensitivity obtained for the SMAT treated AISI304 stainless steel sheet are recalculated.(3) A multi-layers elastic plastic damageable constitutive model is proposed to model SMAT treated AISI304 stainless steel sheet. The parameters are identified using tensile testing results. The elastic plastic behavior is curve fitted with a simple Ludwig hardening model. However, the damage parameters should be identified using an inverse method on the basis of numerical simulation of these tensile tests. In order to validate this multi-layer elastic plastic damageable constitutive model, indentation/piercing tests on SMAT treated AISI304 stainless steel sheet are performed. Numerical simulation of this indentation/piercing tests is also realized. It is found that the identified multi-layer elastic plastic damageable constitutive model allows for a quite accurate prediction of the experimental piercing tests.(4) In order to evaluate the impact anti-piercing capacity of the SMAT treated AISI304 stainless steel sheet, the impact perforation tests using Hopkinson bar are carried out. Numerical simulation of these impact perforation tests are realized with a similar FEM model as the quasi-static case. As the rate sensitivity of SMAT treated AISI304 stainless steel sheet is experimentally characterized with double shear test, a rate sensitive multi-layer elastic plastic damageable constitutive model is introduced. The numerical results agree well with the experimental ones, which indicates the effectiveness of the numerical model as well as the rate sensitive multi-layer elastic plastic damageable constitutive model.