Reproduction de : Thèse de doctorat : Génie mécanique, productique et transports : Poitiers : 2008. / Titre provenant de l'écran-titre. Bibliogr. 69 réf.

Thèse doctorat : Mécanique : Villeurbanne, INSA : 2005. / Titre provenant de l'écran-titre. Bibliogr. p. 203-213.

Dans ce travail de thèse, une stratégie de calcul multi-échelle en dynamique transitoire basée sur la méthode LATIN est proposée. Ce travail fait suite, entre autre, à la thèse de H. Lemoussu qui a appliqué la méthode LATIN dans sa version mono-échelle au cas de la dynamique, ainsi qu'aux avancées plus récentes concernant l'introduction d'une vision à deux échelles au sein de la méthode LATIN pour des calculs statiques et quasi-statiques. Notre but a été d'étendre cette vision à deux échelles au cas de la dynamique transitoire. Une écriture de la stratégie de résolution multi-échelle est proposée pour le cas de la dynamique. Ce travail a permis de mettre en évidence plusieurs particularités de la méthode au sujet des conditions d'admissibilité des quantités macroscopiques en dynamique. L'introduction de l'approche multi-échelle en dynamique a permis d'obtenir l'extensibilité numérique de la méthode de décomposition de domaine pour le cas des interfaces. L'autre volet de cette thèse concerne la prise en compte de variabilité paramétrique en dynamique transitoire. Le but était ici de s'appuyer sur le savoir faire du LMT-Cachan, en matière de technique de calcul multi-résolution pour des problèmes avec contact pour développer une stratégie de calcul multi-résolution multi-échelle en dynamique, apte à prendre en compte les incertitudes, tout en diminuant de façon drastique le coût de calcul par rapport aux approches conventionnelles.

Les instabilités générées par frottement sont responsables des divers bruits tels que le crissement, le sifflement ou le broutement ... Pour modéliser et comprendre ce phénomène d'instabilités, les analyses temporelle et fréquentielle sont utilisées sur un système modèle constitué de deux poutres en contact. Dans l'analyse fréquentielle, linéaire, l'instabilité se manifeste par la coalescence de deux modes propres du système. Dans l'analyse temporelle, qui tient compte de l'aspect non linéaire d'un contact frottant, l'instabilité est caractérisée par des zones d'adhérence ou de décollement qui apparaissent au niveau de la surface de contact. Les résultats issus des deux analyses sont cohérents et complémentaires, malgré quelques différences de prédiction. Une validation expérimentale a été effectuée et montre une bonne corrélation entre les résultats numériques temporels et expérimentaux. On met en évidence la pertinence de l'analyse temporelle dans l'étude des instabilités de contact, phénomène vibratoire complexe. On montre également que, même si elle n'apparaît pas suffisante pour caractériser le phénomène d'instabilité de contact, l'analyse fréquentielle donne de bons résultats. L'aspect tribologique est également abordé et met en évidence une interaction permanente entre les phénomènes aux échelles micro- et macroscopique.

Thèse doctorat : Mécanique : Lille 1 : 2003. / N° d'ordre (Lille 1) : 3319. Résumé en français et en anglais. Bibliogr. p. 127-131.

Ce travail s'inscrit comme contribution au calcul numérique des problèmes de contact-impact entre solides déformables en grandes transformations. Il comporte trois parties. Dans la première, nous élaborons, grâce notamment à l'usage de "champs de signes" (type Level-Sets), une formulation lagrangienne originale pour les problèmes de contact, qualifiée de stabilisée, permettant de généraliser les formulations lagrangiennes classique et augmentée, tout en clarifiant l'intérêt des paramètres d'augmentation et unifiant les implémentations numériques. La deuxième partie est focalisée sur les problèmes d'impact. Nous y proposons une formulation continue hybride (en vitesse) faible-forte dérivant du modèle de Signorini-Moreau, écrit en équations, moyennant l'introduction de champs de signes inconnus. De la formulation continue sont dérivés des éléments de contact-impact après discrétisations en temps par une variante de Θ-schéma et en espace par la méthode des éléments finis compatibles et une méthode de collocation (points finis). La troisième partie est centrée sur les aspects multi-échelles des problèmes de contact-impact. Nous y proposons, tout particulièrement, un modèle d'interface multi-niveau rendant compte des comportements locaux et globaux des interfaces de contact. Ce modèle permet d'intégrer une réalité physique du contact, tout en atténuant significativement les écueils numériques de conditionnement et d'oscillations parasites. Par ailleurs, l'approche Arlequin a été appliquée aux problèmes multi-échelles (en espace) de contact-impact montrant la possibilité de "zoomer" avec une grande flexibilité les zones critiques de contact-impact et de mixer des modèles et des schémas différents au sein d'une même structure impactée. La pertinence des idées proposées est illustrée par des résultats numériques obtenus suite à des développements réalisés dans Code_Aster.

Ce travail de thèse concerne la mécanique numérique du contact entre solides déformables. Il s'agit de contribuer à l'amélioration de la performance, de la fiabilité et de la robustesse des algorithmes et des modèles numériques utilisés dans les codes éléments finis en particulier Code_Aster qui est un code libre développé par Électricité De France (EDF) pour ses besoins en ingénierie. L'objectif final est de traiter les problèmes industriels de grande dimension avec un temps de calcul optimisé. Pour parvenir à ces objectifs, les algorithmes et formulations doivent prendre en compte les difficultés liées à la mécanique non régulière à cause des lois de Signorini-Coulomb ainsi que la gestion des non linéarités dûes aux grandes déformations et aux comportements des matériaux étudiés.Le premier axe de ce travail est dédié à une meilleure compréhension de la formulation dite de « Lagrangien stabilisé » initialement implémentée dans le code. Il a été démontré l'équivalence entre cette formulation et la formulation bien connue de « Lagrangien augmenté ». Les caractéristiques mathématiques liées aux opérateurs discrets ont été précisées et une écriture énergétique globale a été trouvée. / This work deals with computational contact mechanics between deformable solids. The aim of this work is to improve the performance, the reliability and the robustness of the algorithms and numerical models set in Code_Aster which is finite element code developped by Électricité De France (EDF) for its engineering needs. The proposed algorithms are used to solve high dimensional industrial problems in order to optimize the computational running times. Several solutions techniques are available in the field of computational contact mechanics but they must take into account the difficulties coming from non-smooth aspects due to Signorini-Coulomb laws coupled to large deformations of bodies and material non linearities. Firstly the augmented Lagrangian formulation so-called « stabilized Lagrangian » is introduced. Successively, the mathematical properties of the discrete operators are highlighted and furthermore a novel energetic function is presented. Secondly the kinematical condition with regard to the normal unknowns are reinforced through unconstrained optimization techniques which result to a novel formulation which is so-called « non standard augmented Lagrangian formulation ». Three types of strategies are implemented in the code. The generalized Newton method is developped : it is a method in which all the non linearities are solved in one loop of iterations. The partial Newton method is an hybrid technique between the generalized Newton one and a fixed point method.

Afin d'étudier la technologie du stockage géologique du CO2, l'intégrité dans le temps de la formation hôte doit être assurée, notamment par l'évaluation de l'altération des propriétés mécaniques des géomatériaux en présence de CO2. Le scénario à l'origine des ces travaux de thèse illustre la possibilité d'une remontée de CO2 le long d'une faille de la roche de couverture. Ce problème géologique est complexe en raison du nombre de paramètres à considérer : température, confinement, effet d'échelle, hétérogénéités du géomatériau, composition du géofluide, réactions chimiques... Ces travaux se cantonnent à l'étude du comportement d'une discontinuité dans un géomatériau sous sollicitations chemo-mécaniques à l'échelle du laboratoire. Ils ont été menés en cotutelle entre le laboratoire SIAME de l'Université de Pau et des Pays de l'Adour (France) et le laboratoire de mécanique des roches et de géologie appliquée de l'Université de Sherbrooke (Québec, Canada). Le programme expérimental a été défini pour caractériser d'une part la refermeture d'une fissure sous effort uni-axial cyclique, et d'autre part le comportement en cisaillement d'un joint rocheux dégradé chimiquement. Les données obtenues ont été utilisés pour modéliser le comportement d'une discontinuité. L'étude du comportement d'une fissure sous effort normal a été conduite par un essai cyclique de traction-compression. Il précise l'évolution des déformations inélastiques générées lors de cycles de fermeture-ouverture d'une discontinuité dans un géomatériau de synthèse (béton). Les données obtenues confirment la présence d'une hystérésis lors des cycles de chargement. L'analyse du champ de déplacements par corrélation d'images lors de la refermeture de fissure permet d'en déduire que les déformations inélastiques seraient gouvernées en partie par les frottements générés lors du ré-emboîtement des lèvres de la discontinuité. La non-correspondance des profils rugueux, entraînant la dissipation d'énergie par frottements illustrée par l'hystérésis, serait directement liée à la présence de contraintes internes. Le comportement d'un joint rocheux altéré sous sollicitation tangentielle utilise l'essai de cisaillement direct. Les surfaces rugueuses de roche sont dégradées chimiquement par immersion dans des solutions acides. La numérisation de ces surfaces par un profilomètre laser indique de légères modifications de la rugosité par dissolution. Les essais révèlent une profonde modification du comportement des joints attaqués : le pic de contrainte disparaît, la rigidité du joint diminue avec l'angle de dilatance et la phase de contractance est accrue. Ces évolutions s'expliquent par l'endommagement chimique qui accentue la non-correspondance des profils rugueux des joints et diminue les propriétés mécaniques du matériau de part et d'autre de la discontinuité. D'un point de vue numérique, l'apport de ces travaux de thèse réside dans la modélisation du comportement d'une discontinuité par le couplage d'un modèle élasto-plastique endommageable continu avec une résolution discrète du problème de contact/frottement. Les modélisations sont réalisées avec le code de calcul aux éléments finis Cast3M. Les géométries des lèvres de discontinuités sont directement modélisées à partir des profils de rugosité issus des numérisations expérimentales. Les résultats numériques représentent correctement les phénomènes de frottements constatés expérimentalement. Enfin, une modélisation de l'essai de cisaillement des joints altérés est réalisée en couplant le modèle mécanique avec un modèle d'endommagement chimique.

Résumé : Ces travaux de thèse s'intéressent à l'étude du comportement d'une discontinuité dans un géomatériau sous sollicitations chemo-mécaniques à l'échelle du laboratoire. Des essais de traction-compression cyclique étudient la refermeture d'une fissure. Ils indiquent que les déformations inélastiques seraient gouvernées en partie par les frottements générés lors du ré-emboîtement des lèvres de la discontinuité, non correspondantes du fait des contraintes internes. Des joints rocheux altérés chimiquement sous sollicitation tangentielle sont étudiés au travers d'essais de cisaillement direct : le comportement des joints dégradés est profondément modifié du fait de la diminution des propriétés mécaniques du matériau de part et d'autre de la discontinuité et de l'accentuation de la non-correspondance des profils rugueux. Le comportement d'une discontinuité est modélisée par le couplage d'un modèle élasto-plastique endommageable continu avec une résolution discrète du problème de contact/frottement (code calcul aux éléments finis Cast3M). Les résultats numériques confirment les phénomènes constatés expérimentalement.//Abstract: To probate the technology of CO[indice inférieur 2] geological storage, the integrity of the site must be assure over time. This industruial problematic involves the study of the mechanical properties alteration of geomaterials in the presence of CO[indice inférieur 2]. The scenario at the origin of this thesis illustrates the possibility of a CO[indice inférieur 2] leakage on a fault located in the caprock. This geological problem is complicated by the many parameters to consider: in situ temperature and pressure, scale effect, heterogeneities of the geomaterial, geofluide composition, chemical reactions ... These works focus on the behavior of a discontinuity in a geomaterial solicited chemomechanically at the laboratory scale. They were realised in cotutelle between SIAME laboratory at the University of Pau and Pays de l'Adour (France) and the laboratory of rock mechanics and engineering geology from the University of Sherbrooke (Quebec, Canada). The first part of the experimental program was defined to characterize the reclosing of a crack under cyclic uniaxial stress. The second experimental campaign has studied the shear behavior of a rock joint chemically degraded. The data obtained were used to model the behavior of a discontinuity by the finite element method. The mechanical behavior of a crack under normal stress is assessed with cyclic tension-compression tests. Stress curve showed hysteresis during opening and closing cycles of a discontinuity in concrete, it indicated inelastic deformations The analysis of displacement field by image correlation indicated that theses deformations were partially governed by the friction generated during the closing of the discontinuity lips. Frictional phenomena are due to asperities mismatching induced by the internal stresses in the concrete. The shear behavior of a rock joint chemically damaged was studied through direct shear tests. Rough surfaces were immersed in acid solution during 6 hours at constant pH. Digitalization of these surfaces befor and after immesion, with a lase profilometer, indicates little modifications of the geometry induced by dissolution of material. Results of tests pointed out significant modifications for altered joints illustred by a of the peak shear strength and an increased of contractancy. They are induced by: 1) the mismatch enhancement of the rough profiles of the discontinuity and, 2) the degradation of the mechanical properties of the material on both sides of the discontinuity due to the chemical attack. Numerical contribution of the thesis lies in modeling the behavior of a discontinuity by the coupling of an continuous elastic-plastic damaged model with a discrete resolution of the contact/friction problem. The model is developed with the finite element code Cast3M. Geometries lips discontinuities are modeled directly from the roughness profiles from experimental scans. The numerical results correctly represent the friction phenomena observed experimentally. Finally, a model of the shear test altered joints is performed by coupling the mechanical model with chemical damage model.

Ces travaux de thèse s'intéressent à l'étude du comportement d'une discontinuité dans un géomatériau sous sollicitations chemo-mécaniques à l'échelle du laboratoire. Des essais de traction-compression cyclique étudient la refermeture d'une fissure. Ils indiquent que les déformations inélastiques seraient gouvernées en partie par les frottements générés lors du ré-emboîtement des lèvres de la discontinuité, non correspondantes du fait des contraintes internes. Des joints rocheux altérés chimiquement sous sollicitation tangentielle sont étudiés au travers d'essais de cisaillement direct : le comportement des joints dégradés est profondément modifié du fait de la diminution des propriétés mécaniques du matériau de part et d'autre de la discontinuité et de l'accentuation de la non-correspondance des profils rugueux. Le comportement d'une discontinuité est modélisée par le couplage d'un modèle élasto-plastique endommageable continu avec une résolution discrète du problème de contact/frottement (code calcul aux éléments finis Cast3M). Les résultats numériques confirment les phénomènes constatés expérimentalement. / This PhD dissertation presents a study aiming at a better understanding of cracks behavior in a geomaterial. The study focuses on the behavior of discontinuities under chemo-mechanical solicitations at the laboratory scale. The mechanical behavior under normal stress is assessed with cyclic tension-compression tests. Experimental data indicate that the inelastic deformations could be partially governed by the friction generated during the closing of the discontinuity lips, and the asperities mismatch is related to the internal stresses. The shear behavior of a rock joint chemically degraded was studied through direct shear tests. Results pointed out significant modifications for altered joints induced by: 1) the mismatch enhancement of the rough profiles of the discontinuity and, 2) the degradation of the mechanical properties of the material on both sides of the discontinuity due to the chemical attack. These experimental results have been used as input data to model the behavior of a discontinuity by the coupling of a continuous elastic-plastic damaged model with a discrete solving of the contact/friction problem. The simulations performed under Cast3M correctly represent the phenomena observed during the experimental testing program.