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.

La mousse de polyuréthane est un matériau cellulaire caractérisé par un spectre de propriétés mécaniques intéressant : une faible densité, une capacité à absorber l’énergie de déformation et une faible raideur.Elle présente également des propriétés telles qu’une excellente isolation thermique et acoustique, une forte absorption des liquides et une diffusion complexe de la lumière. Ce spectre de propriétés fait de la mousse de polyuréthane un des matériaux couramment utilisés dans de nombreuses applications phoniques, thermiques et de confort. Pour contrôler la vibration transmise aux occupants des sièges, plusieurs dispositifs automatiques de régulation et de contrôle sont actuellement en cours de développement tels que les amortisseurs actifs et semi-actifs. La performance de ces derniers dépend bien évidemment de la prédiction des comportements de tous les composants du siège et en particulier la mousse. D’une façon générale, il est indispensable de modéliser le comportement mécanique complexe de la mousse de polyuréthane et d’identifier ses propriétés quasi-statique et dynamiques afin d’optimiser la conception des systèmes incluant la mousse en particulier l’optimisation de l’aspect confort. Dans cette optique, l’objectif principal de cette thèse consiste à implémenter des modèles mécaniques de la mousse de polyuréthane fiables et capables de prévoir sa réponse sous différentes conditions d’essais. Dans la littérature, on retrouve les divers modèles développés tels que les modèles de mémoire entier et fractionnaire. L’inconvénient majeur de ces modèles est lié à la dépendance de leurs paramètres vis-à-vis des conditions d’essais, chose qui affecte le caractère général de leur représentativité des comportements quasi-statique et dynamique de la mousse polyuréthane. Pour pallier à cet inconvénient, nous avons développé des modèles qui, grâce à des choix judicieux de méthodes d’identification, assurent une représentativité plus générale des comportements quasi-statique et dynamique de la mousse polyuréthane. En effet, nous avons démontré qu’on peut exprimer les paramètres dimensionnels des modèles développés par le produit de deux parties indépendantes ; une regroupant les conditions d’essais et une autre définissant les paramètres adimensionnels et invariants qui caractérisent le matériau. Ces résultats ont été obtenus à partir de plusieurs études expérimentales qui ont permis l’appréhension du comportement quasi-statique (à travers des essais de compression unidirectionnelle) et dynamique (à travers des tests en vibration entretenue). La mousse, sous des grandes déformations, présente à la fois un comportement élastique non linéaire et un comportement viscoélastique. En outre, une discrimination entre les modèles développés particulièrement en quasi-statique a été effectuée. Les avantages et les limites de chacun y ont été discutés. / Polyurethane foam is a cellular material characterized by an interesting mechanical spectrum of properties: low density, capacity to absorb the deformation energy and low stiffness. It presents also several other properties, such as excellent thermal and acoustic insulation, high absorption of fluids and a complex scattering of light. This spectrum of properties makes polyurethane foam commonly used in many thermal, acoustic and comfort applications. To control the vibration transmitted to the seat occupants, several automatic devices for regulation and control are currently outstanding developments like active and semi-active dampers. The performance of these devices depends, of course, on the prediction of the behaviour of all the seat components and especially foam. Generally, it is essential to model the complex mechanical behaviour of polyurethane foam and identify its quasi-staticand dynamicproperties in order to optimize the design of systems with foam particularly the optimization of the comfort aspect. In this mind, the main goal of this thesis is to implement mechanical models of polyurethane foam reliable and able to provide its response under different test conditions. Several models has been developed in literature such as memory fractional and integer models. The main disadvantage of these models is the dependence of their parameters against the test conditions. It affects the general character of their representativeness to the quasi-static and dynamic behaviours of polyurethane foam. To solve this problem, we developed models with specific identification methods to ensure broader representation of the quasi-static and dynamic behaviour of polyurethane foam. Indeed, we have demonstrated that we can express the dimensional parameters of the developed models by the product of two independent parts; the first contain only the test conditions and the second define the dimensionless and invariant parameters that characterize the foam material. The developed models have been establish after several experimental studies allowing the apprehension of the quasi-static behaviour (through unidirectional compression tests) and the dynamic behaviour (through harmonic vibration tests). The polyurethane foam, under large deformations, exhibits a non linear elastic behaviour and viscoelastic behaviour. In addition, discrimination between the models developed especially in quasi-static case has been conducted. The advantages and limitations of each model have been discussed.

La mousse flexible de polyuréthane est couramment utilisée dans nombreuses applications telles que acoustiques, thermiques et de bâtiment grâce à sa faible densité et à son pouvoir d’isolation thermique et acoustique. Elle est également utilisée dans les applications de confort pour les sièges tels que véhicule, train, avion etc. grâce sa faible raideur et à son pouvoir à absorber l’énergie de déformation. Pour optimiser le confort des systèmes d’assise, il est nécessaire de modéliser le siège et en particulier la partie flexible, c’est-à-dire la mousse de polyuréthane. Les objectifs principaux de cette thèse consistent à identifier puis à modéliser le comportement quasi-statique de la mousse de polyuréthane sous différentes conditions d’essais sous grandes déformations. Des essais de compression/décompressions unidirectionnels monocycle et multicycle à différentes vitesses de déformations ont été réalisés sur trois types de mousse de polyuréthane, afin de comprendre le comportement du matériau. Ces essais ont permis de déduire que les mousses de polyuréthanes sous grandes déformations présentent à la fois un comportement hyperélastique et un comportement viscoélastique. Ils ont également montrés que les mousses de polyuréthanes présentent un phénomène d’assouplissement appelé ‘effet de Mullins’ lors que les essais de compression/décompressions multicycle, c’est-à-dire que les contraintes dans 1er cycle sont moins faibles que les contraintes dans les cycles suivants pour une même déformation. Sur la base des résultats d’expérimentaux et afin de modéliser le comportement quasi-statique de la mousse de polyuréthanne, nous avons développé trois modèles visco-hyperélastiques qui se composent de deux éléments à savoir la partie modèles énergétiques hyperélastiques, utilisés généralement pour des matériaux à comportement caoutchoutique, et la partie modèle à mémoire entier qui tient compte de l’historique et permettant de décrire le comportement viscoélastique. Les paramètres des modèles ont été identifiés en utilisant la méthode d’identification et la méthode d’optimisation appropriée. Les résultats des modélisations du comportement mécanique de la mousse sur les essais monocycles et multicycles ont été comparés aux résultats expérimentaux, monteront à la fois une très bonne capacité à simuler le 1er cycle de charge/décharge, ainsi que les cycles suivant. Nos modèles ont prouvé leur capacité à modéliser l’effet de Mullins sur les mousses de polyuréthane souple. Ces modèles ont été validés sur les trois types de mousse et pour trois vitesses de sollicitation, permettent de conclure leurs efficacités et de leurs représentativités. / Flexible polyurethane foam is widely used in numerous applications such as acoustic, thermal and building due to its low density and its ability to absorb thermal and acoustic energy. It is also used for the comfort of the seats such as the vehicle, train, plane due to its low stiffness and its ability to absorb deformation energy. In order to optimize the comfort of the car seat, it is necessary to model the behavior of seat system, particularly the flexible component - polyurethane foam. The main objective of this study is to identify and model the quasi-static behavior of polyurethane foam under different test conditions in large deformations. Compression / decompression uniaxial unicycle and multicycle tests were carried out on three types of polyurethane foam at different strain rates, which allows us to understand the behavior of the material. The results of the tests indicate that the polyurethane foams exhibit a hyperelastic behavior and a viscoelastic behavior under large deformations. They also showed that the polyurethane foams have a stress softening phenomenon which is called 'Mullins effect' during the compression / decompression multicycle tests. In other words, the stress in first cycle is lower than the stresses in the subsequent cycles in the same deformation. ‘Mullins effect’ for the polyurethane foam is also an important study in this dissertation. Based on the experimental results and the goal of modeling quasi-static behavior of the polyurethane foam, three visco-hyperelastic models were developed. These models consist in two elements: hyperelastic models, which is normally used for description the behavior of rubber materials, and entire memory model which takes into account the history and describing the viscoelastic behavior. Model parameters were identified using appropriate identification and optimization methods. The results of modeling the mechanical behavior of the foam on the unicycle and multicycle tests were compared with experimental results. The models showed a very good competence to simulate the first cycle and the following cycles during the charge / discharge tests. Our models have proven its ability to model Mullins effect on flexible polyurethane foams. These models have been validated on three types of foam in order to present a comparative study of their effectiveness and their representativeness.

Dans le cadre de la protection du bâti face au risque sismique, les techniques d'analyse simplifiées, basées sur des calculs quasi-statiques en poussée progressive, se sont fortement développées au cours des deux dernières décennies. Le travail de thèse a pour objectif d'optimiser une stratégie d'analyse simplifiée proposée par Chopra et al. (2001) et adoptée par les normes américaines FEMA 273. Il s'agit d'une analyse modale non linéaire découplée, dénommée par les auteurs UMRHA qui se caractérisent principalement par : des calculs de type pushover selon les modes de vibration dominants de la structure, la création de modèles à un degré de liberté non linéaire à partir des courbes de pushover, puis le calcul de la réponse temporelle de la structure en recombinant les réponses temporelles associées à chaque mode de vibration. Dans ce travail, la méthode UMRHA a été améliorée en investiguant les points suivants. Tout d'abord, plusieurs modèles à un degré de liberté non linéaire déduits des courbes de pushover modal sont proposés afin d'enrichir la méthode UMRHA originelle qui emploie un simple modèle élasto-plastique : autres modèles élasto-plastiques avec des courbes enveloppes différentes, le modèle de Takeda prenant en compte un comportement hystérétique propre aux structures sous séismes, et enfin, un modèle simplifié basé sur la dégradation de fréquence en fonction d'un indicateur de dommage. Ce dernier modèle à un degré de liberté privilégie la vision de la chute de fréquence au cours du processus d'endommagement de la structure par rapport à une description réaliste des boucles d'hystérésis. La réponse totale de la structure est obtenue en sommant les contributions non linéaires des modes dominants aux contributions linéaires des modes non dominants. Enfin, la dégradation des déformées modales, due à l'endommagement subi par la structure au cours de la sollicitation sismique, est prise en compte dans la méthode M-UMRHA proposée dans ce travail, en généralisant le concept précédent de dégradation des fréquences modales en fonction d'un indicateur de dommage : la déformée modale devient elle-aussi dépendante d'un indicateur de dommage, le déplacement maximum en tête de l'ouvrage ; l'évolution de la déformée modale en fonction de cet indicateur est directement identifiée à partir des calculs de pushover modal. La pertinence de la nouvelle méthode M-UMRHA est investiguée pour plusieurs types de structures, en adoptant des modélisations éprouvées dans le cadre de la simulation des structures sous séismes : portique en béton armé modélisé par des éléments multifibres pour le béton et les armatures, remplissage en maçonnerie avec des éléments barres diagonales résistant uniquement en compression, bâti existant contreventé (Hôtel de Ville de Grenoble) avec des approches coques multicouches. Les résultats obtenus par la méthode simplifiée proposée sont comparés aux résultats de référence issus de l'analyse temporelle non linéaire dynamique.

Cette thèse vise à étudier le comportement de revêtement articulé du tunnel en développant une nouvelle approche numérique à la Méthode de Réaction Hyperstatique (HRM) et la production des modèles numériques en deux dimensions et trois dimensions à l'aide de la méthode des différences finies (FDM). L'étude a été traitée d'abord sous charges statiques, puis effectuée sous charges dynamiques. Tout d'abord, une étude bibliographique a été effectuée. Une nouvelle approche numérique appliquée à la méthode HRM a ensuite été développée. En même temps, un modèle numérique en deux dimensions est programmé sur les conditions de charge statique dans le but d'évaluer l'influence des joints, en termes de la distribution et des caractéristiques des joints, sur le comportement du revêtement articulé de tunnel. Après cela, des modèles complets en trois dimensions d'un seul tunnel, de deux tunnels horizontaux et de deux tunnels empilés, dans lesquels le système des joints est simulé, ont été développés. Ces modèles en trois dimensions permettent d'étudier le comportement non seulement du revêtement du tunnel, mais encore le déplacement du sol entourant le tunnel lors de l’excavation. Un modèle numérique en trois dimensions simplifié a ensuite été réalisé afin de valider la nouvelle approche numérique appliquée à la méthode HRM.Dans la dernière partie de ce mémoire, la performance du revêtement articulé du tunnel sous chargements dynamiques est prise en compte par l’analyse quasi-statique et dynamique complète en utilisant le modèle numérique en deux dimensions (FDM). Un modèle HRM a également été développé prenant en compte des charges quasi-statiques. Les différences de comportement de tunnel sous chargements statiques et sismiques sont mises en évidence et expliquées. / This PhD thesis has the aim to study the behaviour of segmental tunnel lining by developing a new numerical approach to the Hyperstatic Reaction Method (HRM) and producing two-dimensional (2D) and three-dimensional (3D) numerical models using the finite difference method (FDM). The study first deals with under static loads, and then performs under dynamic loads. Firstly, a literature review has been conducted. A new numerical approach applied to the HRM has then been developed. At the same time, a 2D numerical model is programmed regarding static loading conditions in order to evaluate the influence of the segmental joints, in terms of both joint distribution and joint stiffness characteristics, on the tunnel lining behaviour. After that, full 3D models of a single tunnel, twin horizontal tunnels and twin tunnels stacked over each other, excavated in close proximity in which the joint pattern is simulated, have been developed. These 3D models allow one to investigate the behaviour of not only the tunnel lining but also the displacement of the ground surrounding the tunnel during the tunnel excavation. A simplified 3D numerical model has then been produced in order to validate the new numerical approach applied to the HRM. In the last part of the manuscript, the performance of the segmental tunnel lining exposed to dynamic loading is taken into consideration through quasi-static and full dynamic analyses using 2D numerical models (FDM). A new HRM model has also been developed considering quasi-static loads. The differences of the tunnel behaviour under static and seismic loadings are highlighted.

Ce travail s'inscrit dans le cadre du développement de matériaux composites nanostructurés à propriétés électriques et magnétiques inédites. Afin de répondre à des besoins technologiques pour, l'électronique de puissance intégrée : le stockage ou la transmission de l'énergie, les télécommunications (antenne intégrée...), le stockage de l'information par enregistrement magnétique et le marquage biologique, le composite doit présenter globalement une polarisation magnétique élevée ainsi qu'un comportement isolant permettant de pousser les limites fréquentielles, minimiser les pertes dynamiques et découpler les grains entre eux. Les matériaux composites élaborés sont constitués d'une matrice d'accueil - magnétique (ferrite spinelle) ou non-magnétique (diélectrique = silice) - dans laquelle sont dispersées des particules métalliques (Fer-Nickel ou Cobalt). Ces matériaux sont novateurs dans la mesure où le matériau final peut bénéficier d'un couplage des propriétés magnétiques des deux phases constitutives. L'holographie électronique en transmission a mis en évidence une ocnfiguration de spins de type « vortex » dans les nanoparticules de Fe30Ni70. Les mesures holographiques ont été comparées au profil de l'aimantation, dans un vortex, modélisé par une approche micromagnétique. Des analyses physico-chimiques approfondies nous ont permis de confirmer les topologies visées : pour le composite métal-diélectrique, l'épaisseur de la couche d'enrobage a pu être contrôlée à l'échelle nanonométrique. Pour le composite métal-oxyde obtenu par croissance directe du ferrite sur la phase métallique, on a démontré une bonne dispersion des particules métalliques. Les propriétés magnétiques et structurales des différents composites, en poudre ou compactés par SPS (compactage-frittage flash), ont été caractérisées et discutées. Les propriétés fonctionnelles ont été aussi étudiées et sont très prometteuses pour les applications visées. L'enrobage des nanoparticules par la silice a permis la préparation de leur surface dans la perspective d'une fonctionnalisation par des entités biologiques.

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. Une réflexion a été mise en oeuvre en vue de renforcer de manière faible la condition cinématique sur la normale dans la zone de contact via les techniques d'optimisation sans contraintes. La nouvelle formulation est dite de " Lagrangien augmenté non standard ". Trois nouvelles stratégies basées sur le Lagrangien augmenté ont été implémentées. Il s'agit de la méthode de Newton généralisée : c'est une méthode d'optimisation qui permet de résoudre toutes les non linéarités du problème en une seule boucle d'itérations. La méthode de Newton partielle est une méthode hybride entre la méthode historique du code appelée méthode de point fixe et la méthode de Newton généralisée. Enfin, on a mis en place une façon originale de traiter le cyclage : ce phénomène apparaît souvent au cours de la résolution du problème et il entraîne la perte de convergence des algorithmes. La stratégie nouvelle permet d'améliorer la robustesse des algorithmes.

Les substrats Silicium-sur-Isolant (SOI) représentent la meilleure solution pour obtenir des dispositifs microélectroniques ayant de hautes performances. Des méthodes de caractérisation électrique sont nécessaires pour contrôler la qualité SOI avant la réalisation complète de transistors. La configuration classique utilisée pour les mesures du SOI est le pseudo-MOFSET. Dans cette thèse, nous nous concentrons sur l'amélioration des techniques autour du Ψ-MOFSET, pour la caractérisation des plaques SOI et III-V. Le protocole expérimental de mesures statiques ID-VG a été amélioré par l'utilisation d'un contact par le vide en face arrière, permettant ainsi d'augmenter la stabilité des mesures. De plus, il a été prouvé que ce contact est essentiel pour obtenir des valeurs correctes de capacité avec les méthodes split-CV et quasi-statique. L'extraction des valeurs de Dit avec split-CV a été explorée, et un model physique nous a permis de démontrer que ceci n'est pas possible pour des échantillons SOI typiquement utilisés, à cause de la constante de temps reliée à la formation du canal. Cette limitation a été résolue un effectuant des mesures de capacité quasi-statique (QSCV). La signature des Dit a été mise en évidence expérimentalement et expliquée physiquement. Dans le cas d'échantillons passivés, les mesures QSCV sont plus sensibles à l'interface silicium-BOX. Pour les échantillons non passivés, un grand pic dû à des défauts d'interface apparait pour des valeurs d'énergie bien identifiées et correspondant aux défauts à l'interface film de silicium-oxyde natif. Nous présentons des mesures de bruit à basses fréquences, ainsi qu'un model physique démontrant que le signal émerge de régions localisées autour des contacts source et drain. / Silicon-on-insulator (SOI) substrates represent the best solution to achieve high performance devices. Electrical characterization methods are required to monitor the material quality before full transistor fabrication. The classical configuration used for SOI measurements is the pseudo-MOSFET. In this thesis, we focused on the enrichment of techniques in Ψ-MOSFET for the characterization of bare SOI and III-V wafers. The experimental setup for static ID-VG was improved using a vacuum contact for the back gate, increasing the measurement stability. Furthermore, this contact proved to be critical for achieving correct capacitance values with split-CV and quasi-static techniques (QSCV). We addressed the possibility to extract Dit values from split-CV and we demonstrated by modeling that it is impossible in typical sized SOI samples because of the time constant associated to the channel formation. The limitation was solved performing QSCV measurements. Dit signature was experimentally evidenced and physically described. Several SOI structures (thick and ultra-thin silicon films and BOX) were characterized. In case of passivated samples, the QSCV is mostly sensitive to the silicon film-BOX interface. In non-passivated wafers, a large defect related peak appears at constant energy value, independently of the film thickness; it is associated to the native oxide present on the silicon surface. For low-frequency noise measurements, a physical model proved that the signal arises from localized regions surrounding the source and drain contacts.

Cette étude a pour objectifs la caractérisation expérimentale des frottements saccadés dans les matériaux granulaires modèles constitués des billes de verre monodisperses en compression triaxiale drainée et l'identification des paramètres de contrôle. Cinq paramètres macroscopiques caractérisent ces frottements saccadés : la chute de déviateur et la contraction volumique, l'intermittence de déformation axiale, le module d'Young et le coefficient de Poisson. Les frottements saccadés affectent simultanément le déviateur et la déformation volumique. Le comportement macroscopique est globalement contractant tandis que le matériau tend vers un état limite critique en grandes déformations, à la manière des sables lâches. Cependant il présente localement, dans les phases de blocage qui suivent immédiatement les ruptures temporaires, le comportement dilatant des sables denses, qui obéit à une relation contrainte-dilatance linéaire et unique. Les frottements saccadés disparaissent au-delà d'une vitesse critique d'écrasement axial, qui dépend de la contrainte de confinement et de la taille des grains. Le module d'Young dynamique par propagation d'ondes varie avec la contrainte de confinement selon une loi de puissance. Le module d'Young quasi-élastique au départ des phases de blocage est constant à l'intérieur du domaine élastique, de même que le coefficient de Poisson. Le suivi par granulométrie laser et par analyse d'images des matériaux après un ou plusieurs essais triaxiaux permet de suivre l'évolution de la taille moyenne et de la forme des grains. Tandis que les instabilités par saccade disparaissent suite à un certain nombre d'essais, on observe, simultanément à une légère diminution du volume moyen, l'apparition progressive de populations d'objets non sphériques par une fusion des grains analogue au frittage. / The objectives of this study are the experimental characterisatino of the stick-slip instabilities in a model granular material and the identification of relevant control parameters. As monodisperse glass beads are subjected to drained triaxial compression tests, five macroscopic parameters characterize the stick-slip phenomenon: the deviator drop, the jump in volumetric contraction, the intermittence of the axial strain, Young’s modulus and Poisson’s ratio. The stick-slip events simultaneously affect the deviator stress and the volumetric strain. While the global material behavior is that of a loose sand, gradually contracting and hardening as it approaches its large strain critical state, its response in the “stick” phases immediately following the “slip” instabilities is similar to that of dense, dilatant sands, with a unique, linear stress-dilatancy relationship. Stick-slip events disappear beyond a critical axial strain rate, depending on the confining stress and on the grain diameter. The Young modulus associated to wave propagation varies with the confining stress according to a power law. The quasi-elastic modulus measured at the beginning of the stick phase is constant inside the elastic domain, as well as the Poisson ratio. The evolution of grain size and shape after one or several triaxial tests is monitored by laser granulometry and image analysis. The gradual vanishing of stick-slip events, on repeating the tests, is likely related to the global decreasing trend of average particle volume and to the formation of non-spherical objects, apparently by some phenomenon analogous to sintering.

Dans le cadre de la protection du bâti face au risque sismique, les techniques d’analyse simplifiées, basées sur des calculs quasi-statiques en poussée progressive, se sont fortement développées au cours des deux dernières décennies. Le travail de thèse a pour objectif d’optimiser une stratégie d’analyse simplifiée proposée par Chopra et al. (2001) et adoptée par les normes américaines FEMA 273. Il s’agit d’une analyse modale non linéaire découplée, dénommée par les auteurs UMRHA qui se caractérisent principalement par : des calculs de type pushover selon les modes de vibration dominants de la structure, la création de modèles à un degré de liberté non linéaire à partir des courbes de pushover, puis le calcul de la réponse temporelle de la structure en recombinant les réponses temporelles associées à chaque mode de vibration. Dans ce travail, la méthode UMRHA a été améliorée en investiguant les points suivants. Tout d’abord, plusieurs modèles à un degré de liberté non linéaire déduits des courbes de pushover modal sont proposés afin d’enrichir la méthode UMRHA originelle qui emploie un simple modèle élasto-plastique : autres modèles élasto-plastiques avec des courbes enveloppes différentes, le modèle de Takeda prenant en compte un comportement hystérétique propre aux structures sous séismes, et enfin, un modèle simplifié basé sur la dégradation de fréquence en fonction d’un indicateur de dommage. Ce dernier modèle à un degré de liberté privilégie la vision de la chute de fréquence au cours du processus d’endommagement de la structure par rapport à une description réaliste des boucles d’hystérésis. La réponse totale de la structure est obtenue en sommant les contributions non linéaires des modes dominants aux contributions linéaires des modes non dominants. Enfin, la dégradation des déformées modales, due à l’endommagement subi par la structure au cours de la sollicitation sismique, est prise en compte dans la méthode M-UMRHA proposée dans ce travail, en généralisant le concept précédent de dégradation des fréquences modales en fonction d’un indicateur de dommage : la déformée modale devient elle-aussi dépendante d’un indicateur de dommage, le déplacement maximum en tête de l’ouvrage ; l’évolution de la déformée modale en fonction de cet indicateur est directement identifiée à partir des calculs de pushover modal. La pertinence de la nouvelle méthode M-UMRHA est investiguée pour plusieurs types de structures, en adoptant des modélisations éprouvées dans le cadre de la simulation des structures sous séismes : portique en béton armé modélisé par des éléments multifibres pour le béton et les armatures, remplissage en maçonnerie avec des éléments barres diagonales résistant uniquement en compression, bâti existant contreventé (Hôtel de Ville de Grenoble) avec des approches coques multicouches. Les résultats obtenus par la méthode simplifiée proposée sont comparés aux résultats de référence issus de l'analyse temporelle non linéaire dynamique. / In the context of building’s protection against seismic risk, simplified analysis techniques, based on quasi-static analysis of pushover, have strongly developed over the past two decades. The thesis aims to optimize a simplified method proposed by Chopra and Goel in 2001 and adopted by American standards FEMA 273. This method is a nonlinear decoupled modal analysis, called by the authors UMRHA (Uncoupled Modal for Response History Analysis) which is mainly characterized by: pushover modal analysis according to the dominant modes of vibration of the structure, setting up nonlinear single degree of freedom systems drawn from modal pushover curves, then determining the history response of the structure by combining of the temporal responses associated with each mode of vibration. The decoupling of nonlinear history responses associated with each mode is the strong assumption of the method UMRHA. In this study, the UMRHA method has been improved by investigating the following points. First of all, several nonlinear single degree of freedom systems drawn from modal pushover curves are proposed to enrich the original UMRHA method, in which a simple elastic-plastic model is used, other elastic-plastic models with different envelope curves, Takeda model taking into account an hysteretic behavior characteristic of structures under earthquakes, and finally, a simplified model based on the frequency degradation as a function of a damage index. The latter nonlinear single degree of freedom model privileges the view of the frequency degradation during the structure damage process relative to a realistic description of hysteresis loops. The total response of the structure is obtained by summing the contributions of the non linear dominant modes to those of linear non dominant modes. Finally, the degradation of the modal shapes due to the structure damage during the seismic loading is taken into account in the new simplified method M-UMRHA (Modified UMRHA) proposed in this study. By generalizing the previous model of frequency degradation as a function of a damage index: the modal shape becomes itself also dependent on a damage index, the maximum displacement at the top of the structure; the evolution of the modal shape as a function of this index is directly obtained from the modal pushover analysis. The pertinence of the new method M-UMRHA is investigated for several types of structures, by adopting tested models of structures simulation under earthquakes: reinforced concrete frame modeled by multifibre elements with uniaxial laws under cyclic loading for concrete and steel, infill masonry wall with diagonal bars elements resistant only in compression, existing building (Grenoble City Hall) with multilayer shell elements and nonlinear biaxial laws based on the concept of smeared and fixed cracks. The obtained results by the proposed simplified method are compared to the reference results derived from the nonlinear response history analysis.