Ce travail concerne la mise en oeuvre expérimentale de la méthode des champs virtuels (MCV) pour identifier des paramètres pilotant une loi de comportement élastoplastique. La méthode a déjà été validée sur des données simulées, mais plusieurs aspects restaient à analyser pour une application expérimentale. La méthode d'identification est basée sur la minimisation d'une fonction coût traduisant l'écart à l'équilibre de la structure testée. Cette fonction coût est construite à partir du principe des travaux virtuels (pour plusieurs champs virtuels choisis arbitrairement), des champs cinématiques mesurés à la surface de l'éprouvette et de la force de traction appliquée. Dans un premier temps, un essai de traction sur une éprouvette plane à section variable, conduisant à un état de contrainte quasiment uniaxial, a été mis en place. Cette essai a conduit à l'identification des six paramètres d'une loi élastoplastique isotrope. Les valeurs des paramètres identifiés coïncident avec les valeurs de référence identifiées lors d'essais de traction classiques. La stabilité de l'identification vis à vis des paramètres influant a été examinée. Dans un deuxième temps, des essais mécaniques générant des champs de contraintes multiaxiales ont été conduit sur des éprouvettes planes à entailles circulaires. Les paramètres identifiés à partir de ces essais se sont révélés d'avantage éloignés des valeurs de référence. Plusieurs hypothèses ont été évoquées pour expliquer ces écarts, comme une construction non optimale de la fonction coût ou des erreurs lors du choix du modèle décrivant le comportement du matériau.

L'allègement de la masse des structures est au cœur des développements des industries du transport. Pour y parvenir, les pièces intègrent de plus en plus de fonctions, leurs géométries deviennent alors très complexes. Pour réaliser ces produits, les industriels ont recours à des stratégies d'optimisation utilisant des méta-modèles numériques prenant en compte l'enchaînement des étapes de conception et de fabrication des pièces. Cependant les temps de simulation sont encore très élevés. Ces travaux se sont intéressés à l'amélioration des protocoles de simulation numérique de procédés de mise en forme d'emboutis profonds. Dans un premier temps, ils ont permis de développer des méthodes numériques pour diminuer les temps de simulation du procédé de formage incrémental. Un schéma de résolution en dynamique explicite est utilisé avec une vitesse de l'outil adaptée. Cette méthode à permis de réduire significativement le temps CPU, tout en conservant une bonne prédiction des géométries et des épaisseurs. Les travaux se sont également intéressés à l'utilisation de la théorie de la déformation incrémentale pour diminuer le temps de résolution du calcul élasto-plastique. De bons résultats sont observés mais le gain de temps est faible (4.5%). Une approche simplifiée du contact a également été développée. De bons résultats sont obtenus avec un gain de temps d'un facteur 2.13 en comparaison de celui obtenu avec une vitesse adaptée. Dans un deuxième temps, une étude sur la réalisation d'emboutis profonds par formage superplastique a été menée de manière à pouvoir à terme enchaîner les deux procédés. Elle a porté sur la détermination d'une loi de pression optimale. Deux algorithmes ont été développés en fonction du protocole pour entrer la loi de pression dans les machines industrielles. Les algorithmes développés permettent un bon contrôle du domaine superplastique. Un temps de calcul réduit d'un facteur 3 est observé, en comparaison avec des méthodes issues de la littérature.

L'étude de sûreté d'un stockage souterrain de déchets radioactifs passe, du point de vue de la géomécanique, par la caractérisation de la zone endommagée induite par le creusement du puits d'accès aux cavités de stockage. L'objet de cette thèse est donc de simuler la réponse hydro-mécanique des argilites de l'Est au creusement du puits d'accès au laboratoire souterrain de Meuse/Haute-Marne, en vue d'une comparaison avec les mesures qui seront acquise in situ en 2005 (expérimentation REP).<br/> Nous avons d'abord analysé l'ensemble des essais mécaniques qui ont été réalisés entre 1995 et 2001. Nous observons que chaque série d'essais a sa propre cohérence, que ce soit en termes de paramètres élastiques, de résistance à la rupture ou de capacité au fluage. Mais il existe de fortes différences entre les séries. Ces différences, qui sont dues au protocoles expérimentaux et non au matériau, nous ont permis de mettre en évidence trois résultats importants : les phases de resaturation imposées par plusieurs laboratoires ont dégradé les propriétés mécaniques des argilites de l'Est, l'existence d'une contrainte effective au sens de Biot n'est pas évidente pour ces argilites, et leur résistance à la rupture augmente nettement avec le séchage.<br/>D'après ces essais, nous avons distingué trois mécanismes de déformation irréversible : compaction, pré-rupture, rupture. Nous les avons décrits chacun par une loi de comportement élasto-plastique simple, fondée sur le critère de Mohr-Coulomb ou de Drücker-Prager, avec un écrouissage linéaire. On obtient alors une loi de comportement complète pour les argilites de l'Est en considérant les trois mécanismes simultanément (plasticité multi-critère).<br/>Enfin, les différentes modélisations que nous avons réalisées permettent de tirer les conclusions suivantes :la prise en compte des déformations irréversibles avant-rupture (compaction et pré-rupture) a pour effet de réduire considérablement l'extension de la zone en rupture qui apparaît autour du puits (par rapport au cas où le comportement est élastique jusqu'à la rupture) ; la diffusion des pressions interstitielles modifie peu la réponse mécanique du massif (en particulier, les déformations plastiques créées lors du déchargement mécanique n'évoluent pas au cours du drainage du massif) ; l'anisotropie des contraintes naturelles horizontales, dont par ailleurs nous confirmons l'existence, influence fortement la zone endommagée autour du puits (polarisation en azimut) ; enfin nous montrons qu'il n'est pas indispensable de simuler le détail des opérations de creusement du puits.

Ce travail de recherche, qui s'inscrit dans le cadre d'un projet Européen, se propose d'étudier le comportement des tôles métalliques lors des opérations de mise en forme. Les applications visées sont la prédiction, au moyen de la simulation numérique du procédé, des défauts majeurs tels que le retour élastique et la localisation des déformations. A cette fin, une modélisation phénoménologique du comportement a été retenue. Un modèle de comportement élasto-plastique a été considéré, dans lequel deux modèles d'écrouissage sont introduits pour reproduire l'anisotropie induite lors de trajets complexes de chargement. En particulier, le modèle microstructural de Teodosiu-Hu, basé sur l'évolution des structures de dislocations, a été sélectionné afin de bien prendre en compte les phénomènes de transition liés aux changements de trajets de déformation. Ce modèle a été enrichi par la suite, en le couplant au modèle d'endommagement de Lemaitre-Chaboche, afin de reproduire la chute de la contrainte d'écoulement à l'approche de la rupture du matériau. De plus, pour modéliser la limite de ductilité des tôles métalliques, le critère de localisation par bandes de cisaillement de Rice a été introduit, dans le cadre d'une formulation complètement tridimensionnelle en grandes déformations. La modélisation développée tient ainsi compte de nombreux phénomènes tels que l'anisotropie plastique initiale, celle induite par écrouissage, l'adoucissement généré par l'endommagement et enfin la localisation des déformations. Pour permettre l'exploitation de cette modélisation dans le cadre de la simulation numérique des procédés de mise en forme, une implantation numérique robuste et efficace a été réalisée dans un code de calcul par éléments finis. Différents schémas d'intégration ont été testés et comparés afin de retenir le meilleur compromis entre précision et temps de calcul. A des fins de validation, cette modélisation a d'abord été appliquée à la simulation d'essais rhéologiques à trajets directs et séquentiels. En comparaison avec les courbes expérimentales, les prédictions obtenues ont bien montré toute la capacité du modèle élasto-plastique à écrouissage microstructural à reproduire les phénomènes de transition observés sur des trajets séquentiels. Par la suite, l'application de cette modélisation à l'étude du retour élastique a mis en évidence les effets de plusieurs paramètres (rhéologiques, numériques ou liés au procédé). En particulier, il a été montré, à travers une application spécifiquement déterminée, que l'influence d'un modèle d'écrouissage avancé, tel que le modèle microstructural, pouvait se révéler importante dans certaines applications. Enfin, une analyse de la localisation des déformations a été menée en utilisant le critère de Rice conjointement au modèle élasto-plastique couplé à l'endommagement. Cette étude a montré la capacité de cette approche à prédire les courbes limites de formage, ainsi que l'orientation des bandes de localisation, sur des trajets de chargement directs et séquentiels.

Ce travail de thèse présente la mise au point d'une méthodologie de mesure de champs de déformations planes à l'échelle micrométrique. La mesure cinématique est faite par une méthode de grille qui analyse les déformations d'un réseau périodique attaché à la surface à étudier. Les grilles ont été réalisées par photolithographie interférentielle directe. Pour l'observation et la numérisation des réseaux, le choix s'est porté sur la microscopie interférométrique en lumière blanche qui donne accès à la topographie 3D de la surface d'intérêt. Une étude de bruit à permis d'optimiser la réalisation des réseaux pour l'application mécanique auxquels ils étaient destinés. Par ailleurs, l'application de déplacements de corps rigide a permis de mettre en évidence un problème de corrélation spatiale de ce bruit qui ne pouvait être ignoré. Ainsi, et en apportant un soin particulier sur la réalisation expérimentale des essais mécaniques, la résolution en déformation obtenue est de l'ordre de 1 à 2×10−3 pour une résolution spatiale d'environ 20 μm. Une première application concernant l'étude des déformations élasto-plastiques d'un acier ferritique a permis de valider la présente méthode, de l'étendre à l'étude de régions d'intérêt de plus grandes tailles et d'en souligner les points forts mais aussi les limites.

Lors de la conception des structures du génie civil et du génie mécanique il est très important de pouvoir prévoir comportement près de leur état ultime. Or, les matériaux utilisés sont souvent fortement hétérogènes, comme par exemple le béton et les métaux composites, et sont en conséquence difficiles à modéliser en cas de sollicitations extrêmes. Dans ce travail nous proposons une nouvelle stratégie multi échelles adaptée pour ce type des problèmes. Selon l'approche on couple un modèle basé sur la méthode des éléments finis (MEF) à l'échelle de la structure (macro) avec un modèle très fin à l'échelle de la microstructure (micro), lui aussi basé sur la MEF. Cela nous permet de modéliser l'échelle micro beaucoup plus fidèlement que ce qui est possible avec les approches phénoménologiques ou avec des méthodes d'homogénéisation analytique. L'inconvénient principal d'une telle approche étant en général le coût de calcul très élévé, nous considérons trois aspects rendant possible une application aux situations réelles. Premièrement, la formulation variationnelle permet de n'utiliser le modèle fin que dans les parties restreintes de la structure, là où les échelles sont vraiment fortement couplées, tout en gardant un modèle homogénéisé ailleurs. Deuxièment, nous adaptons la MEF sur l'échelle micro en employant l'approche dite structurée, avec laquelle on peut réduire le nombre de degrés de liberté d'une manière importante en gardant la même précision. Le troisième aspect porte sur l'implantation de la méthode dans un code de calcul, adaptée pour des machines parallèles, ce qui nous a permis de lancer les analyses d'une taille comparable à celle des situations réelles. Enfin, nous montrons de quelle manière cette approche peut être utilisée dans l'optimisation de forme des inclusions d'un matériau composite.

Ce travail de thèse présente la mise au point d'une méthodologie de mesure de champs de déformations planes à l'échelle micrométrique. La mesure cinématique est faite par une méthode de grille qui analyse les déformations d'un réseau périodique attaché à la surface à étudier. Les grilles ont été réalisées par photolithographie interférentielle directe. Pour l'observation et la numérisation des réseaux, le choix s'est porté sur la microscopie interférométrique en lumière blanche qui donne accès à la topographie 3D de la surface d'intérêt. Une étude de bruit à permis d'optimiser la réalisation des réseaux pour l'application mécanique auxquels ils étaient destinés. Par ailleurs, l'application de déplacements de corps rigide a permis de mettre en évidence un problème de corrélation spatiale de ce bruit qui ne pouvait être ignoré. Ainsi, et en apportant un soin particulier sur la réalisation expérimentale des essais mécaniques, la résolution en déformation obtenue est de l'ordre de 1 à 2×10-3 pour une résolution spatiale d'environ 20 µm. Une première application concernant l'étude des déformations élasto-plastiques d'un acier ferritique a permis de valider la présente méthode, de l'étendre à l'étude de régions d'intérêt de plus grandes tailles et d'en souligner les points forts mais aussi les limites.

La simulation numérique du comportement et l'analyse de stabilité des ouvrages de géotechnique renforcés par inclusions linéaires "rigides" ou "souples", c'est-à-dire capables ou non de reprendre des efforts de flexion et cisaillement, reste un problème difficile, en raison de la très forte hétérogénéité du sol renforcé combinée au grand nombre d'inclusions employées, ainsi qu'à leurs dimensions caractéristiques très petites vis à vis de celles de l'ouvrage dans son ensemble. Une approche alternative, fondée sur une modélisation dite "multiphasique", a été récemment développée, qui permet de s'affranchir des difficultés précédentes et d'aboutir ainsi à des méthodes de dimensionnement accessibles à l'ingénieur. On s'intéresse dans le cadre de ce travail aux développements théorique et numérique de cette approche en prenant en compte non seulement les effets de flexion et de cisaillement des inclusions, mais également les différents types d'interaction entre ces dernières et le sol environnant. La première partie du mémoire est consacrée à la construction des approches du calcul à la rupture et à la mise en œuvre numérique élasto-plastique du modèle multiphasique des ouvrages renforcés par inclusions "souples", qui est alors appliqué à l'analyse de la stabilité d'ouvrage en "terre armée". L'extension de l'approche du calcul à la rupture et du code éléments finis en élastoplasticité, développés dans la première partie, à la prise en compte des effets de flexion et de cisaillement fait l'objet de la seconde partie du mémoire. La version la plus complète du modèle multiphasique élaborée est appliquée à la simulation du comportement d'une fondation renforcée par inclusions "rigides", soumise à un séisme induisant un chargement latéral. Partant de l'idée d'optimiser le schéma de renforcement de l'ouvrage traité à la deuxième partie, en inclinant les inclusions, un modèle "triphasique" est formulé et développé. Ce modèle vise à décrire le comportement macroscopique d'un sol renforcé par deux familles d'inclusions ayant des orientations différentes.

Ces travaux de thèse portent sur l'étude expérimentale des propriétés physiques et de la biomécanique des agrégats cellulaires embryonnaires. Le but de cette thèse était d'une part de mieux comprendre l'origine biologique de la viscosité et de la tension de surface tissulaire, d'autre part d'étudier quantitativement en détail l'élasticité cellulaire par des nouvelles mesures de rhéologie en cisaillement. Un premier chapitre concerne les mesures de tension de surface tissulaire par la méthode de compression et de viscosité tissulaire par analyse de la cinétique de fusion de deux agrégats en faisant varier comme paramètre principal la contractilité cellulaire que certains suspectent comme étant la principale origine biologique de ces paramètres. Nous utilisons le formalisme du DITH (Haris, 1976: Differential Interfacial Tension Hypothesis) pour interpréter les données. Le deuxième chapitre concerne les mesures rhéologiques en cisaillement à l'aide d'un rhéomètre commercial plan-plan sur plusieurs centaines ou milliers d'agrégats cisaillés ensembles. Nous montrons que les cellules deviennent moins rigides pour une déformation minimale d'environ 4%, mais sur l'échelle de l'heure les cellules sont capables de se rigidifier à nouveau. Ces expériences sont analysées à l'aide d'un modèle de ressorts qui cassent sous contrainte puis se ressoudent à contrainte nulle.

Nous présentons deux approches différentes pour décrire le couplage hydromécanique des géomatériaux. Dans une approche de type phénoménologique nous traitons le milieu poreux comme un milieu continu équivalent dont les interactions entre la phase fluide et le squelette solide constituent le couplage du mélange à l'échelle macroscopique. En caractérisant le comportement de chaque phase nous arrivons à décrire le comportement couplé du milieu couplé saturé.Nous utilisons cette approche pour modéliser des essais expérimentaux faits sur un cylindre creux pour une roche argileuse (argile de Boom). Les résultats expérimentaux montrent de façon claire que le comportement de cette roche est fortement anisotrope. Nous avons choisi de modéliser ces essais en utilisant une lois de comportement élasto-plastique pour laquelle la partie élastique est transversalement isotrope.Le problème aux conditions aux limites étudié met en évidence des déformations localisées autour du forage intérieur. Afin de décrire de façon objective le développement de ces bandes de cisaillement nous avons considéré un milieu continu local de type second gradient qui permet d'introduire une longueur interne. De ce fait nous avons pu étudier le problème d'unicité en montrant qu'un changement de la discrétisation temporelle du problème aux limites peut conduire à des solutions différentes.Dans la deuxième approche étudiée nous caractérisons la microstructure du matériau avec des grains et un réseau de canaux pour la phase fluide. À l'aide d'un processus numérique d'homogénéisation nous arrivons à calculer numériquement la contrainte du mélange et le flux massique. Cette méthode d'homogénéisation numérique a été implémentée dans un code aux éléments finis afin d'obtenir des résultats macro. Une validation de l'implentation est proposée pour des calculs en mecanique pure et en hydromécanique.