L'objectif principal de cette étude consiste à réduire les temps de calcul des simulations de procédés de mise en forme continus sous le logiciel Forge3®. Ces procédés, tel que le laminage et le tréfilage, sont caractérisés par des pièces dont la longueur est très importante en comparaison des dimensions de la section ainsi que de la zone de contact. Une approche incrémentale générale implique des temps de calcul conséquents allant de quelques heures à plusieurs jours. En se concentrant sur le régime quasi permanant de ces procédés, une formulation stationnaire est développée pour accélérer leur simulation. Le domaine de calcul correspond initialement à une estimation de la forme de l'écoulement solution au voisinage des outils. Une étape de correction du domaine est ajoutée en plus du calcul stationnaire de l'écoulement. Comme les conditions aux limites sont modifiées, ces deux étapes sont répétées jusqu'à convergence.L'étude s'est concentrée principalement sur l'étape de correction du domaine correspondant à la résolution d'un problème de surface libre par la méthode des éléments finis. Le caractère purement convectif du problème ainsi que la prise en compte du contact nécessite l'utilisation de formulations faibles faisant apparaître un décalage amont (SUPG). Deux nouvelles formulations basées sur la méthode des moindres carrés sont développées avec succès (MC_supg et MC_lc). Pour appliquer la méthode à des géométries complexes, différentes méthodes de généralisations sont développées où un 2ème degré de liberté est ajouté aux nœuds de surface. La méthode la plus performante (CSL_dif) consiste à utiliser ce 2ème degré de liberté pour le calcul de surface libre uniquement sur les nœuds appartenant à une arête géométrique, pour les autres une régularisation du maillage dans la direction tangente y est effectuée. Des résultats excellents ont été observés sur un grand nombre de cas tests analytiques. Le contact est appliqué par une méthode de pénalisation aux nœuds. Afin de renforcer le couplage entre cette étape et celle du calcul de l'écoulement, un contact bilatéral glissant est attribué aux nœuds en compression alors que pour les autres nœuds un contact unilatéral est employé. Un algorithme spécifique est développé pour déterminer avec précision la zone de contact.Cette formulation itérative pour la recherche de l'état stationnaire a été appliquée avec succès sur un grand nombre de cas tests de mise en forme. Des accélérations comprises entre 10 et 60 ont été obtenues par rapport à Forge3®.

Les AMF à base fer se distinguent par la présence d'un couplage entre la transformation de phase et le glissement plastique. Pour caractériser leur comportement thermomécanique, des essais cycliques de chargement mécanique effectués à différentes températures et à différents niveaux de chargement, suivis par chauffage, ont été effectués. Ceux ci nous ont permis de distinguer les spécificités des AMF à base fer comme la non linéarité des interactions inter et intragranulaires, l'évolution des contraintes critiques avec la température, l'effet de la déformation plastique sur celle de transformation et l'activation de la transformation inverse uniquement par chauffage. En s'inspirant des résultats de ces essais, une loi de comportement, intégrant toutes ces spécificités, est proposée. Elle dérive d'une expression de l'énergie libre de Gibbs issue de considérations micromécaniques. Elle présente deux variables internes, la fraction volumique de martensite pour décrire la transformation de phase et le taux de plasticité pour le comportement plastique. Cette loi a été implémentée dans le code éléments finis Abaqus via la subroutine UMAT. Elle a été validée par comparaison des simulations numériques avec les résultats expérimentaux. Suite à cette validation, des applications en AMF à base fer ont été étudiées. Les résultats obtenus ont montré la capacité du modèle à être utilisé comme outil de dimensionnement de structures en AMF à base fer / The Fe-based shape memory alloys (SMAs) present a specific thermomechanical behaviour compared with classical SMAs. In this PhD thesis, experimental thermomechanical tests were performed in order to study such behavior. The applied loading is a tension followed by a significant heating. The loading cycle is repeated at different constant temperatures and maximum stresses. The experimental results show a coupling between two non linear inelastic mechanisms: phase transformation and plasticity. The reverse transformation activated only during heating, the effect of plastic strain on the transformation one, were also analysed. Taking into account all these specificities, a finite element numerical tool adapted to Fe-based SMA structural analysis is proposed. It is based on a developed constitutive model which describes the effect of phase transformation, plastic sliding and their interactions on the thermomechanical behavior. Two scalar internal variables were considered to describe phase transformation and plastic sliding effects. This model was derived from an assumed expression of the Gibbs free energy taking into account, in addition to mechanical and chemical quantities, the non linear interaction quantities related to inter- and intra-granular incompatibilities. The numerical tool derived from the implicit resolution of the non linear partial derivative constitutive equations was implemented in the Abaqus finite element code via the UMAT subroutine. After verification tests for homogeneous and heterogeneous thermo-mechanical loadings, two examples of Fe-based SMA applications were studied. They correspond to Fe-based SMA tightening systems: a fish plates for crane rails and a ring for tubes connection

La fissuration à chaud est un défaut majeur rencontré en solidification des alliages d'aluminium. Elle est liée à l'incapacité du liquide de s'écouler dans les zones où des porosités sont présentes, ne permettant pas de les refermer avant qu'elles gagnent en volume. Pour comprendre la fissuration à chaud, il est crucial de développer nos connaissances du comportement mécanique de la zone pâteuse. Pour cela, il est très utile d'effectuer des expériences de microtomographie aux rayons X et des simulations mécaniques sur des volumes élémentaires représentatifs. Dans cette thèse, nous proposons de coupler les deux approches en initialisant une simulation par éléments finis grâce à des données de microtomographie issues d'un test de traction isotherme d'un alliage d'aluminium-cuivre à l'état pâteux. Cette approche originale nous donne directement accès à la réalité expérimentale et permet des comparaisons des évolutions numérique et expérimentale de l'éprouvette. Nous expliquons dans un premier temps comment obtenir la représentation numérique à l'aide de l'algorithme des marching cubes et de la méthode d'immersion de volume. Nous présentons ensuite notre modèle numérique qui s'appuie sur une résolution monolithique des équations de Stokes. Une fois le champ de vitesse obtenu dans l'ensemble des phases solide, liquide et gazeuse, nous utilisons une méthode level set dans un formalisme eulérien afin de faire évoluer la morphologie de notre échantillon numérique. Malgré la simplicité du modèle, les résultats expérimentaux et numériques montrent un accord raisonnable en ce qui concerne la propagation de l'air à l'intérieur de l'échantillon. / Hot tearing is a major defect arising during solidification of aluminium alloys. This defect is associated with the inability of liquid to feed areas where voids have started to appear, not allowing to heal small defects before they grow bigger. To understand hot tearing, it is mandatory to develop a good knowledge of the semi-solid mechanical behaviour. It is thus very useful to carry out X-ray microtomographies experiments and mechanical simulations on representative elementary volumes. In this work, we couple the both approaches by initialising a finite element simulation with the help of microtomography data obtained during an isothermal tensile testing of an aluminium-copper alloy in the mushy state. This innovative approach gives a direct access to the experimental reality and allows comparisons of numerical and experimental evolutions of the sample. We explain in a first time how to get the numerical representation thanks to a marching cubes algorithm and the immersed volume method. Then, we present our numerical model for which we solve the Stokes equations in a monolithic way. Once the velocity computed in all the solid, liquid and gaseous phases, we use a level set method in a Eulerian formalism to obtain the morphological evolution of our numerical sample. Despite the model simplicity, numerical and experimental results show a reasonable agreement concerning the air propagation inside the sample.

Les vibrations des systèmes de forage pétrolier sont à l'origine de nombreux dysfonctionnements (ruptures des tiges par une fatigue accélérée, réduction des performances, endommagement des outils de mesures, endommagement des parois du puits, etc.). Face à la complexité des puits forés aujourd'hui, la maîtrise des vibrations des systèmes de forage est plus que jamais un enjeu majeur dans la réussite économique d'un projet pétrolier. Durant l'opération de forage, les tiges en rotation entrent en interaction avec les parois du puits (tubage et/ou formation) et encaissent dans certains cas des vibrations sévères. On distingue généralement trois modes de vibrations suivant le plan de leur occurrence : axiales, latérales et de torsion. Nous ne nous intéressons dans ce mémoire qu'aux vibrations latérales des tiges de forage. Si les vibrations latérales ont fait l'objet de nombreuses études dans le passé, il reste néanmoins des axes d'amélioration possible, tant sur la compréhension des phénomènes (contact garniture-puits par exemple) que sur la recherche de méthodes numériques permettant de réduire les temps de calcul. Dans le cadre de cette thèse, un modèle a été développé pour étudier les vibrations latérales des garnitures de forage dans des forages à trajectoires complexes. Ce modèle permet de prédire les vibrations latérales des tiges pour des paramètres opératoires donnés (vitesse de rotation, poids sur l'outil de forage). Un modèle numérique en éléments finis a été développé pour résoudre les équations du mouvement et analyser ainsi la sensibilité des vibrations aux paramètres opératoires du forage en particulier la vitesse de rotation et l'effort axial dans les tiges. Le modèle permet en outre d'analyser la réponse dynamique d'une garniture en cours du forage (conception). En outre, cette étude a permis de mieux élucider le phénomène sévère de précession des tiges (whirling), très nuisible à l'intégrité mécanique des systèmes de forage. Un nouveau banc d'essais a été mis au point par le Centre de Géosciences de Mines ParisTech pour reproduire les vibrations latérales, mieux comprendre le phénomène du whirling et valider les résultats numériques du modèle. Par ailleurs, des mesures dynamiques en surfaces et en fond de puits au cours de forages réels ont été analysées afin de mettre en évidence les vibrations latérales les plus sévères et en particulier le whirling. Ces données de terrain ont permis de comparer les fréquences propres du système mesurées et celles fournies par le modèle numérique. / Drillstring vibrations are commonly observed during oil & gas well drilling operations. Vibrations are a major cause of drilling tools dysfunction (drillstring breaking because of fatigue, reduced drilling efficiency, measurement-while-drilling tools failure, damaging of drill bits, etc.). Because of the increasing complexity of oil & gas wells drilled nowadays, operators need to mitigate efficiently the drillstring vibrations in order to successfully achieve the drilling process. During the drilling operation, rotating drillstrings are in interaction with the well borehole (casing and/or rock) which may lead to severe vibrations. Different vibrations modes occur simultaneously while drilling, we identify mainly three modes: axial, torsional and lateral. This work deals only with lateral vibrations. Literature survey papers show numerous experimental and numerical studies carried out on drillstring dynamics. The developed models don't take into account sufficiently the complex drillstring-borehole interactions or the efficient numerical methods needed to reduce the computation time. A new drillstring dynamics model has been developed within this thesis in order to compute the lateral vibrations of drillstrings in a complex well trajectory. Given the operating parameters (rotary speed, weight on bit) the model predicts the dynamics response of the drillstring in terms of lateral vibrations. A finite element model has been implemented to solve for the equations of motion of the dynamics model and study the dependence of the lateral vibrations on some operating parameters of the system, mainly the rotary speed and the axial load on the drillstrings. The finite element model can be used to compute and enhance the dynamic response of a given drillstring configuration for design issues. Besides, the model can be used to understand some dynamic phenomena encountered while drilling (post-analysis). Moreover, this study was useful to better understand the “whirling” phenomenon which is very harmful for the drilling system components. A new lateral vibrations simulator was built at Mines ParisTech in order to understand the whirling phenomenon and validate the numerical results provided by the dynamics model. Surface and downhole fieldmeasurements have been analyzed in order to understand the occurrence of whirling. The eigenfrequencies evaluated from the field data have been found very close to those provided by the dynamics model.

L'application du procédé SLM est limitée par la difficulté à contrôler le procédé. Son application aux céramiques est particulièrement difficile en raison de leur faible absorption au laser et de leur faible résistance au choc thermique. La maîtrise de ce procédé nécessite une compréhension complète du transfert de chaleur, de la dynamique des fluides et de la mécanique des solides. Dans ce travail, nous proposons un modèle numérique pour la simulation du procédé SLM appliqué aux céramiques. Le modèle est développé à l'échelle du cordon et avec l'hypothèse d'un lit de poudre continu. Il est basé sur la méthode level set et l'homogénéisation multiphasique, avec laquelle nous sommes capables de suivre l'évolution de l'interface gaz/matière et les transformations de phase. La simulation dévelopée permet d'étudier l'influence des propriétés du matériau et des paramètres du procédé sur la température, la forme du bain liquide, la dynamique des fluides et la mécanique des solides. En dehors de la puissance du laser et de la vitesse de balayage, l'absorption du matériau est également importante pour la thermique et la forme du bain liquide. Avec la dynamique des fluides, la forme convexe du cordon est obtenue sous tension de surface. Les gouttelettes liquides se forment lors de la fusion de la poudre et créent une instabilité du bain. Ceci entraîne une irrégularité du cordon après solidification. L'effet Marangoni, provoqué par le gradient surfacique de la tension de surface, est étudié. Son influence sur la répartition de la température, la forme du bain liquide et la régularité du cordon est évoquée. Cet effet peut lisser la surface du cordon avec ∂γ/∂T négatif. En augmentant la vitesse de balayage, la surface du cordon devient plus irrégulière. L'effet de « balling » est reproduit avec une vitesse de balayage élevée. Cela peut être utile pour trouver le régime donnant une forme de cordon régulière étant données la puissance et la vitesse du laser. Le défaut de fissuration est délétère dans la fabrication additive. L'utilisation d'un laser auxiliaire peut aider à éviter ce défaut en diminuant la contrainte de traction maximale. Le mode de fonctionnement de ce laser auxiliaire reste un sujet intéressant à étudier et quelques pistes ont été données par les simulations présentées. Le modèle est validé par la comparaison de la forme du bain liquide avec des expériences dans différentes conditions de procédé. Les simulations peuvent également révéler la tendance de variation de la surface du cordon dans certains cas. Par la simulation de la déposition de cordons multiples, l'influence de taux de recouvrement sur la surface d'une couche, la température et l'évolution de contrainte est soulignée. / The application of SLM process is limited by the difficulty of process control. Its application to ceramics is especially challengeable due to their weak absorption to laser and weak resistance to thermal shock. The mastery of this process requires a full understanding of heat transfer, fluid dynamics in melt pool and solid mechanics. In this work, we propose a numerical model for the simulation of SLM process applied to ceramics. The model is developed at the track scale and with the assumption of continuous powder bed. It is based on level set method and multiphase homogenization, with which we are able to follow the evolution of gas/material interface and phase transformation. Simulations are performed to study the influence of material properties and process parameters on temperature, melt pool shape, fluid dynamics and solid mechanics. Apart from the laser power and scanning speed, material absorption is also found to be important to the thermal behavior and the melt pool shape. With the fluid dynamics, convex shape of track cross section is achieved under surface tension. Besides that, liquid droplets collapsing formed by the melting of powder create melt pool instability when falling, thus leading to track irregularity after solidification. The Marangoni effect, caused by surface tension gradient at gas/material interface, is investigated. Its influence on temperature distribution, melt pool shape and track regularity is recognized. One interesting finding is the smoothing effect of track surface with negative ∂γ/∂T. When combine surface tension with scanning speed, track surface becomes more irregular with the increase of scanning speed. The well-known balling effect is reproduced with high scanning speed. This can be helpful to find the regime for regular track shape with given laser power and scanning speed. Cracking defect is deleterious in additive manufacturing. The use of an auxiliary laser can help to avoid this defect by decreasing the maximum tensile stress. The process mode of this auxiliary laser remains an interesting subject to be studied and some guidelines have been given by the presented simulations. The model is validated by the comparison of melt pool shape with experiments under different process conditions. Simulations can also reveal the tendency of track surface variation for certain cases. By the application to multi-track deposition, the influence of hatch distance on layer surface, temperature and stress evolution is emphasized.

Ce travail présente une méthode mixte expérimentale-numérique traitant des mesures mécaniques réalisées sur structures pour l’identification de propriétés matériau viscoélastiques linéaires. Les caractéristiques de ces matériaux sont complexes et dépendent de la fréquence et nous cherchons à les identifier par une approche directe et non paramétrique. Pour cela, deux méthodes d’identification adaptées au traitement des mesures de champs de réponses en fréquences sont présentées. Cellesci utilisent les données fréquentielles, sans analyse modale expérimentale. Les champs expérimentaux sont obtenus par vibromètrie laser et les champs numérique par la méthode des éléments finis. Nous utilisons un modèle volumique adapté à l’étude des structures minces et épaisses pour l’analyse numérique. Les applications se basent sur des plaques simples ou multicouches. Les résultats sont corrélés avec des mesures DMA et des identifications à partir de données modales identifiées. / This work deals with a mixed experimental numerical method for linear viscoelastic material parameter identification. These parameters are complex and frequency dependent. Identication is performed using a direct and non parametric approach. In this goal, two identification methods adapted to frequency response fields measurements treatment are presented. These datas are durectly exploited, without any experimental modal analysis. Experimental fields are obtained using a vibrometer and numerical ones with finite element method. A brick element dedicated to thin/thick structures analysis is employed. Appications are performed using simple plates or sandwich plates. Results are also compared to DMA measurements and identified modal parameters.

Dans de nombreuses applications industrielles, les structures en vibration à dimensionner sont en contact avec un fluide (fluide autour des coques des bateaux, réservoirs, échangeurs de chaleur dans les centrales, l’industrie automobile, etc). Cependant, le comportement dynamique de la structure peut être modifié de façon importante par la présence du fluide. Le dimensionnement doit donc prendre en compte les effets de l’interaction fluide-structure.Ces applications nécessitent un couplage efficace. En outre, l’analyse dynamique des systèmes industriels est souvent coûteuse du point de vue numérique. Pour les modèles éléments finis des problèmes couplés fluide-structure, l’importance de la réduction de la taille devient évidente car les degrés de liberté du fluide seront ajoutés à ceux de la structure. Des méthodes de réduction du modèle seront utilisées pour réduire la taille des matrices obtenues.Traditionnellement, l’étude de ces systèmes couplés est fondée sur une démarche déterministe dans laquelle l’ensemble des paramètres utilisés dans le modèle prennent une valeur fixe.Par contre, il suffit d’avoir procédé à quelques expérimentations pour se rendre compte des limites d’une telle modélisation, d’où la nécessité de la prise en compte des incertitudes sur les paramètres du système couplé.Ce travail de thèse s’articule autour de trois études principales. La première consiste à mener une étude déterministe numérique et analytique des problèmes vibro-acoustiques sans réduction de modèles. Cette dernière est basée sur une formulation non symétrique déplacement/pression et une formulation symétrique déplacement/pression et potentiel des vitesses. Dans la deuxième étude, on propose deux méthodes de réduction du modèle : analyse et synthèse modales pour la résolution des problèmes vibro-acoustiques des grandes tailles des systèmes couplés modélisés par la méthode des éléments finis. La méthode de synthèse modale développée couple une méthode de sous-structuration dynamique de type Craig et Bampton et une méthode de sous domaines acoustiques.Enfin, pour tenir compte des incertitudes sur les paramètres du système couplé, on a développé dans la troisième étude une méthode numérique stochastique de synthèse modale étendue à une étude de fiabilité basée sur les approches FORM et SORM pour la résolution de ces problèmes. Ces démarches vont nous permettre de résoudre les problèmes vibro-acoustiques, sans utiliser les méthodes classiques, qui consistent à faire un calcul modal direct allié à la simulation de Monte Carlo demandant un coup de temps très élevé.Plusieurs exemples académiques et industriels ont été traités pour valider les approches proposées.L’étude numérique est conduite en utilisant un code élaboré sous MATLAB couplé au code commercial ANSYS afin d’évaluer la fiabilité du système couplé. La confrontation des résultats numériques, analytiques et expérimentaux nous permet de valider conjointement le processus de calcul et les méthodes proposées dans le domaine de l’analyse fréquentielle et l’étude fiabiliste des structures immergées. D’un point de vue industriel, ces méthodes visent à promouvoir l’introduction de la culture de l’incertain dans les métiers de la conception et encouragent la construction d’un modèle fiable et robuste pour les problèmes d’interaction fluide-structure. / In several industrial applications, the vibrating structures are in contact with a fluid (fluid around the hulls of a boats, reservoirs, heat exchangers in power plants, ...), but the dynamic behavior of the structure can be significantly modified by the presence of the fluid. The sizing must take into account the effects of fluid-structure interaction. These applications require an effective coupling. In addition, the dynamic analysis of the industrial systems is often expensive from the numerical point of view. For the coupling fluid structure finite elements models, the importance of the size reduction becomes obvious because the fluid’s freedom degrees will be added to those of the structure. A proposed condensation method will be used to reduce the matrixes size. Traditionally, the study of the fluid-structure interaction is based on a deterministic approach where all the parameters used in the model have a fixed value. But it suffices having conducted a few experimentations to realize the limitations of such modeling. Hence it need to take into accounts the uncertainty on the parameters of mechanical systems. In this thesis, we deal with the simulation of vibro-acoustic problems. The first part presents a numerical and analytical study of deterministic problems without model’s reduction, based on a non-symmetric formulation displacement/pressure and on a symmetric formulation displacement/pressure and velocity potential. In the second part of this work, two methods are proposed to reduce the model : modal analysis and modal synthesis for solving vibro-acoustic problems of large sizes modeled by finite elements method. The developed modal synthesis method is coupling dynamic substructure of Craig and Bampton type and acoustic subdomain based on a pressure formulation. To take into account the parameter’s uncertainties of the coupled system, we have developed a numerical stochastic method of the modal synthesis and modal analysis extended to reliability study, based on the FORM and SORM approaches. These approaches will allow us to solve the vibro-acoustic problems without using classical procedure. It may become prohibitive in terms of computation time. Several academic and industrial examples are studied to validate the proposed methods. The numerical study is performed using a code developed with MATLAB coupled with the commercial code ANSYS in order to evaluate the reliability of systems. The comparison of numerical, analytical and experimental results enables us to jointly validate the calculation process and the proposed methods in the domain of frequency analysis and reliability study of submerged structures. From the industrial point of view, our research work aim is to promote the introduction of the uncertainty’s culture during modeling in the context of design processes.

Dans le cadre des procédés de chauffage et de trempe réalisés lors d’opération de mise en forme des matériaux, le rayonnement thermique joue un rôle prépondérant. Lors de l’élaboration de modèles numériques permettant la simulation de ces procédés, il est donc nécessaire de disposer d’outils performants pour simuler ce phénomène.La simulation numérique de tels procédés soulèvent de nombreuses problématiques, comme la représentation d’un environnement complexe impliquant plusieurs composants (pièces, bruleurs, buses d’injection, parois), la gestion des divers phénomènes physiques couplés (écoulement, transfert thermiques, ébullition, rayonnement). Dans cette perspective, les méthodes dites « d’immersion », permettant un traitement généraliste de ces divers problèmes, rencontrent depuis quelques années un intérêt grandissant dans la communauté scientifique.C’est dans ce contexte que s’inscrit le projet Thost, au sein duquel est réalisée cette thèse. L’objectif est donc de développer des outils pour la modélisation du rayonnement dans le contexte d’immersion de volume propre au logiciel Thost. Deux approches sont développées : l’une consistant en l’adaptation d’une méthode existante au contexte de l’immersion de volume, l’autre explorant l’élaboration d’une formulation pour un modèle particulier de rayonnement. Les outils développés sont ensuite mis à l’épreuve sur des simulations de cas industriels fournis par nos partenaires. / For heating and quenching operations occurring during material forming processes, thermal radiation is the the predominant physical phenomenon. Hence, when one tries to simulate such processes, it is important to have at disposal powerful tools for the numerical modelling of thermal radiation.The numerical simulation of these processes often rises numerous problems and questions, as the representation of a complex environment, involving several components ( ingots, burners, nozzles, walls), to deal with different coupled physical phenomena ( flow, heat transfer, boiling, thermal radiation). In this regard, some “immersed” numerical methods, allows a generalist treatment of these different problems, have gained popularity and drag interest of the scientific community in the recent years.The Thost project, aiming to produce a software for heat transfer during material forming processes, fits in the framework, and this PhD is part of this project. The goal is therefore to design tools for numerical modelling of thermal radiation within the immersed volume method of the Thost software. Two approaches are presented: one consisting in the adaptation of an existing method to the context of the immersed volume method, another concerning the development of a formulation for a specific model of radiation. These methods are then tested on industrial applications provided by our partners.

De nombreuses structures industrielles soumises à des chargements à long terme statique (fluage) ou cyclique (fatigue) sont constituées de matériaux polymères semi-cristallins. C’est le cas notamment des canalisations et réservoirs sous pression. Il est donc essentiel de traiter les problématiques de durabilité pour être capable d'anticiper et de contrôler leur fin de vie. Par ailleurs, elles présentent généralement des formes complexes et sont soumises à des états de contrainte multiaxiaux.Le matériau de l'étude est un polymère semi-cristallin : le Polyamide 6. Il est caractérisé par la coexistence d'une phase cristalline et d'une phase amorphe qui s'arrangent selon une microstructure sphérolitique.Dans un premier temps, les liens entre comportement mécanique à l'échelle globale de l'éprouvette et les micro-mécanismes de déformation sous-jacents conduisant à la rupture sont établis expérimentalement pour des sollicitations en traction monotone et en fluage qui présentent des résultats similaires puis en fatigue. L'influence de la multiaxialité de l'état de contrainte est étudiée à partir d’éprouvettes axisymétriques entaillées de différents rayons de fond d'entaille et d'éprouvettes « Compact Tensile ». Les phénomènes de cavitation sont caractérisés grâce aux techniques de tomographie et laminographie à rayonnement X synchrotron qui permettent d'observer et de quantifier les distributions spatiales de taux de porosité volumique et le caractère anisotrope des cavités. Et l'analyse des faciès de rupture a permis de mettre en évidence que les mécanismes de croissance et de coalescence de cavités étaient à l'origine de l’amorçage ductile de la rupture.Ensuite, un modèle poro-visco-plastique à deux mécanismes (permettant de différencier le comportement des phases amorphe et cristalline) a été utilisé. Ce modèle permet de reproduire à la fois le comportement global (courbes de chargement) en traction monotone et en fluage mais aussi les distributions spatiales de taux de porosité obtenues expérimentalement. De plus, les calculs par éléments finis permettent d'étudier les distributions spatiales du champ de contrainte et d'établir l'influence de l'état de contrainte sur l'état de cavitation. Les évolutions temporelles en cours de déformation de la pression hydrostatique (ou contrainte moyenne) ont été reliées aux distributions spatiales de taux de porosité volumique. Et l'anisotropie de cavitation (et donc la morphologie et les facteurs de forme des cavités) a été reliée aux évolutions des composantes du tenseur des contraintes de Cauchy. Enfin, la définition d'un critère de rupture en taux de porosité critique a permis de simuler l'amorçage et la propagation de fissures en traction monotone et fluage. / Many industrial structures subjected to quasi-static (creep) or cyclic (fatigue) long-term loadings are made of semi-crystalline polymers. Such is the case, for instance, of pressure vessels and pipes. It is therefore considered critical to study the issues related to their durability in order to be able to anticipate and control their end of life. Furthermore, they generally have complex designs and are subjected to multiaxial stress states.The material which has been studied was a semi-crystalline Polyamide 6. Its structure consisted of amorphous and the crystalline phases and a spherolitic microstructure.As a first step, the links between the mechanical behaviour at the global scale of the specimens and the underlying micro-mechanisms of deformation that lead to failure have been established experimentally for monotonic and creep loadings that show similar results and then for fatigue loadings. The influence of the multiaxiality of the stress state has been studied using circumferentially notched round bars with different notch root radii and Compact Tensile specimens. The cavitation phenomena were characterized using synchrotron radiation tomography and laminography techniques that enabled the observation and quantification of the spatial distributions of the voids and the anisotropy of the cavities. An analysis of the fracture surfaces has shown that the initiation of ductile failure resulted from void growth and coalescence mechanismsA poro-visco-plastic model with two mechanisms (that allow the behaviours of the amorphous and crystalline phases to be distinguished) has been used. Thanks to this model, the global behaviour (loading curves) under steady strain rates and steady loads but also the spatial distributions of the void volume fraction could be reproduced numerically. In addition finite element calculations have permitted the spatial distributions of the stress field to be studied and the influence of the stress state on the cavitation state to be investigated. The temporal evolutions during the deformation of the hydrostatic pressure have been linked to the spatial distributions of void volume fraction. The void anisotropy (and thus the void morphology and shape factors) has been related to the evolutions of the components of the Cauchy stress tensor. Finally, the definition of a rupture criterion based on a critical value of the void volume fraction has enabled crack propagation under steady strain rate and steady load to be simulated.

Cette thèse porte sur le développement d’une méthode multigrilles parallèle visant à réduire les temps de calculs des simulations éléments finis dans le domaine de la mise en forme de pièces forgées en 3D. Ces applications utilisent une méthode implicite, caractérisées par une formulation mixte en vitesse/pression et une gestion du contact par pénalisation. Elles impliquent de grandes déformations qui rendent nécessaires des remaillages fréquents sur les maillages tétraédriques non structurés utilisés. La méthode multigrilles développée suit une approche hybride, se basant sur une construction géométrique des niveaux grossiers par déraffinement de maillage non emboîtés et sur une construction algébrique des systèmes linéaires intermédiaires et grossiers. Un comportement asymptotique quasi-linéaire et une bonne efficacité parallèle sont attendus afin de permettre la réalisation de simulations à grand nombre de degrés de liberté dans des temps plus raisonnables qu’aujourd’hui. Pour cela, l’algorithme de déraffinement de maillages est compatible avec le calcul parallèle, ainsi que les opérateurs permettant les transferts de champs entre les différents niveaux de maillages partitionnés. Les spécificités des problèmes à traiter ont mené à la sélection d'un lisseur plus complexe que ceux utilisés plus fréquemment dans la littérature. Sur la grille la plus grossière, une méthode de résolution directe est utilisée, en séquentiel comme en calcul parallèle. La méthode multigrilles est utilisée en tant que préconditionneur d’une méthode de résidu conjugué et a été intégrée au logiciel FORGE NxT et montre un comportement asymptotique et une efficacité parallèle proches de l’optimal. Le déraffinement automatique de maillages permet une compatibilité avec les remaillages fréquents et permet à la méthode multigrilles de simuler un procédé du début à la fin. Les temps de calculs sont significativement réduits, même sur des simulations avec des écoulements particuliers, sur lesquelles la méthode multigrilles ne peut être utilisée de manière optimale. Cette robustesse permet, par exemple, de réduire de 4,5 à 2,5 jours le temps de simulation d’un procédé. / A parallel multigrid method is developed to reduce large computational costs involved by the finite element simulation of 3D metal forming applications. These applications are characterized by a mixed velocity/pressure implicit formulation with a penalty formulation to enforce contact and lead to large deformations, handled by frequent remeshings of unstructured meshes of tetrahedral. The developed multigrid method follows a hybrid approach where the different levels of non-nested meshes are geometrically constructed by mesh coarsening, while the linear systems of the intermediate and coarse levels result from the algebraic approach. A close to linear asymptotical behavior is expected along with parallel efficiency in order to allow simulations with large number of degrees of freedom under reasonable computation times. These objectives lead to a parallel mesh coarsening algorithm and parallel transfer operators allowing fields transfer between the different levels of partitioned meshes. Physical specificities of metal forming applications lead to select a more complex multigrid smoother than those classically used in literature. A direct resolution method is used on the coarsest mesh, in sequential and in parallel computing. The developed multigrid method is used as a preconditioner for a Conjugate Residual algorithm within FORGE NxT software and shows an asymptotical behavior and a parallel efficiency close to optimal. The automatic mesh coarsening algorithm enables compatibility with frequent remeshings and allows the simulation of a forging process from beginning to end with the multigrid method. Computation times are significantly reduced, even on simulations with particular material flows on which the multigrid method is not optimal. This robustness allows, for instance, reducing from 4.5 to 2.5 days the computation of a forging process.