Les travaux de cette thèse ont porté sur le développement d’une nouvelle loi de comportement hyperélastique, isotrope et incompressible permettant de modéliser les matériaux caoutchouteux en grande déformation et en grand déplacement. Cette nouvelle loi combine une approche moléculaire et une approche phénoménologique, ce qui permet de couvrir un spectre large de sollicitations. Elle est constituée par la superposition de quatre termes :– un terme lié à la contrainte d’entrelacement des chaînes macromoléculaires observée avec le phénomène de cristallisation. Ce terme est modélisé par une fonction logarithmique provenant de l’énergie phénoménologique de Gent-Thomas,– un terme lié à l’hypothèse des déformations affines observées avec le raidissement final de certaines chaînes macromoléculaires des élastomères. Ce terme provient de la probabilité non-Gaussienne de Langevin. Nous l‘avons modélisé par la loi moléculaire 8-chaines d’Arruda-Boyce avec un aménagement qui consiste à utiliser une approximation originale de la fonction de Langevin inverse,– un terme lié à la contrainte des chaînes ayant des déformations non-affines. Ce terme est modélisé par une fonction Gaussienne sous forme intégrale. Il s’agit de l’une des contributions originale de ce travail de thèse,– une partie volumique standard permettant de prendre en compte l’incompressibilité du matériau.Les deux principales originalités de la thèse concernent donc l’élaboration d’une approximation inédite de la fonction de Langevin inverse ainsi que la construction d’une nouvelle densité d’énergie hyperélastique isotrope, incompressible et hybride.Afin d’étudier la pertinence du modèle proposé, des comparaisons ont été réalisées avec plusieurs jeux de données expérimentales disponibles dans la littérature. Ces comparaisons ayant été couronnées de succès, l’implémentation numérique du modèle que nous proposons a été effectuée dans le code universitaire aux éléments finis FER. / This thesis concerns the development of a new incompressible isotropic hyperelastic behavior law allowing the modeling of rubber materials with large strain and large displacement. This new law mixes a molecular approach with a phenomenological one and therefore covers a wide range of loading. It has been built by a sum over four terms:– a term related to the interleaving macromolecular chains observed with the crystallization phenomenon. This term is modeled by a logarithmic function coming from the phenomenological energy of Gent-Thomas,– a term related to the assumption of affine deformations observed with the final stiffening of a part of macromolecular elastomeric chains. This term comes from the non Gaussian probability of Langevin. We have modeled it by the 8-chains molecular law of Arruda-Boyce but with an original approximation of the inverse of the Langevin function,– a term related to the stress occurring with non affine strains. This term has been modeled by a Gaussian function adopting an integral form. This is one of the original contribution of this thesis work,– a classical volumetric term taking into account the incompressibility of the material.The two main originalities of the thesis are therefore the introduction of a new approximation of the inverse of the Langevin function and the development of a new hyperelastic energy density which is isotropic, incompressible and hybrid.In order to study the efficiency of the proposed model, comparisons were made with several experimental data available in the literature. These comparisons have been successful and we have implemented our model in the university finite element software FER.

L'identification du comportement des polymères et notamment des élastomères reste un problème délicat. Dans ce travail, nous proposons une méthode d'identification qui associe la mesure de champ de déformation par analyse d'images avec l'optimisation d'un champ de contraintes adapté à l'essai. L'essai retenu est un étirage biaxial réalisé sur une éprouvette en forme de croix. L'approche proposée transforme l'inconvénient de l'hétérogénéité en avantage puisqu'il permet de réaliser l'identification simultanée sur plusieurs états de déformation : typiquement traction uniaxiale, plane et biaxiale. Le champ de contrainte est approché par la somme d'un champ homogène et d'un champ complémentaire vérifiant les conditions de bords libres et qui décroît lorsqu'on pénètre dans l'échantillon. La longueur caractéristique de la décroissance est optimisée de telle sorte que le champ approché vérifie au mieux les équations d'équilibre. En combinant l'analyse d'images avec le champ de contrainte optimisé, on identifie le potentiel hyperélastique en calculant explicitement les deux dérivées f=∂W/∂I1 et g=∂W/∂I2 où et sont les deux 1er invariants du tenseur de Cauchy droit. Enfin, un algorithme spécifique est mis en oeuvre par éléments finis pour une simulation 2D des matériaux hyperélastiques incompressibles. Cet algorithme est utilisé pour valider l'identification en comparant les résultats de la simulation et ceux de l'expérience / The identification of the polymer's behaviour and especially rubber-like materials remains a challenging task. In this work, we propose a method for the identification which combines strain field obtained by digital image analysis and the optimisation of an approximated stress field adapted to the specimen geometry. A biaxial stretching test is performed on a crosshair rubber specimen. With the proposed approach, heterogeneity of the strain field during this equi-biaxial tension test becomes an advantage. It allows the simultaneous identification of several strain states: uniaxial, biaxial and planar elongations as well as shear. The stress field is approximated by the sum of a homogeneous field and an additional field. The latter, checking the boundary conditions on the free edge is decreasing when entering the sample. The characteristic decreasing length is optimized so that the approximate field verifies the equilibrium equations. Combining image analysis with an optimized stress field, we manage the identification the hyperelastic potential by calculating explicitly the two derivatives f=∂W/∂I1 and g=∂W/∂I2 and conclude on their dependence on I1 and I2 the two first invariants of the rigth Cauchy-Green tensor. Finally, a specific finite element algorithm has been developed to similate a 2D-incompresible hyperelastic material. This algorithm is used to validation the identification potentiel by comparing simulation results and experimental data