Ce travail a été effectué dans le cadre du projet Rem3D® dans lequel participent plusieurs entreprises avec l'objectif de développer un logiciel d'injection en 3D par éléments finis. L'objectif est de développer une méthode numérique pour modéliser le comportement viscoélastique des polymères de l'état solide à l'état liquide à travers une approche multiphasique qui est largement utilisé pour traiter le problème de l'interaction fluide-structure (IFS). La philosophie est d'utiliser une formulation mixte de trois champs (u, v, p) (déplacement, vitesse, pression), où u et v représentent les principales variables de déformation et de vitesse de déformation. Nous sommes amenés au problème de Navier-Stokes compressibles avec l'extra-contrainte, qui est résolu en utilisant la méthode des éléments finis mixte. Le présent travail contribue aussi certains éléments de stabilisation pour la simulation numérique des problèmes multiphasiques par l'approche monolithique.Comparaison entre la littérature et l'expérience est accompli par la validation du cas élastique et cas modèle viscoélastique de Kelvin-Voigt dans le lagrangien approche ainsi qu'eulérien approche. L'extension de la méthodologie au modèle visco-hyper-élastique est débuté par la modélisation et la validation au point matériel, puis l'implémentation dans la bibliothèque des éléments finis CimLib®. Enfin, un schéma stabilisation de résolution du type EVSS est adopté pour le modèle viscoélastique de Kelvin-Voigt, le modèle visco-hyper-élastique de Néo-Hookean, et aussi le modèle visco-hyper-élastique qui propose une prometteuse porte ouverte dans la simulation et modélisation, non seulement pour la viscoélasticité, mais aussi pour les applications dynamique complexes. / This work concerns the simulation of viscoelastic behavior of polymer at different states. Viscoelastic modeling of polymer was performed from the solid state to the liquid state via a multiphase approach which is largely used to deal with the fluid structure interaction. To ensure the appreciation of the FSI, viscoelasticity is considered in two parts: an elastic one and viscous other where the main idea is to use a mixed formulation in three fields (u, v, p) (displacement, velocity, pressure), with u and v, represented the primary variables of a strain and a strain rate formulation. We are led to the Navier-Stokes compressible problem with extra-stress, which is solved by using the Mixed Finite Element. The present work contributes some stabilization elements to the numerical simulation of multiphase problem by the monolithic approach.Comparison between the literature and experiments was performed through the validation of an elastic case and the viscoelastic Kelvin-Voigt model in the context of Lagrangian framework as well as Eulerian framework. The extension of the methodology to a visco-hyper-elastic is given through the modeling and validation on material point on the finite elements library CimLib®. Finally, a stabilization scheme of the EVSS type is adopted for viscoelastic Kelvin-Voigt model, hyper-elastic Neo-Hookean model, and also visco-hyper-elastic model which proposed an open door in computational modeling, not only with viscoelasticity but also complex dynamic application.

Le procédé de déroulage est un procédé de coupe orthogonale du bois vert dont l'arête de coupe est parallèle à la fibre et dont la valeur ajouté est apportée au copeau. L'enjeu actuel est de d'augmenter la déroulabilité du Douglas français ainsi que de déterminer celle des bois guyanais tout en évitant de trop lourdes expérimentations. Une modélisation de ce procédé est donc nécessaire. Le noeud scientifique est la loi de comportement du bois durant le déroulage, comme ce procédé est caractérisé par une grande vitesse de déformation et d'importantes et complexes déformations mêlant cisaillement et compression. Afin comprendre les phénomènes physiques mis en jeux, nous avons réalisé plus de 700 essais quasistatiques de traction, compression et torsion avec différents types de chargement. Nous avons aussi réalisé une centaine d'essais de déroulage sur une micro-dérouleuse instrumentée à laquelle nous avons ajouté des mesures de températures. Nous proposons une modélisation analytique adaptée à la description des régimes continus de coupe. La description de la zone de déformation est lagrangienne et la résolution se fait par les puissances mises en jeu. Nous proposons une loi de comportement hyper-élastique isotrope compressible de Rivlin généralisée d'ordre 4 dont nous avons déterminé la forme grâce aux essais quasi-statiques et les coefficients grâce à une résolution de problématique inverse de Levemberg-Marquardt afin de tenir compte de la vitesse de déformation. Pour la première fois, l'angle d'inclinaison de la zone de déformation prédit correspond aux observations expérimentales. De plus, au moins 75% des efforts calculés sont à moins de 15% des efforts expérimentaux.