Le Bouclier Baltique montre dans sa partie orientale deux ceintures de haut degré de métamorphisme (granulites de Laponie et amphibolites de la Tana), contigües et .1 situées entre deux domaines archéens (cratons d'Inari-Kola et de Laponie méridionale-Karélie). La ceinture de la Tana peut être considérée comme un ensemble volcanique (séries fractionnées tholéiitiques et calco-alcalines) et plutonique (suite anorthositique notamment). La ceinture des Granulites est interprétée comme une série méta-sédimentaire à caractère "flysch" (suite khondalitique), contenant des niveaux de métavolcanites, recoupée par un cortège d'intrusions variées d'affinité tholéiitique et calco-alcaline (complexe charnockitique). Ces deux ceintures d'âge protérozoïque inférieur (r.9 - 2.0 Ga), constituent l'axe de la zone mobile bélomoride. Les gneiss archéens (> 2.5 Ga) portent les traces d'une intense réactivation au Protérozo1"que inférieur, témoignant de leur participation à la tectogénèse bélomorienne. Enfin, le caractère épicontinental des formations volcano-sédimentaires du jatulien, leur disposition symétrique accusée par rapport aux Bélomorides, leur âge (r.9 - 2.4 Ga) et leur implication dans la tectonique karélienne montrent leur appartenance au même édifice géotectonique. Nous proposons de ré- interpréter les Bélomorides comme une chaîne par collision (cycle de Wilson), comprenant un bassin géosynclinal (ceintures de la Tana et des Granulites), une plateforme . (Jatulien) et un socle archéen ré-activé (gneiss bélomorides). Cette chaîne, la plus ancienne dans laquelle une suture ait pu être localisée est un exemple de choix pour discuter de la signification géodynamique des ceintures granulitiques, en ce qu'elle constitue probablement un jalon entre les orogènes ensialiques et les marges actives.

Ce travail traite d'une formulation pour l'étude du couplage symétrique entre éléments finis et éléments de frontière en mécanique. La formulation ainsi décrite, est basée sur l'adaptation des méthodes intégrales variationnelles de Galerkin entièrement régularisées en élasticité tridimensionnelle et s'appuie sur la construction de la matrice de rigidité d'un macro-élément représentant un domaine élastostatique tridimensionnel sans forces de volume. La procédure de couplage mise en oeuvre présente l'avantage d'être symétrique et a été implantée avec succès dans les environnements éléments finis de CAST3M. Dans le cadre du traitement des double intégrales de surface, un nouveau schéma d'intégration singulière respectant parfaitement la symétrie des formulations intégrales variationnelles mises en jeu, a été développé et présenté dans le cadre de ce travail. L'application directe de la méthode des éléments de frontière traditionnelle dite des collocations pour modéliser les solides fissurés conduit à la dégénéréscence de la formulation numérique. Une des alternatives consiste à appliquer l'approche variationnelle. Une formulation intégrale variationnelle adaptée à l'étude d'un domaine élastique tridimensionnel présentant une fissure interne non débouchante a été développée dans le cadre de ce travail. Le développement et l'implantation de la matrice de rigidité du macro-élément fissuré en termes de degrés de liberté en déplacement sur la frontière externe du domaine et en saut de déplacement sur la surface propre à la fissure, nous ont permis de procéder à une validation numérique portant sur un barreau cylindrique présentant une fissure circulaire, centrée, en mode I et non débouchante. Nous proposons également une formulation intégrale variationnelle pour l'étude d'un domaine semi-infini à surface partiellement libre et contenant une zone de non linéarités. L'utilisation des solutions fondamentales de Mindlin permettent, dans le cas de la modélisation d'un espace semi-infini, de réduire les supports géométriques à discrétiser dans le domaine éléments de frontière. Une perspective intéressante serait son implantation numérique afin de réaliser une application réelle sur un problème de géotechnique.

A la suite d'une synthèse bibliographique sommaire sur le rôle de la vibration dans le compactage des matériaux granulaires, on examine les bases théoriques relatives à la schématisation du comportement de l'ensemble sol-vibrateur à l'aide des modèles à paramètres concentrés du type masse-ressort-amortisseur ; l'analyse du rôle des caractéristiques mécaniques du matériau dans le comportement vibratoire de l'ensemble sol-rouleau permet d'interpréter plusieurs constatations expérimentales antérieures (essais sur le Vibrex). Cette approche met par ailleurs en évidence l'importance de la variabilité de l'aire de contact sol-rouleau. Les variations de l'aire de contact en fonction des caractéristiques mécaniques du matériau et de la charge statique appliquée fait l'objet d'une étude expérimentale (essais en fosse) et théorique au cours de laquelle la recherche de modèles aux éléments finis élasto-plastiques pouvant expliquer les résultats expérimentaux permet de mieux connaître le rôle des principaux paramètres dont dépend le phénomène réel. La dernière partie est consacrée à l'étude, par la méthode des Eléments Finis, des sollicitations dynamiques induites dans le massif de sol (supposé élastique ou visco-élastique). Afin de simuler les conditions aux limites réelles, on est amené à appliquer la méthode des "bords à amortisseurs visqueux" proposée par Lysmer et Kuhlemeyer.

On développe un modèle viscoélastique non linéaire en grandes déformations, appelé modèle pseudo-linéaire. Ce modèle est basé sur un choix particulier d'une mesure de déformation permettant d'écrire le potentiel d'énergie libre sous forme quadratique et la loi de comportement sous forme d'une convolution simple. Afin de valider ce modèle, on compare le potentiel pseudo-linéaire à un potentiel d'énergie libre écrit sous la forme générale d'un développement de Fréchet d'ordre 4. Les fonctions de relaxation de ce potentiel sont écrites sous forme d'exponentielles décroissantes à un temps caractéristique par décade. Cette comparaison montre une bonne concordance entre les deux descriptions. La formulation proposée permet de calculer simplement la dissipation dans le matériau. Une illustration est donnée dans le cas d'oscillations harmoniques. On développe alors une formulation variationnelle associée au modèle pseudo-linéaire basée sur la minimisation de l'énergie totale. On effectue une approximation linéaire en temps et on propose une procédure numérique de résolution. On construit alors un code de calcul par éléments finis et on le valide dans le cas de comportements simples en comparant les solutions analytiques aux résultats numériques.

On propose une modélisation numérique des plaques multicouches en matériaux composites lorsqu'elles sont soumises à l'impact de charges dynamiques. Cette modélisation permet de déterminer les efforts engendrés par l'impact, le but étant de pouvoir étudier certains endommagements comme le délaminage. La nouvelle modélisation utilisée appelée modélisation multiparticulaire des matériaux multicouches (M4), permet de modéliser l'objet tridimensionnel par une géométrie bidimensionnelle tout en augmentant le nombre de champs cinématiques. Elle fait apparaître de manière naturelle des efforts au niveau des interfaces entre les couches. Ces efforts peuvent être reliés au phénomène de délaminage. L'analyse des résultats de la littérature permet de montrer que cette modélisation est pertinente dans le cas des chocs relativement rapides. Un outil numérique a été réalisé sur la base de la méthode des éléments finis en dynamique. Le code éléments finis obtenu est ensuite validé, puis appliqué à l'analyse des efforts dans une plaque multicouche impactée. Nous montrons sur un exemple, la pertinence de cette modélisation et sa capacité à fournir de manière simple et rapide, des informations utiles pour l'étude des différents aspects du délaminage (initiation, propagation, étendue).

La thèse a pour thème la modélisation numérique et la prédiction de l'endommagement (l'usure ou la fissuration) dû au contact avec frottement avec petits débattements (fretting). Elle comporte deux parties. Dans la première partie, une nouvelle méthode de calcul éléments finis est développe pour les chargements cycliques. Contrairement aux algorithmes classiques incrémentaux, cette méthode est basée sur un algorithme original qui permet d'obtenir directement l'état stabilisé mécanique éventuel. La méthode est itérative et consiste à chercher sur tout le cycle de chargement, en un premier temps des solutions statiquement et cinématiquement admissibles, et puis après des solutions plastiquement admissibles. A la fin de chaque itération, la condition de périodicité est imposée. Et ainsi, lorsque l'algorithme converge, on obtient des champs de contraintes et de déformations périodiques et vérifiant l'équilibre et la loi de comportement sur tout le cycle. La comparaison des temps de calculs par rapport à la méthode incrémentale a démontré la performance de cette méthode. Dans la deuxième partie, cette méthode est mise en oeuvre dans le cas du chargement cyclique mobile de fretting. Ainsi des cartes d'adaptation et d'accommodation sont établies en fonction de divers paramètres. Les résultats de prédiction de l'initiation des fissures de fatigue faite par l'application du critère multiaxial de DANG VAN sont en très bon accord avec les résultats expérimentaux. L'usure est obtenue pour une réponse mécanique de type accommodation et elle est donc le résultat d'un mécanisme de fatigue oligocyclique.

En mécanique des matériaux, la volonté actuelle est de chercher à mieux comprendre certains phénomènes macroscopiques en étudiant la microstructure et les phénomènes physiques à l'échelle microscopique. Cette approche est rendue possible par les nombreux développements dans les techniques d'homogénéisation en mécanique multi-échelles. Dans la présente thèse, on travaille sur les élastomères chargés pour lesquels de nombreuses propriétés mécaniques sont étroitement liées à l'arrangement des particules et agrégats de noir de carbone dans la matrice élastomère, à l'échelle microscopique. La démarche adoptée s'articule autour de deux objectifs principaux. Le premier consiste à modéliser la morphologie de la microstructure du matériau. Pour cela, on met en place un modèle aléatoire de morphologie mathématique décrivant chaque échelle de la microstructure. On propose une méthode d'identification de ce modèle à partir de l'exploitation d'images de microscopie à transmission. Cette méthode trouve son originalité dans le fait d'optimiser le modèle en simulant des images de microscopie à transmission, sur lesquelles on mesure les moments statistiques d'ordre deux et trois que l'on compare aux moments expérimentaux. Cette méthode permet finalement de simuler des microstructures dont la morphologie est très proche de celle du matériau réel. Le second objectif consiste, à partir des microstructures ainsi simulées, à déterminer les propriétés effectives du matériau par le calcul par éléments finis à travers la notion de Volume Elémentaire Représentatif (VER). L'idée est de déterminer la taille du VER par une méthode numérique et statistique en cherchant à estimer les propriétés effectives par une approche de type Monte-Carlo, pour des simulations de microstructures de tailles croissantes. La détermination du VER porte sur les modules élastiques et la conductitivité électrique. De nombreux outils tels que le maillage par éléments finis ou encore le calcul parallèle appliqué aux matériaux présentant un fort contraste sur les propriétés entre les phases, ont été explorés afin de mener à bien ce dernier objectif.

Lors de la fabrication d'une bouteille par injection-soufflage, le conditionnement thermique de la préforme joue un rôle essentiel. Nous proposons une procédure d'optimisation numérique permettant de calculer automatiquement les paramètres de réglage du four infrarouge. Dans un premier temps, l'algorithme d'optimisation de Nelder-Mead est couplé avec des simulations éléments finis de l'étape de soufflage, réalisées avec ABAQUS®. L'objectif est de calculer la distribution de température optimale de la préforme, permettant d'uniformiser l'épaisseur de la bouteille. Dans un second temps, un algorithme de programmation quadratique séquentielle est couplé avec un modèle numérique de chauffage infrarouge 3D développé au laboratoire. Cette méthode permet de calculer les paramètres optimaux pour le réglage du four. Des mesures expérimentales réalisées sur une machine de soufflage semi-industrielle ont permis de valider qualitativement notre approche pour une bouteille de géométrie simple. Les propriétés radiatives du PET sont mesurées par spectrométrie infrarouge. Ces mesures sont exploitées pour calculer l'absorption spectrale du rayonnement. Le modèle de chauffage est validé à l'aide de mesures thermographiques. Un capteur a été développé pour mesurer la résistance thermique de contact entre le polymère et le moule. Le débit d'air injecté dans la préforme a été mesuré, puis appliqué en tant de donnée d'entrée. La pression de soufflage est alors automatiquement calculée à chaque itération par un modèle thermodynamique. Les cinématiques de mise en forme, ainsi que les distributions d'épaisseurs calculées par le modèle sont conformes à celles mesurées sur un pilote de laboratoire.

Une bonne compréhension du comportement des radiochronomètres est essentielle pour mener des études de thermochronologie visant à contraindre des processus géodynamiques. Cependant, le comportement des radiochronomètres est complexe car dicté par un grand nombre de facteurs intrinsèques et extrinsèques au système étudié. Au cours de cette thèse, nous avons cherché à déterminer quels étaient les facteurs critiques lors des processus de rééquilibrage isotopique survenant durant un événement métamorphique de haut grade. Ainsi, nous avons étudié le comportement des systèmes isotopiques Rb/Sr, Ar/Ar et Sm/Nd lors de l'éclogitisation partielle (contrôlée par la déformation et la circulation de fluides) au Calédonien de roches granulitiques grenvilliennes dans l'Arc de Bergen en Norvège. Nous montrons, grâce à une étude géochronologique fine associée à différentes observations (pétrologiques, géochimiques in situ et sur fractions minérales, diffraction X et isotopes de l'oxygène) permettant de contraindre avec précision le flux des éléments d'un site textural vers un autre au cours des réactions d'éclogitisation, que le facteur principal garantissant un rééquilibrage isotopique complet pour les systèmes Rb/Sr et Sm/Nd lors d'un événement métamorphique de haut grade est la redistribution homogène des éléments à l'échelle de l'échantillon. Cependant, la circulation de fluides et la déformation qui peuvent être associées aux processus de recristallisation ne sont pas nécessairement suffisant pour garantir une telle réhomogénéisation. La mise en solution des éléments d'un système ou la fusion partielle apparaissent comme deux mécanismes de réhomogénéisation bien plus efficaces. Pour le système Ar/Ar, le facteur critique est l'évacuation totale de l'argon radiogénique précédemment accumulé dans le système. Cette évacuation est favorisée, mais non garantie, par la circulation de fluides en raison de la plus forte solubilité de l'argon dans les fluides. Cependant un grand nombre de facteurs vont contrôler ce processus : quantité d'argon radiogénique initialement présente (fonction de la composition chimique de la roche et de son âge), composition chimique du fluide, et circulation en système ouvert ou fermé.

Les éléments transposables sont des séquences génomiques capables de se répliquer et de se déplacer dans les génomes. Leur capacité à s'insérer près des gènes et à produire des réarrangements chromosomiques par recombinaison entre copies, font des éléments transposables des agents mutagènes. Les éléments transposables sont de plus capables de modifier l'expression des gènes voisins grâce aux régions promotrices qu'ils possèdent. Les éléments transposables ont été trouvés dans la plupart des génomes dans lesquels ils ont été recherchés. Ils forment ainsi 45 % du génome de l'homme et peuvent représenter jusqu'à 90 % du génome de certaines plantes. Dans la première partie de ma thèse, je me suis penché sur les facteurs qui déterminent la distribution de ces éléments. Je me suis intéressé à un facteur particulier, qui est la fonction des gènes dans le voisinage des insertions d'éléments transposables. Dans la deuxième partie, j'ai essayé de déterminer l'impact de l'altération des modifications épigénétiques (modifications d'histones plus précisément) associées aux différents composants géniques, dont les éléments transposables, sur la variation de l'expression des gènes en condition tumorale.