Les lois d'endommagement incrémentales présentent de nombreux avantages en terme de modélisation thermomécanique. S'affranchissant de la notion de cycle en fatigue, elles s'appliquent naturellement aux chargements complexes, anisothermes. Des extensions à l'anisotropie du modèle d'endommagement isotrope de Lemaitre sont proposées et identifiées sur le Haynesl88, superalliage à base cobalt utilisé pour la réalisation de chambres de combustion de turbomachines, permettant de déterminer les conditions d'amorçage de fissure par la méthode des éléments finis. Le cas des hautes températures est traité via un couplage viscoplasticité/endommagement dans le cadre de la thermodynamique des matériaux solides. Le seuil d'endommagement en énergie stockée et l'extension du modèle aux conditions unilatérales de refermeture des micro défauts sont présentés. Différents schémas numériques de résolution sont proposés, dans le cas de calculs couplés, ou de post traitements de calculs 3D viscoplastiques sans endommagement. Des applications anisothermes sont simulées: fluage à température variable, chargements biaxiaux, non proportionnels, aléatoires... Les apports de l'anisotropie de l'endommagement et des conditions unilatérales de refermeture des microdéfauts sont étudiés. Enfin, la robustesse du modèle est testée sur un essai original de fatigue thermique structural sur éprouvette multiperforée, reproduisant les sollicitations subies par une pièce réelle. La corrélation entre résultats numériques et expérimentaux en terme de comportement et de durée de vie permet la validation de la méthodologie d'identification et d'utilisation du modèle incrémental, en vue de son industrialisation.

.

Durant un tremblement de terre, la rupture grandit et se propage sur la faille. Lorsqu'elle s'arrête, le séisme atteint sa taille finale. Comprendre ce qui, dans la nature, détermine la capacité à se propager ou à s'arrêter de la rupture est fondamental en sismologie, puisqu'il existe des séismes de toutes les tailles, mais que seuls les plus grands sont dévastateurs. Dans cette thèse, nous étudions l'impact de différentes façons d'arrêter ou de perturber la propagation de la rupture, à travers des études numériques dynamiques. <br /><br />Tout d'abord, nous avons inclus une limite à l'élasticité du milieu entourant la faille, afin de simuler la propagation de la rupture dans un milieu fracturé, qui dissipe une partie de l'énergie libérée. Nous avons, pour la première fois, inclus et étudié l'impact de cette dissipation dans un modèle de rupture 3D. La rupture, dans ces conditions, est beaucoup plus sensible aux barrières, et s'arrête plus facilement. La cinématique de la rupture est remarquablement modifiée (vitesse de propagation plus lente, vitesse de dislocation maximale limitée). Les mouvements engendrés en surface sont atténués. La plasticité est, de fait, un phénomène crucial à prendre en compte dans la modélisation de la rupture sismique.<br /><br />Dans une deuxième partie, nous avons étudié l'impact d'une hétérogénéité spatiale de résistance à la rupture sur la faille. L'introduction de l'hétérogénéité permet d'obtenir des profils de glissement dont la forme se rapproche des formes observées dans les cas naturels. La propagation et l'arrêt de la rupture perdent leur caractère prédictible lorsque ce sont de petites barrières qui arrêtent la rupture progressivement. On peut obtenir une grande variété de tailles d'événements pour une même statistique de taille des barrières. L'obtention d'une loi puissance de type Gutemberg-Richter sur toute la gamme des tailles d'événements est conditionnée par l'augmentation progressive de l'énergie de fracturation avec la taille de la rupture, d'une façon similaire à ce qui est obtenu en considérant un comportement plastique du milieu. <br /><br />Enfin, l'étude des relations glissement final - taille de l'aspérité rompue, dans des modèles lisses de type aspérité/barrière, a montré que la dynamique contrôlait une partie de la loi d'échelle du glissement maximum, et que la segmentation des failles modifie sensiblement la loi d'échelle.

Nous étudions les propriétés de modèles de séries chronologiques vectorielles à composantes binaires en vue de la modélisation du phénomène climatologique de la succession des jours selon leur caractère sec ou humide dans un réseau de stations de mesures. Nous menons une étude expérimentale par simulations des modèles considérés et envisageons leur utilisation dans des essais de modélisation concernant certains sous-réseaux francais de stations météorologiques

Cette thèse expérimentale porte sur l'évaporation d'une gouttelette mouillante : il s'agit d'un problème de ligne de contact mobile dans lequel deux processus de relaxation aux effets antagonistes (étalement et évaporation) s'opèrent de façon concomitante. <br /> Le travail qui est rapporté fait suite à la thèse de C. Poulard soutenue en 2005 au sein de notre groupe. Il s'est dégagé de ce premier travail de grandes tendances : notamment, la phase de rétraction est caractérisée par l'existence de lois de puissance pour le rayon et l'angle de contact. Ce caractère particulier a donné lieu à un modèle élaboré en collaboration avec M. Ben Amar. Le domaine de validité de ce modèle peut toutefois être amélioré. Le manuscrit prend ainsi la forme d'une discussion illustrée de faits expérimentaux, où les hypothèses sont reprises, les processus invoqués, testés et au cours de laquelle l'accent est mis sur les points importants en vue d'une meilleure description de la dynamique de ces objets volatils.<br /> Soulignons qu'en raison de la forme particulière du taux d'évaporation, la ligne de contact recule alors que l'angle de contact n'est pas nul. Or, cet angle est particulièrement significatif dans les dynamiques. Un accent particulier est porté sur la forme de l'interface liquide au niveau de la ligne de contact, gouvernée par un réel problème d'hydrodynamique. Notamment, des tests précis sur le couplage entre l'évaporation et les écoulements hydrodynamiques à l'extension maximale sont menés et les résultats de cette étude sont appliqués avec succès à la remise à l'échelle des dynamiques.<br /> Nous abordons dans une seconde partie, la question peu typique de l'évaporation d'une goutte de mélange. Dans un tel système, l'évaporation provoque un défaut du composé le plus volatil au niveau de la ligne de contact. Nous verrons dans cette partie que les hétérogénéités ainsi créées sont à l'origine de gradients de tension de surface qui contrôlent les dynamiques.

Nous étudions des méthodes d'approximation particulaire de solutions d'équations aux dérivées partielles décrivant l'état macroscopique de certains systèmes physiques. Elles consistent en l'introduction d'un grand nombre N de particules fictives évoluant selon des équations différentielles couplées, ordinaires ou stochastiques, dans un sens plus simple à résoudre que l'équation macroscopique; l'état de ce système de particules est décrit par une mesure de probabilité, dite mesure empirique. La validité de la méthode est donnée par la convergence, quand N tend vers l'infini, de cette mesure empirique vers la solution macroscopique originale, appelée limite de champ moyen. Nous cherchons principalement à en donner des estimations explicites, quantifiant ainsi la précision de l'approximation.<br /><br />Dans ce cadre nous étudions l'approximation des équations de transport de Vlasov et d'Euler par des systèmes de particules déterministes en interaction. Le problème de la convergence de la méthode se ramène à un problème de stabilité de solutions que nous traitons par des propriétés de type contraction pour des distances (de Wasserstein) liées à la théorie du transport optimal de mesures. Nous établissons aussi une propriété analogue de contraction pour des lois de conservation scalaires. <br /><br />Nous étudions également l'approximation d'équations de diffusion de McKean-Vlasov par des systèmes de particules stochastiques. Nous en donnons l'erreur de manière quantitative à l'aide de techniques de couplage, d'estimations de propagation du chaos et d'inégalités de concentration ou de déviation.<br /><br />De façon plus systématique nous nous intéressons à de telles inégalités de concentration pour des mesures de probabilité et à leurs relations avec des inégalités de transport (liant distances de Wasserstein et entropie) et de Sobolev logarithmiques. En particulier nous établissons de telles inégalités pour certaines classes de lois de variables dépendantes.

Le but de cette thèse est de trouver une méthodologie numérique pour simuler le comportement et la fragmentation des matériaux quasi-fragiles, tels que les roches et le béton, sous sollicitations dynamiques. Une telle méthodologie aidera à une meilleure compréhension des mécanismes de fragmentation des roches à l'explosif et, par conséquent, contribuera dans l'optimisation des performances du tir à l'explosif. Dans cette étude, on considère que le processus de fragmentation est une extension naturelle de celui de la rupture. Afin de décrire ce dernier processus, un modèle de comportement phénoménologique adapté aux matériaux considérés est développé et implémenté dans un code de calcul par éléments finis. Quant à la fragmentation dynamique, elle est traitée par une analyse en post-traitement basée sur l'historique de l'état thermodynamique du matériau. Pour reproduire le comportement macroscopique des matériaux quasi-fragiles incluant l'anisotropie induite par le chargement, la méthode des variables internes basée sur la thermodynamique des milieux continus est employée. Une variable interne de type scalaire est introduite pour modéliser, au niveau macroscopique, l'adoucissement du matériau suite à des chargements de compression. Sous des chargements de traction, un tenseur symétrique du second ordre est utilisé pour décrire l'endommagement anisotrope induit. Sous des chargements complexes, ces deux modèles sont couplés et l'effet de fermeture-réouverture des fissures est également traité. L'intégration du modèle développé dans le code éléments finis VIPLEF3D a conduit à l'élaboration d'une méthode de relaxation caractérisée par une actualisation explicite des variables d'état. Concernant la fragmentation dynamique, à partir d'essais de fragmentation au laboratoire, une formulation générale permettant de prédire la distribution des fragments est fournie. Dans cette formulation une taille moyenne de fragments est liée à une grandeur mécanique, supposée être l'origine du processus de fragmentation, par une fonction intrinsèque qui peut être identifiée en utilisant des essais adéquats de fragmentation. Cette grandeur mécanique est donnée par la résolution du problème aux limites où le modèle rhéologique est utilisé. Enfin, l'approche complète est alors appliquée à la modélisation des essais de rupture et de fragmentation en chambre, d'échantillons cylindriques de roches.

Ce mémoire présente une étude des propriétés de polymères à petite (nano) échelle, par les modes de modulation d'amplitude (« Tapping ») et de modulation de fréquence (« Non-contact » résonant) de la Microscopie à Force Atomique. Il débute par la détermination des origines physiques du contraste d'images obtenues sur une surface de copolymère tri- séquencé : la reconstruction d'une section d'image grâce aux courbes d'approche- retrait permet d'en attribuer le contraste aux seules variations de la réponse viscoélastique du copolymère (i.e. de montrer l'absence de contribution topographique au contraste). L'étape suivante est une évaluation des constantes physiques des deux nanophases du même copolymère tri- séquencé par l'analyse de courbes d'approche- retrait en mode modulation d'amplitude. Pour cela nous modélisons analytiquement et numériquement l'interaction non-linéaire et dissipative entre nano-pointe oscillante et polymère. Puis un résultat expérimental charnière montre qu'en mode modulation d'amplitude avec un facteur de qualité Q suffisamment élevé un matériau de raideur locale très faible peut apparaître dur et sa topographie être mesurée sans déformation. Cela signifie que, sur un échantillon suffisamment mou, la variation contrôlée du facteur de qualité Q permet d'accéder sélectivement aux propriétés mécaniques ou topographiques. La dernière partie concerne la dépendance d'échelle de la viscoélasticité d'un polymère fondu enchevêtré. On utilise le mode modulation de fréquence avec un facteur Q très élevé : lorsque l'indentation de la pointe dans le polymère et leur temps d'interaction atteignent des valeurs seuils, la puissance nécessaire pour maintenir une amplitude d'oscillation constante présente une augmentation abrupte. On montre que cette augmentation de la dissipation est due au passage d'une relaxation moléculaire rapide (<10-8 s, mouvement de type Rouse à faible indentation) à une relaxation beaucoup plus lente (>10-5s, désenchevêtrement lorsque la taille de l'objet indenté dépasse celle de la maille du réseau de chaînes). Cette transition correspond à une variation de profondeur d'indentation de quelques angströms.

La Simulations des Grandes Echelles d'écoulements turbulents anisothermes a été utilisée dans le cadre de problèmes d'interaction thermique fluide/solide. Le but de ce travail réside non seulement dans l'identification des divers éléments susceptibles de mésestimer les fluctuations de température à la paroi mais aussi à l'introduction de modélisations de parois adaptées. Le choix d'un schéma de convection "optimisant" le transport de scalaire a conduit à l'adoption d'un schéma décentré avec limiteur de pente : le schéma QUICK. Deux nouvelles approches de reconstruction de fluctuations de température en proche paroi sont proposées et testées. La première repose sur une résolution complète des équations de Navier-Stokes sur une grille fine en proche paroi permettant de reconstruire le champ fluctuant. Une seconde méthode repose sur la résolution 1D et simultanée d'une équation de température moyenne et d'une équation de transport pour la variance de température en proche paroi.

Nous étudions l'approximation cinétique discrète de lois de conservation scalaires quasi-linéaires dans le quart d'espace positif. Cette approximation est obtenue par l'introduction de systèmes de type BGK relaxant la loi scalaire. Nous démontrons la convergence des systèmes semi-linéaires vers la loi scalaire. Nous discrétisons ces modèles pour obtenir une gamme de schémas numériques adaptés au problème avec bord. Dans une troisième partie, nous appliquons ces schémas à un certain nombre de cas test numériques.