Plusieurs milliers de citoyens haïtiens se sont établis au Québec après la survenue du séisme ayant frappé Haïti le 12 janvier 2010. Beaucoup d’entre eux ont été parrainés au Québec par leur famille. L’immigration au Canada en provenance d’Haïti a débuté à la fin des années 1950 en se concentrant plus particulièrement à Montréal. Le contexte politique, social et économique en Haïti est souvent cité dans la littérature pour expliquer le départ des Haïtiens et leur installation au Québec où la « communauté haïtienne » est estimée à environ 120 000 personnes. Nous avons mené une recherche qualitative exploratoire afin de mieux comprendre la situation de nouveaux arrivants haïtiens établis à Montréal suite au séisme. Pour cela, nous nous sommes inspirés des notions d'expérience (Dubet, 1994) et de représentation (Hall, 1997). Entre novembre 2012 et janvier 2013, nous avons réalisé neuf entrevues dont sept avec des femmes et deux avec des hommes âgés de 23 à 60 ans. Les résultats suggèrent que les participants accordent beaucoup d'importance à la vie dans le pays d'origine avant la migration dans l'élaboration de leurs expériences. L'expérience actuelle « ici » est élaborée en comparaison à l'expérience passée « là-bas ». Les participants insistent également sur l'importance du réseau familial lors du processus migratoire, à la fois dans les démarches d'émigration et dans l'accompagnement au pays d'accueil. Des émotions douloureuses ont été exprimées face à la séparation familiale, vécue comme une des conséquences de la migration. Les participants ont tous entamé des démarches d'insertion en emploi mais ils ont fait face à des difficultés dont la plupart d'entre eux qualifient de « structurelles » et que certains qualifient de difficultés « personnelles ». Enfin, les résultats indiquent que les participants sont porteurs d'un imaginaire migratoire qui s'exprime par des aspirations et des attentes par rapport à la société d'accueil et qui rend compte d’une tension entre ces aspirations et les expériences vécues. / Thousands of Haitian citizens have settled in Quebec after the occurrence of the earthquake that struck Haiti on January 12, 2010. Many of them have been sponsored in Quebec by their families. Immigration to Canada from Haiti began during the late 1950s focusing more particularly in Montreal. The political, social and economic context in Haiti is often cited in the literature to explain the departure of Haitians and their installation in Quebec where the 'Haitian community' is estimated at about 120,000 people. We conducted an exploratory qualitative research to better understand the situation of Haitian newcomers established in Montreal following the earthquake. To do this, we inspired by notions of experience (Dubet, 1994) and representation (Hall, 1997). Between November 2012 and January 2013, we carried out nine interviews with seven women and two men aged from 23 to 60 years. The results suggest that participants place great importance to life in the country of origin prior to migration in the development of their experiences. Experience current "here" is developed in comparison to experience past 'out there '. Participants also emphasized the importance of the family network when the migration process, both in the processes of emigration and in the accompaniment to the host country. Painful emotions were expressed against the separation family, lived as one of the consequences of the migration. Participants began efforts to integrate employment but they were faced with difficulties that most of them describe as "structural" and that some "personal". Finally, the results indicate that participants are carriers of an imaginary migration which is expressed by aspirations and expectations compared to the host society and which reflects a tension between these aspirations and experiences.

Les catastrophes naturelles sont de plus en plus fréquentes et sont à l'origine de conséquences socio-économiques souvent graves et coûteuses. La prévision de ces risques, notamment en ce qui concerne les glissements de terrain, est dès lors indispensable. Dans cette problématique, la modélisation numérique doit être pertinente et prédictive. Actuellement, les glissements sous une faible pente (inférieure à 14°) voire très faible (inférieure à 8°), où la rupture est souvent de type diffuse, ne peuvent être expliqués ni par les méthodes empiriques classiques ni par l'analyse limite traditionnelle (théorie de plasticité). De nouvelles approches sont donc nécessaires pour rendre compte de ces mécanismes de glissement. La condition suffisante de stabilité de Hill (1958) basée sur le signe du travail du second ordre pourrait être un critère pertinent pour analyser ces modes de rupture diffus.<br /><br />Ce travail de thèse porte d'une part sur l'analyse et l'implémentation du critère de stabilité de Hill dans deux codes de calculs par éléments finis (LAGAMINE et PLAXIS) et d'autre part sur leur application, dans le cadre des modélisations numériques de problèmes aux limites non linéaires, aux glissements de terrain.<br /><br />L'étude du travail du second ordre, menée sur des modèles de comportement incrémentaux non linéaires, montre que certains sols exhibent de larges domaines potentiellement instables à l'intérieur strict du critère limite de plasticité. Ces analyses sont corroborées par les résultats expérimentaux.<br /><br />Nous montrons, dans ces travaux, par des approches analytique et numérique dans le cadre de la théorie de l'élasto-plasticité, que le critère de Hill est toujours vérifié avant la condition limite de plasticité (critère de Mohr-Coulomb) et le critère de localisation de déformation (critère de Rice).<br /><br />L'utilisation du critère de Hill nous permet d'analyser le glissement de la côtière de Trévoux après une période de fortes précipitations en 1983. Ces travaux de modélisations ont été menés dans le cadre d'un couplage hydromécanique en milieu non saturé.<br /><br />La même démarche a été menée pour l'analyse du glissement catastrophique du versant de Las Colinas (El Salvador) provoqué par un séisme en 2001. La définition du chargement considéré est basée sur une méthode de type pseudo-statique.

Les observations sismologiques montrent l'existence d'une phase d'initiation qui précède l'arrivée des ondes P associées à la rupture.Par ailleurs de nombreux travaux expériementaux ont permis ont permis de préciser la phénoménologie du frottement sec et décrivent également une phase de nucléation précédant l'instabilité dynamique. Nous étudions les implications du formalisme en glissement tant sur l'initiation que sur le déclenchement des séismes. La méthode proposée est fondée sur une analyse aux valeurs propres couplée à la modélisation des processus en différences finies et éléments finis. Elle permet la prise en compte des échelles de temps caractéristiques de la nucléation et de la propagation de la rupture (0.01 à 100 s). De plus cette méthode est valide pour les failles homogènes comme pour les failles hétérogènes. Nous présentons lci principaux résultats observationnels et expérimentaux, ainsi que les possibles modélisations de la nucléation. L'analyse aux valeurs propres met en lumière un domaine de stabilité de la faille finie , limité par la première valeur propre statique &0. La mécanique propre de la faille finie est contrôlée par le facteur & = a x a', où a est la demi longueur de la faille et a' l'affaiblissement de la loi de frottement. Lorsque & < &0 la faille est stable lorsque &>=&0 la faille est instable. La méthode utilisée permet de relier le taux de croissance de l'instabilité y à& lorsque la faille est instable. Nous montrons que la durée de la phase d 'initiation Tc peut être extrêmement longue lorsque & tend vers &0. La faille est un objet géologique à replacer dans son contexte géologique. La mécanique intrinsèque de la faille est alors confrontée aux différents termes de forçage. Nous nous intéressons à l' interaction entre failles, composée d 'une partie statique et d 'une partie dynamique. Nous montrons que le terme dynamique, actuellement peu étudié, peut déclencher la rupture . L'application au cas réel du séisme d'Irpinia permet de préciser les valeurs des paramètres du frottement, et de préciser les influences respectives de la dynamique et du statique sur le déclenchement. Nous étudions ensuite la faille finie hétérogène. Dans un premier temps, nous montrons l'équivalence spectrale entre un système de failles et une faille homogène équivalente associée à un affaiblissement moindre. Nous montrons l'existence d'un domaine de confiance sur lequel il est possible de mesurer l'affaiblissement effectif du système hétérogène. Cette mesure sert à contraindre le temps d'occurence du prochain événement. La modélisation en différences finies permet de définir une loi de frottement effective non linéaire qui correspond à l'homogénéisation du système hétérogène. L'introduction de l'hétérogénéité pose la question de la résistance de la faille et de l'arrêt de la rupture. Nous présentons une définition possible de la résistance fondée sur les expériences de frottement et de rupture des matériaux. La loi constitutive présente une partie durcissante avant l'affaiblissement. L'introduction du durcissement supprime les singularités de contrainte qui apparaîssent avec les barrières rigides. L'arrêt de la rupture se fait sur une zone qui correspond à l'endommagement de la barrière en durcissement. Ce modèle d 'arrêt de la rupture est à même d 'expliquer la croissance des failles au travers de la répétition des séismes .

La sismicité associée à l'activité d'un volcan intraplaque tel que le Piton de la Fournaise se développe, a priori, sans aucun contrôle tectonique régional. La com- plexité des structures internes et le fonctionnement des structures de stockage et de transfert magmatique déclenchent une sismicité partiellement récurrente. Cette récurrence de la sismicité s'observe grâce aux formes d'ondes des séismes, quasiment identiques, et témoigne de la réactivation d'une même source. Notre étude exploite la partie récurrente de la sismicité volcano-tectonique du Piton de la Fournaise entre 2000 et 2008, dans le but d'identifier les sources et de déterminer la structure interne du volcan et son fonctionnement. En utilisant les données du réseau sismologique de l'Observatoire Volcanologique du Piton de la Fournaise, nous avons regroupé par familles les séismes ayant les mêmes formes d'ondes. Puis, nous avons mesuré des différences de temps de trajet entre les séismes d'une même famille. Ces derniers sont relocalisés par la méthode de double différence et leur mécanisme au foyer est déter- miné synthétiquement en conservant l'hypothèse de l'activité d'une même source. Les déterminations et les relocalisations sont utilisées pour extraire les orientations du plan de faille et du glissement cosismique de chaque famille. Nous avons ainsi pu cartographier les zones sismogènes et observer leur comportement lors des périodes de préparation des éruptions. Ces observations ont permis de construire un modèle descriptif des phénomènes intrusifs sismogènes associés aux éruptions distales du Pi- ton de la Fournaise. Le champ de contraintes induit par chaque objet de ce modèle a été modélisé numériquement en utilisant la méthode des éléments finis pour vérifier la validité du modèle descriptif. Notre étude est limitée par la non-exhaustivité des jeux de données utilisés. Notre méthode d'analyse de la sismicité a donc été com- binée à une méthode de détection et d'estimation précise des temps d'arrivée des ondes P et S, ce qui a permis d'appliquer les apports de cette thèse à la surveillance en temps réel des volcans.

Le but de ce travail est de préciser la déformation en mer à l'est de Taïwan et de notamment caractériser les failles actives majeures en particulier au niveau de l'avant-arc des Ryukyus. Le problème majeur associé à l'étude de cette région est que les séismes sont localisés en dehors des réseaux permanents et que cela conduit à des localisations peu précises limitant les analyses statistiques, tectoniques et donc l'estimation correcte de l'aléa sismique. Ce travail est inclus dans le projet ANR ACTS (Active Tectonics and Seismic Hazard in Taiwan). Pour atteindre cet objectif, nous travaillons à trois échelles de temps différentes pour lesquelles une méthode de localisation différente a été proposée pour obtenir une image précise de la déformation sismique et mettre en évidence les failles majeures en mer à l'est de Taïwan. D'abord, à l'échelle du siècle (à partir de l'analyse d'un catalogue de sismicité homogène en magnitude), les séismes historiques instrumentaux (1897-2007) montrent que la marge des Ryukyus a été affectée par 6 séismes de magnitude supérieure à 7 dont 4 pour lesquels la faille responsable n'est pas connue. Nous proposons une méthode de localisation relative qui permet de relocaliser les évènements historiques à partir des différences de temps d'arrivées P et S (Ts-Tp) à chaque station. La localisation est obtenue en recherchant les séismes analogues dans le catalogue instrumental récent (1991-2008) pour lesquels Ts-Tp à chaque station est proche. Cette méthode a été appliquée au plus gros séisme (en terme de magnitude) jamais enregistré à Taïwan, le séisme du 5 juin 1920 (M7.7 $\pm$ 0.2). Ce séisme a certainement pour origine une rupture sur l'interplaque sismogène de la subduction des Ryukyus avec une possible nucléation à la base d'une faille hors-séquence. Ensuite, à l'échelle de la dizaine d'années, l'amélioration de la méthode de localisation absolue du MAXimum d'Intersection EDT (MAXI, EDT: Equal Differential Time) permet de mieux extraire les temps d'arrivées erronés et empêche les "trade-off" entre profondeur et temps d'origine et, entre profondeur et position épicentrale lorsque la couverture azimutale est faible (gap azimutal > 180°). Des tests synthétiques montrent l'efficacité de MAXI (en utilisant les ondes P uniquement) à déterminer les paramètres x, y et z même lorsque le gap azimutal est important. Dans ce dernier cas, la qualité des résultats reste dépendant du modèle de vitesse à représenter la structure 3D de la Terre. Nous proposons une démarche associant l'utilisation d'un modèle de vitesse des ondes P 3D a priori avec la méthode MAXI pour localiser les séismes qui sont latéraux et éloignés du réseau. Une application au cas de Taïwan et en Équateur valide cette approche. Enfin, à l'échelle de quelques mois, la campagne de géophysique marine RATS (Ryukyu Arc: Tectonics and Seismology) a été menée en deux temps. Une expérience de sismologie passive de juillet à octobre 2008 (RATS1) a été menée au-dessus de l'avant-arc des Ryukyus puis une expérience de sismique active (réfraction et réflexion) a été menée en mai 2009 sur une ligne NNE-SSW à travers la marge des Ryukyus. Ces deux expériences combinées permettent d'améliorer notre connaissance de la structure crustale de la marge. Au niveau de l'avant-arc, le socle de l'avant-arc est caractérisé par une extrémité formant un butoir vertical et une base très déformée certainement associée à une importante déformation hors-séquence. En profondeur, la plaque plongeante est certainement affectée par une déchirure qui contrôle la sismicité dans cette région de transition entre la subduction et la collision.

Simuler les mouvements du sol produits par un séisme est une étape essentielle pour évaluer l'aléa sismique d'une région. L'objectif de cette étude est de fournir un code simple d'utilisation permettant de simuler de manière efficace les mouvements du sol. Pour s'affranchir de la connaissance du milieu de propagation, nous utilisons la technique des fonctions de Green empiriques. Les mouvements du sol sont alors reconstruits par la sommation des enregistrements de plusieurs petits séismes décalés dans le temps de manière à reconstituer le processus de rupture du gros séisme, La principale problématique de cette étude est la mise en place d'un schéma de sommation efficace permettant de générer des accélérogrammes réalistes dans le domaine temporel et en bon accord avec les lois d'échelles (loi de similitude et modèle en w~2) dans le domaine fréquentiel. Grâce à l'intégration de méthodes statistiques, nous évitons les problèmes majeurs rencontrés classiquement dans ce type d'étude. Nous aboutissons finalement à deux nouvelles méthodes de simulations: une méthode purement stochastique nécessitant la connaissance de seulement deux paramètres concernant le séisme cible (moment sismique et chute de contrainte) et une méthode stochasto-cinématique permettant de reproduire certaines propriétés cinématiques de la rupture sismique sur le plan de faille. Ces deux méthodes sont testées avec succès sur des données réelles (Mexique et Italie). La méthode stochastique est ensuite appliquée pour simuler les mouvements du sol susceptibles d'être produits par deux-séismes de magnitude modérée dans la région niçoise: un séisme de magnitude 5.7 à terre (10 km au nord de Nice) et un séisme de magnitude 6.3 en mer (25 km au sud de Nice). Ces simulations mettent en évidence trois points importants: (1) les forts effets d'amplifications sur certains sites à Nice, (2) l'importance des mouvements du sol susceptibles d'être produits par un séisme de magnitude 5.7 à terre et (3) le rôle majeur joué par la valeur de la chute de contrainte du séisme potentiel.

Nous présentons une analyse des variations du taux de sismicité liées aux doublets de forts séismes. Cette étude se concentre sur la zone des répliques du 1er séisme, et cherche à estimer l'impact du 2nd choc principal sur cette activité. Nous nous sommes penchés sur deux cas de doublets de séismes : la séquence turque des séismes d'Izmit et de Düzce en 1999, et le doublet sud-islandais de juin 2000. Ces deux exemples nous apprennent que la réactivation de sismicité induite par le 2nd séisme d'un doublet n'est pas forcément corrélée avec les régions d'augmentation de la contrainte statique, et qu'elle s'établit avec prédilection sur les sites d'une activité géothermale. Sa décroissance temporelle n'est pas systématiquement conforme à la loi d'Omori, puisque des apaisements de ce taux peuvent apparaître avec un délai de quelques jours à quelques mois. Nous avons également investigué les propriétés des doublets de séismes à l'échelle globale, et nous montrons que ceux-ci peuvent être simplement décrits par les lois habituelles gouvernant la sismicité, la loi d'Omori-Utsu et la loi de Gutenberg-Richter. Enfin, contrairement aux résultats obtenus à partir de catalogues régionaux, il semblerait que la production moyenne de répliques consécutive aux 2nds séismes des doublets soit équivalente à la superposition des productions de deux séismes isolés de magnitudes similaires. Cette dernière remarque mérite toutefois d'être étayée par une analyse plus fine des catalogues régionaux, au vu du manque de résolution d'un tel catalogue global pour les séismes de petite magnitude.

A partir des sismogrammes, une phase de nucléation des séismes, peu radiative, a été détéctée. En laboratoire, des expériences de rupture, entre deux blocs rocheux, décrivent un processus de nucléation et donnent une loi d'affaiblissement de la résistance de l'interface en fonction de son glissement relatif. Nous interprétons une partie de ces observations par un modèle d'initiation: une instabilité de frottement en élastodynamique. Dans le cas unidimensionnel, alors qu' une loi en vitesse aboutit à un problème mal posé (multiplicité de solution, chocs et nécessité de la convention de retard maximal), une loi en glissement fournit une solution unique, continue, avec un temps d'instabilité non nul. Nous étudions les propriétés de l'initiation dans des problèmes de cisaillement bi et tridimensionnels en élasticité complète avec la loi en glissement. Quand une linéarisation est possible, nous menons une analyse spectrale et trouvons une approximation analytique (partie dominante). Pour des processus complexes (affaiblissement non-linéaire, hétérogénéité et propagation de la rupture), nous simulons numériquement en différences finies. Le temps d'initiation est fonction de la géométrie et de la taille de la faille, en rapport avec le taux d'affaiblissement initial. Proches de la stabilité (taux faible) les failles finies s'initient très lentement. Aussi, le signal émis en surface par l'initiation en profondeur est très dépendant de ce taux. Enfin nous montrons que l'hétérogénéité rallonge le temps d'initiation. Ces propriétés ont été vérifiées à la fois dans les cas bi et tridimensionnels.

La variation spatiale du mouvement sismique a des effets significatifs sur la réponse dynamique des structures de génie civil de grandes dimensions. Dans la pratique courante, l'excitation du mouvement sismique le long de la fondation de la structure est considérée uniforme, approche cependant inadéquate pour les structures de large portance au sol localisées à proximité des failles ou sur des sites présentant une structure du sous-sol latéralement hétérogène. Cette thèse se propose donc de comprendre les facteurs clefs contrôlant localement la variabilité spatiale du mouvement sismique, avec en ligne de mire la mise en place de recommendations en vue d'incorporer ces effets dans l'estimation de l'aléa sismique et le dimensionnement des structures. La première partie de cette thèse s'intéresse à la composante intra-évènement de l'écart-type de la distribution du mouvement sismique en champ proche à l'aide de simulations numériques du mouvement sismique pour des sources étendues présentnat une cinématique de rupture réaliste. Les résultats suggèrent que la variabilité intra-évènement dépend significativement du type de rupture, cette variabilité augmentant avec la distance pour les ruptures bilatérales et diminuant pour les ruptures unilatérales. La seconde partie traite de la caractérisation de la composition du champ d'onde dans la vallée de Koutavos-Argostoli, qui est une vallée de petite dimension et d'épaisseur sédimentaire faible, située sur l'île – sismiquement active - de Céphalonie en Grèce. Les champs d'onde générés par 46 séismes, ayant des magnitudes variant entre 2 et 5 et des distances épicentrales jusqu'à 200 km, ont été analysés à partir de l'enregistrement par deux réseaux denses de capteurs sismologiques. L'algorithme de traitement d'antenne MUSIQUE est utilisé pour extraire la vitesse, l'azimut, le type et la polarisation des ondes dominantes se propageant à travers le réseau. Les résultats montrent clairement d'importantes diffractions d'ondes de surface aux bords de vallée au-delà de la fréquence de résonance de la vallée. Tandis que les ondes de Love dominent clairement le champ d'ondes proche de la fréquence de résonance, les ondes de Rayleigh dominent à plus haute fréquence dans des gammes de fréquences étroites. Par ailleurs, un excellent accord est observé entre les champs d'onde de surface diffractés localement et les caractéristiques d'amplification du site. La “cohérence décalée” de la partie la plus énergétique du signal a été quantifiée pour chaque paire de stations du réseau. En général, la cohérence calculée sur les composantes horizontales diminue avec la distance entre stations et la fréquence. La cohérence sur la composante verticale indique des valeurs relativement fortes à haute fréquence. Les valeurs de cohérence apparaissent très faiblement corrélées à la magnitude, l'azimut et la distance épicentrale du séisme, mais sont au contraire liées aux caractéristiques géométriques de la vallée. La coherence est systématiquement plus forte pour les couples de capteurs orientés selon la direction parallèle à l'axe de la vallée, et moins forte pour des couples de capteurs orientés dans la direction perpendiculaires. Cette observation est en accord avec les résultats du traitement d'antenne: la prédominance d'ondes de surface diffractées sur les bords de la vallée conduit à des mouvements en phase le long de la direction parrallèle à l'axe de la vallée. Les résultats de cette thèse apporte des elements de compréhension sur la variabilité spatiale du movement sismique et ouvrent de nombreuses perspectives d'application. / Spatial variation of earthquake ground motion over short distances significantly affects the dynamic response of engineered structures with large dimensions. In current practices, the ground motion excitation across the foundation of a structure is assumed to be spatially uniform, which becomes inadequate for spatially extended structures in the near-fault region or on sites with lateral inhomogeneity. This PhD thesis seeks to understand the key parameters that locally control the ground motion spatial variability with the intent of putting forth practical propositions for incorporating such effects in seismic design and hazard assessment. The first part of the thesis addresses the within-event component of the standard deviation of ground-motion distribution in near source region by means of numerical simulation of ground motions for extended sources with realistic rupture kinematics. The results suggest that the within-event variability significantly depends on the rupture type, depicting an increase with distance for bilateral ruptures and a decrease for unilateral ruptures. The second part deals with the characterization of seismic wave field at the Koutavos-Argostoli site, a small-size, shallow, alluvial valley located in the seismically active Cephalonia Island in Western Greece. The seismic wave field was investigated from the recordings of a dense seismological array for a set of 46 earthquakes, with magnitude 2 to 5 and epicentral distance up to 200 km. The MUSIQUE array analysis algorithm was used to extract the phase velocity, back-azimuth, type and polarization of the dominant waves crossing the array. The results clearly indicate dominant scattering of seismic surface waves, mainly from the valley-edge directions, beyond the fundamental frequency of the valley. While Love surface waves clearly dominate the wave field close to the resonance frequency, Rayleigh waves strongly dominate only in relatively narrow frequency bands at higher frequency. Besides, an excellent consistency is observed between the dominance of the identified surface wave type in the wave field and the site amplification. The "lagged coherency" of the most energetic part of the ground motion has been quantified for each station-pair within the array. In general, spatial coherency estimated from the horizontal components exhibit decays with frequency and interstation distance. Estimates from the vertical component exhibit rather larger values at some higher frequencies. Although coherency does not show any consistent trend indicating dependence on the magnitude, back-azimuth or site-to-source distance of the event, it seems to be primarily controlled by the site geometry. Larger coherency is systematically observed when the station pair is oriented parallel to the valley axis, while lower values are observed in the perpendicular direction. This observation proves to be consistent with the MUSIQUE analysis results: the predominance of scattered surface waves propagating across the valley implies an in-phase motion along valley-parallel direction and out-of-phase motion along valley-perpendicular direction. The findings of the present research are expected to contribute in enhancing our understanding of spatial variability of ground motion and improving the coherency models used in engineering. This work also opens up new insights and many questions in need of further investigation.

Trois stratégies de modélisation de structures sont proposées dans ce document. La première consiste à l'utilisation des éléments poutres multifibres de cinématique Timoshenko et des lois de comportement basées sur la mécanique de l'endommagement pour le béton et la plasticité pour les aciers. Après une présentation succincte de plusieurs éléments poutres multifibres existants dans la littérature, un nouvel élément Timoshenko avec des fonctions d'interpolation d'ordre supérieure est proposé. La deuxième stratégie concerne la prise en compte de l'Interaction du Sol avec la Structure (ISS) et le développement d'un élément d'interface 3D modélisant une fondation rigide superficielle reposant sur un massif de sol. Basé sur le concept des macro-éléments, ce nouvel outil est formulé en variables globales (forces et déplacements) et selon la théorie de plasticité, permettant ainsi de diminuer le temps de calcul. Enfin, la troisième stratégie traite de problème de la localisation des déformations et de l'objectivité des résultats des calculs dans un milieu enrichi (milieu avec microstructure). L'évolution de la zone de localisation est étudiée avec des modèles second gradient locaux pour différents types de lois simples basées sur la théorie de plasticité ou de l'endommagement. Des comparaisons avec des résultats expérimentaux et/ou des exemples numériques prouvent la performance des stratégies de modélisation proposées.