La prédiction des déformations différées est d'une très grande importance pour l'étude de la durabilité et de l'aptitude au fonctionnement à long terme des structures en béton (ponts, enceintes de confinement de bâtiments réacteurs des centrales nucléaires, etc.). En effet, elles peuvent être à l'origine de la fissuration, de pertes de précontrainte, d'une redistribution des contraintes et même, plus rarement, de la ruine de l'ouvrage.<br /><br />L'objectif de ce travail est alors de développer des outils de calcul numérique, capable de prédire le comportement différé de structures en béton. Pour cela, un nouveau modèle hydromécanique du béton est développé, intégrant la description des phénomènes de séchage, de retrait, de fluage et de fissuration. La modélisation du retrait de dessiccation est basée sur une approche unifiée du fluage et du retrait. Les modèles de fluage propre et de fluage de dessiccation sont basés sur des mécanismes physico-chimiques plausibles, se produisant à différentes échelles d'observation de la pâte de ciment. Le modèle de fluage propre est associé à la micro-diffusion de l'eau adsorbée entre la porosité interhydrates et intrahydrates et la porosité capillaire, et au glissement des feuillets de C-S-H à l'échelle des nanopores. Le fluage de dessiccation est induit par la micro-diffusion de l'eau adsorbée à différentes échelles de porosité sous l'effet d'une sollicitation mécanique et hydrique combinée. Le retrait de dessiccation résulte, en effet, de la déformation élastique et différée du squelette solide, sous les effets de la pression capillaire et de la pression de disjonction. Le comportement mécanique du béton fissuré est modélisé en utilisant le formalisme de l'élastoplasticité endommageable orthotrope. La combinaison de ces phénomènes est effectuée dans le cadre de la mécanique des milieux poreux non saturés, en s'appuyant sur le concept des contraintes effectives.<br /><br />Ce modèle a été incorporé dans un code de calcul aux éléments finis. L'analyse du comportement différé d'éprouvettes et de structures en béton et en béton précontraint, soumises à des sollicitations hydriques et mécaniques combinées, est alors présentée.

L'objet de cette thèse est l'étude de modèles mathématiques pour les phénomènes de couplage entre l'écoulement de fluides et la déformation mécanique du sol lors de l'extraction d'hydrocarbures en milieu poreux. Dans la partie théorique, on considère deux modèles de couplage, d'une part (1) entre les déformations du sol et un écoulement linéaire compressible, et d'autre part (2) entre les déformations du sol et un écoulement diphasique non linéaire. Pour le modèle (1), on prouve l'existence et l'unicité d'une solution faible par la méthode de Galerkin. Le modèle (2) est fortement couplé et comporte une équation parabolique dégénérée; pour démontrer l'existence de solution, on considère une suite de problèmes uniformément paraboliques associés et on démontre qu'ils admettent une solution classique l'aide du théorème de point fixe de Schauder. On s'appuie ensuite sur le théorème de Fréchet-Kolmogorov pour prouver la compacité relative des suites de solutions et établir la convergence d'une sous-suite vers une solution faible du problème initial. Dans une seconde partie, on aborde l'étude numérique. On compare deux algorithmes pour les modèles de couplage. Le premier, utilisé par les ingénieurs du pétrole, est basé sur une méthode de point fixe; le second, que nous proposons et qui est plus robuste que le premier, s'appuie sur la méthode du gradient conjugué préconditionné.

Le premier chapitre présente une synthèse bibliographique sur la rigidité des sols en très petites déformations ainsi qu'une revue des principales techniques de mesure in situ et en laboratoire. Les principaux paramètres qui influencent le module maximum sont discutés. La dépendance rigidité-déformation des sols est présentée avec l'exemple de la courbe en "S". Dans le deuxième chapitre, on examine, avec la méthode des éléments finis, la stabilité d'un anneau en maçonnerie placé dans un massif élastique linéaire semi-infini et soumis à un chargement d'origine géostatique. On met ainsi en évidence un domaine de stabilité limité par deux valeurs du coefficient de pression horizontale des terres. Le troisième chapitre est consacré à l'étude du comportement en service des collecteurs enterrés. On étudie, par la méthode des éléments finis, l'impact d'une charge linéique appliquée à la surface d'un massif hyperbolique sur une conduite faiblement enterrée. On montre que le comportement non linéaire du massif a un effet très limité sur les déformations de la conduite, même lorsque le chargement approche la charge limite du sol, car les déformations du sol autour de la conduite restent presque partout inférieures à 0,01%. Le quatrième chapitre est consacré aux expérimentations in situ sur des collecteurs enterrés. On présente l'essai de vérinage interne ainsi que sa méthodologie d'interprétation. Cet essai mis en œuvre à l'intérieur des conduites permet de déterminer par calcul inverse et moyennant certaines hypothèses, le module en petites déformation du sol encaissant et le module de la conduite.

Ce travail a pour but la mise au point d'une méthodologie de formage virtuel qui consiste à simuler numériquement tous procédés de mise en forme de pièces minces ou massives avec la prise en compte d'un couplage thermo-mécanique-endommagement. Ceci permettra de prévoir l'apparition de l'endommagement ductile au cours du formage des pièces. On pourra ainsi agir sur les paramètres technologiques pertinents du procédé afin soit de retarder cet endommagement pour obtenir des pièces saines, soit au contraire de favoriser celui-ci afin de simuler des opérations de coupe. Une formulation théorique générale du couplage "fort" thermo-élasto-visco-plastique-endommagement est introduite, prenant en compte une loi d'évolution de l'endommagement ductile isotrope. Sur le plan numérique, un soin particulier a été apporté à l'intégration locale des équations d'évolution couplées. Pour résoudre le système global couplé, deux schémas ont été utilisés : (1) un schéma itératif statique implicite associé au calcul analytique de la matrice tangente consistante, (2) un schéma dynamique explicite associé à un contrôle automatique du pas de temps. Une procédure de remaillage adaptatif en grandes déformations a été utilisée en connexion avec le code de calcul par éléments finis ABAQUS. Après une validation sur des essais de traction, une large gamme de procédés, dont certains sont issus de l'industrie, tels que le découpage de tôle, la compression de cylindres, la coupe orthogonale par enlèvement de copeaux,... est présentée et comparée avec des essais expérimentaux quand cela est possible. La capacité de cette méthodologie de formage virtuel à prédire correctement l'initiation et la propagation de l'endommagement en mise en forme est clairement démontrée.

Ce travail porte sur des essais de traction réalisés sur des éprouvettes d'aluminium lithium entaillées, avec des observations in-situ du champ de déformations.La réduction de consommation de carburant est actuellement une priorité pour l'aviation. Les alliages d' Al-Li sont d'excellents candidats pour réduire le poids des avions puisqu'ils combinent de très bonnes propriétés mécaniques avec une densité inférieure à celle des alliages conventionnels. Cependant, les propriétés mécaniques de ces matériaux sont fortement anisotropes, et il est essentiel de contrôler ce phénomène afin d'utiliser le matériel à bon escient. À cette fin, il est nécessaire d'avoir le maximum d'informations sur le matériau tant sur ses propriétés microstructurales (taille et forme des grains...) que mécaniques (élasticité, module de Young, résistance à la fissuration ...). Mais la connaissance de ces propriétés mécaniques globales n'est pas suffisante pour maximiser la résistance à la fissuration (critique dans l'application aéronautique). Dans ce contexte, ce travail vise à quantifier l'influence de la microstructure locale (orientation et la taille des grains) sur la fissuration.La propagation des fissures dans des tôles d'épaisseur 2 mm a été analysée. Les essais de traction ont été effectués sur des échantillons entaillés avec observation in-situ du champ de déplacement.Ainsi, la méthode de la grille a été adaptée pour déterminer de grandes déformations dans les grains. Le champ de déplacement a été caractérisée par la déformation d'une grille croisée collée (pas 30 um), et les paramètres affectant la qualité des résultats (résolution et résolution spatiale) ont été optimisés. Afin de compenser le mouvement du corps rigide, un nouveau système automatique d'acquisition d'image conduisant à de petites déformations entre chaque image a été développé. Une nouvelle méthode d'extraction de phase a été instauré, permettant un déroulage temporel de la phase.

Ce travail de recherche propose un modèle numérique pour simuler les procédés par infusion de résine en utilisant la méthode des éléments finis. Ce modèle permet de représenter l'écoulement d'une résine liquide dans des préformes poreuses subissant de grandes déformations. Dans cette étude, une modélisation macroscopique est utilisée. Au niveau du procédé, une zone de résine liquide est déposée sur les préformes. Ces dernières étant considérées comme un milieu poreux. Les équations de Stokes et de Darcy sont utilisées pour modéliser l'écoulement de la résine respectivement dans le drainant et dans les préformes. L'originalité du modèle réside dans le fait qu'un seul maillage est utilisé pour les deux milieux. La discrétisation est réalisée avec des éléments mixtes : dans Stokes, des éléments P1+/P1 sont utilisés et dans Darcy, des éléments P1/P1 stabilisés avec une formulation multi-échelle sont employés. Des fonctions distances signées sont utilisées pour représenter l'interface entre Stokes-Darcy et pour représenter le front de résine. Concernant la déformation des préformes, une formulation Lagrangienne réactualisée est utilisée. Dans cette formulation Lagrangienne, le comportement des préformes humides est représenté à l'aide du modèle de Terzaghi dans lequel les préformes sèches ont un comportement élastique non-linéaire. La perméabilité est reliée à la porosité via la relation de Carman-Kozeny. Celle-ci est déterminée à partir de l'équation de conservation de la masse. Ce modèle a été implémenté dans ZéBuLoN. Plusieurs simulations numériques d'infusion de résine sont présentées à la fin de ce manuscrit.

Le J. Boudinar est situé à l'intersection de plusieurs directions rhegmatiques africaines. Il représente en effet un segment d'un des couloirs tectoniques atlasiques (NE-SW) de Tunisie centrale, affecté par l'axe N-S et ses satellites (faille Essouda-Bagra) et les failles de Kasserine (au NE) et de Mezzouna (au SW) de direction NW-SE. L'étude stratigraphique et structurale de ce secteur s'appuie sur une cartographie au 20 000 eme. STRATIGRAPHIE - SEDIMENTOLOGIE : L'examen banc par banc de nombreuses coupes de terrain a permis de préciser l'âge de certaines unités lithostratigraphiques et de caractériser les variations d'épaisseur de chaque Formation du Jurassique à l'Actuel. Ainsi , la base de la Formation Aleg est datée du Turonien terminal, et la Formation Abiod du Campanien à Maastrichtien inférieur. 7 mégaséquences de comblement ont été définies du Tithonique à l'Actuel. Elles sont délimitées par des discontinuités, dont 3 majeures situées à l'Albien inférieur-moyen (?), au Turonien terminal et au Maastrichtien à l'Eocène moyen. Elles sont liées soit à des émersions, soit à une submersion des plates-formes antérieurement exondées, soit encore à une érosion sous-marine. Toutes sont accompagnées d'une induration en milieu aérien ou sous une faible tranche d'eau. HISTOIRE STRUCTURALE : Des variations dans le temps du champ de contraintes appliqué à un bâti structuré aboutissent à une succession et une alternance de stades d'extension et de compression accompagnés à chaque fois par une halocinèse. 5 étapes de déformations ont été distinguées au J. Boudinar : 1 - Une extension N-S caractérise le Crétacé inférieur. Elle provoque le jeu de fractures à toutes échelles selon 4 directions (N 65-85; N100-110; N 130-140; N 165-170). Le couloir tectonique atlasique induit l'effondrement d'un graben kilométrique orienté N 65-85. 2 - Une extension à dominante NE-SW caractérise le domaine de l'Albien supérieur au Turonien. Elle provoque la formation d'aires plus ou moins subsidentes et de blocs basculés orientés N 120-140, N 45-50 et N-S. 3 - L'halocinèse domine du Sénonien à l 'Eocène supérieur. Elle se traduit par la mise en place de rides salifères orientées NE-SW et N-S donnant lieu à une intrusion salifère à leur intersection. Elle pourrait représenter un écho des compressions at1asiques. 4 - Une extension multidirectionnelle caractérise le Miocène moyen, et mobilise des failles orientées N 60 et N 140. 5 - A partir du Miocène supérieur et jusqu'à l'Actuel, le domaine est le Siège de compressions atlasiques provoquant la formation d'un anticlinal à coeur évaporitique et d'un synclinal superposé au collapse salifère Sénonien-Eocène supérieur. Les extrémités du pli sont crochonnées lors du jeu senestre des couloirs tectoniques N-S. Un rejeu dextre est envisagé pour expliquer la rotation horaire de la branche centrale du pli. Le J. Boudinar résulte d'une inversion structurale positive au droit d'un couloir tectonique atlasique (NE-SW). SUBSIDENCE: L'analyse de la courbe des épaisseurs cumulées permet d'individualiser des anomalies (=crises) dans l'évolution de la subsidence. Une corrélation avec l'ensemble la Tunisie montre qu'il s'agit de phénomènes régionaux liés aux évolutions générales africaine, tethysienne-atlantique, et méditerranéenne.

La tectonique des plaques repose sur le postulat que les plaques sont rigides et que l'essentiel des déformations se concentre sur leurs frontières, généralement étroites en domaine océanique. Plusieurs plaques océaniques semblent y déroger. La plaque Indo-Australienne montre par exemple une forte sismicité en dehors de ses frontières, qui traduit une déformation intraplaque intense. Les mécaniques au foyer expriment une déformation compressive dans la partie Est de la plaque, et une déformation extensive sur le flanc Est de la dorsale centrale Indienne et sur le flanc Nord de la dorsale Sud Est Indienne. L'inversion de données cinématiques, qui supposent des plaques rigides, abouti au concept de plaques composites composées d'éléments rigides séparés par des limites larges et diffuses. En utilisant l'exemple de la plaque Indo-Australienne, ce travail cherche à savoir comment se comporte la lithosphère océanique quand elle est soumise à des contraintes horizontales et le rôle joué par les structures préexistantes comme les plateaux sous-marins ou les zones de fracture fossiles sur la déformation. La première partie de ce travail s'est intéressée à la sismicité de faible magnitude détectée par un réseau temporaire de 3 hydrophones déployés dans l'océan Indien Austral. L'objectif était de mieux identifier et localiser les structures actives des zones déformées. En raison de la géométrie du réseau, ces données apportent peu d'informations nouvelles mais confirment le rôle des zones de fracture fossiles dans la déformation des domaines en extension. Par ailleurs, la découverte, sur des zones de fracture actives des dorsales Est Pacifique et Juan de Fuca, de précurseurs sismiques nous a incités à en rechercher sur nos données. Un précurseur sismique est défini comme un événement se produisant dans un rayon de 100 km autour d'un événement majeur, quelques heures à quelques jours avant celui-ci. L'analyse d'une sélection d'évènements détectés par le réseau permanent du traité international de non-prolifération nucléaire (OTICE), dont la localisation est très précise, n'a pas permis de reconnaître de tels précurseurs. En revanche, notre réseau temporaire détecte de nombreux évènements qui pourraient y être assimilés, si leur localisation, basée sur seulement 3 hydrophones, pouvait être améliorée. La combinaison des données de ces 2 réseaux serait nécessaire pour progresser sur ces deux questions. La seconde partie de ce mémoire est consacrée à la modélisation de la déformation de la plaque Indo-Australienne, afin notamment de préciser le rôle de structures préexistantes, reliefs ou zones de fracture fossiles remobilisées, comme en atteste la sismicité. Pour aborder cette question, nous avons testé deux outils " ADELI " (Hassani et al., 1997) et " SHELLS " (Bird, 1989; Bird et Kong, 1994; Kong et Bird, 1995; Bird, 1999). Le second code numérique s'est avéré le plus adapté. Ce modèle décompose les plaques en éléments finis triangulaires sur la sphère, tient compte de l'âge de la lithosphère considérée, et peut prendre en compte des zones de faiblesse (zones de fracture) et des reliefs. A partir de conditions aux limites en vitesse et d'hypothèses sur la rhéologie de la lithosphère, le modèle calcule l'état des contraintes et des déformations de la plaque considérée. La plaque Indo-Australienne présentant des zones déformées disjointes, nous avons d'abord modélisé chacune d'elles séparément : la zone en extension, entre les plaques Inde et Capricorne est découpée par un réseau de zones de fracture oblique par rapport à l'orientation générale Nord-Sud de la dorsale Centrale Indienne et comprend le relief des Chagos ; la zone en extension entre les plaques Capricorne et Australie est découpée par un réseau de zones de fracture perpendiculaires à l'orientation générale Ouest-Nord-Ouest Est-Sud-Est de la dorsale Sud-Est Indienne et comprend le plateau océanique de Saint-Paul et Amsterdam ; enfin, la zone en compression entre les plaques Inde, Capricorne et Australie, est caractérisée par une lithosphère ancienne découpées par des zones de fracture Nord-Sud dans les bassins de Wharton et Indien Central et la ride du 90ème Est. Les conditions aux limites de ces modèles sont issues d'une approche cinématique de plaque rigide. Pour intégrer les limites de plaques diffuses sur les bords des maillages, nous avons utilisé la notion de bord libre, cela nous a permis de reproduire les déformations dans ces trois régions sans imposer de limites particulières entre les plaques. Enfin, nous avons réuni les trois maillages afin de modéliser la plaque Indo-Australienne dans son ensemble. Les modélisations effectuées ont mis en évidence l'importance des structures préexistantes comme le banc des Chagos qui concentre la déformation au niveau de la limite diffuse entre les plaques Inde et Capricorne où des séismes de magnitude 8 se sont produits. Les zones de fracture jouent par leur réactivation un rôle essentiel, notamment à l'Est de la ride du 90ème Est, où elles provoquent des comportements différents entre le bassin Indien Central où la déformation est marquée essentiellement par de la compression dans une direction Nord-Sud et le bassin de Wharton, où du décrochement sénestre se manifeste en plus de la compression. Nos modèles prédisent des taux d'extension et de raccourcissement proches d'études précédentes avec un éloignement de l'Inde par rapport à la plaque Capricorne d'environ 2 mm/an dans une direction Nord-Sud, et un raccourcissement de l'ordre de 11,5 mm/an entre l'Australie et l'Inde dans une direction Nord-Ouest Sud-Est. Au premier ordre, le concept de plaque composite semble valide pour la plaque Indo-Australienne, cependant la répartition des zones diffuses s'avère plus complexe que prévue.

On étudie les surfaces plates à singularités coniques, leur géométrie, leur espaces des déformations et leur paramétrisation. La surface de base est la sphère à trois trous (pantalon). On trouve trois ensembles de paramètres pour le pantalon plat à un point singulier conique et on décrit son espace des déformations. On introduit un flot que l'on appelle flot de Fenchel-Nielsen sur un espace des déformations. On étudie l'injectivité de ce flot en examinant la variation des fonctions longueur de segments géodésiques ou de géodésiques simples fermées le long de ce flot. On étudie également la paramétrisation d'une surface plate à singularités coniques utilisant des longueurs des segments géodésiques joignant des points singuliers ou un point singulier à une composante du bord. A la fin du texte, trois annexes apportent des discussions supplémentaires.

Les comportements mécaniques des films minces polycristallins sont encore mal compris à l'échelle sub-micronique. En particulier des hétérogénéités locales de déformation importantes sont attendues, mais elles restent difficile à quantifier expérimentalement. Les nouvelles possibilités offertes par les micro-faisceaux synchrotron de rayons X ont donc été utilisées dans ce travail pour éclairer cette problèmatique. Une réflexion de Bragg provenant d'un grain unique sub-micronique a été acquise avec une très bonne résolution dans l'espace réciproque en trois dimensions lors d'un cycle thermique. Les propriétés de cohérence du faisceau ont été utilisées pour reconstruire à trois dimensions une composante du champ de déplacement intra-grain avec une résolution d'une vingtaine de nanomètres dans les trois directions. Cette technique est basée sur des algorithmes de reconstruction de phase qui néanmoins connaissent des stagnations dans le cas des échantillons fortement déformés. Une méthodologie basée sur la connaissance de la forme du grain a donc été développée pour contourner ces difficultés. Des analyses complémentaires de diffraction X de laboratoire et de microdiffraction monochromatique ont également mis en évidence des hétérogénéités importantes de déformation entre les différents grains.