Notre travail s'inscrit dans le cadre de l'étude du comportement<br />des milieux granulaires. Dans un premier temps, nous présentons<br />une méthode numérique par Eléments Discrets de type Dynamique des<br />Contacts, qui modélise en 2D le mouvement d'un ensemble de corps<br />rigides, entrant en collision entre eux et avec des parois, et<br />sujets à des forces de frottement lors de ces chocs. L'utilisation<br />du bipotentiel de contact conduit à un algorithme local basé sur<br />un schéma prédicteur-correcteur par projection sur le cône de<br />frottement et à un critère de convergence fondé sur un indicateur<br />d'erreur relative en loi de comportement. La prise en compte<br />exacte des conditions de Signorini et de Coulomb nous oblige à<br />considérer les phénomènes de chocs multiples. Dans ce cadre, nous<br />utilisons le formalisme de la mécanique non-régulière. Nous<br />aboutissons à un algorithme comportant, à chaque itération, une<br />phase de résolution de l'équation de la dynamique, fournissant une<br />nouvelle approximation de la vitesse, et une phase d'utilisation<br />de l'algorithme local, fournissant une nouvelle valeur de<br />l'impulsion. Les simulations numériques tant quasi-statiques que<br />dynamiques mettent en évidence la convergence et la robustesse de<br />l'algorithme. <br />Dans un second temps, pour permettre le passage de l'échelle<br />microscopique à l'échelle macroscopique, nous établissons une<br />expression du tenseur des contraintes moyen, qui prend en compte<br />les efforts volumiques. Nous montrons sur un exemple analytique<br />simple d'un grain rigide cylindrique roulant sur un plan incliné,<br />que les effets dynamiques sont essentiels pour symétriser le<br />tenseur des contraintes moyen.

ETUDE DU COMPORTEMENT D'UNE STRUCTURE MULTICOUCHES EN PETITES PERTURBATIONS, ICI LA POUTRE MULTICOUCHES : IL S'AGIT D'EN PROPOSER LA CONSTRUCTION D'UN MODELE SIMPLIFIE REPOSANT SUR UNE METHODE SYSTEMATIQUE, ET DE PROFITER DES HYPOTHESES SIMPLIFICATRICES REALISEES POUR FORCER LE MODELE A RESTITUER DES INFORMATIONS SUR LES PHENOMENES D'INTERFACES - DANS LA 1ERE PARTIE, LE 1ER CHAPITRE EXPOSE LE POINT DE DEPART, LA MODELISATION MULTIPARTICULAIRE - LE 2D PRESENTE 2 PRINCIPES VARIATIONNELS MIXTES UTILISANT CHACUN UN JEU PROPRE DE VARIABLES INDEPENDANTES DUALES CONTRAINTES-DEFORMATIONS - LE 3EME PRESENTE UNE ETUDE LAGRANGIENNE DE POUTRES DANS LAQUELLE DES CONCEPTS D'EFFORTS GENERALISES CLASSIQUEMENT DEFINIS EN RESISTANCE DES MATERIAUX SONT REECRITS DANS LE CAS DES GRANDES TRANSFORMATIONS SUR LA CONFIGURATION DE REFERENCE - LA 2DE PARTIE EST CONSACREE A LA CONSTRUCTION DU MODELE : PRESENTATION DES HYPOTHESES, EN PARTICULIER DU CHOIX DU PRINCIPE VARIATIONNEL MIXTE QUI SERA UTILISE AINSI QUE DE LA NATURE DES APPROXIMATIONS QUI SERONT INTRODUITES DANS CE PRINCIPE.

Le Golfe de Californie (GOC) est un exemple de rift très oblique, au stade d'accrétion depuis 3,6 Ma au sud. La déformation continentale débute au Miocène dans un contexte arrière-arc, en relation avec la subduction de la plaque Pacifique sous la Plaque Nord-Américaine. L'objectif de cette thèse est de déterminer I'histoire de la déformation à I'origine de la rupture continentale ainsi que les phénomènes post-rupture. Cette étude utilise des données tectoniques, sismologiques et géomorphologiques sur la marge Ouest du GOC en Basse Californie du Sud. Elles sont calées temporellement par des datations d'isotopes radiogéniques et cosmogéniques. Une histoire polyphasée de la déformation, essentiellement post-magmatique, est proposée en termes d'évolution des directions des failles et de déformations relativement à la cinématique des plaques. Je démontre dans cette zone que la déformation qui conduit à l'amincissement et à l'étirement lithosphérique est d'abord fortement oblique et devient transtensive et moins oblique relativement à la cinématique à partir de 7-8 Ma. La marge étudiée est intégrée de manière cohérente à I'ensemble des domaines déformés associés à la formation du GOC. On montre notamment que les marges du GOC sont diachrones et qu'elles se forment par migration vers I'Ouest de la déformation vers la zone de subduction, qui devient inactive vers 12 Ma. La dynamique post-breakup du GOC est interprétée en termes de mobilisation par fluage de la croûte inférieure en relation avec un évènement thermique lié à la rupture du slab. En conclusion, l'évolution du GOC ne peut s'expliquer par les modèles simples de rifting oblique et de formation de marges passives. / The Gulf of California (GOC) is an example of highly-oblique rift. Oceanic accretion started 3.6 Ma ago at its southern end. The earliest continental extension started during the Miocene, in a back-arc setting, in connection with the subduction of the Pacific plate (PAC) beneath North America (NAM). ln this study, I reconstitute the strain evolution along the proximal Baja California margin. For this, I used original tectonic and seismological data which I collected in Baja California Sur (BCS). Those data were time-constrained with absolute dating (radiogenic and cosmogenic isotopes). I outline that the main stretching and thinning of the Late Miocene-Pliocene Baja California margin was highly oblique regarding NAM-PAC kinematic vector, turning less oblique. By integrating the studied margin in the evolution of the GOC, it is proposed that the final break-up mechanism occurred within a broad semi-ductile right-lateral central shear-zone. lt is shown that the two GOC passive margins didn't form at the same time, crustal strain migrating westward during the Miocene in response to a probable retreat of the dying slab. The active post-breakup deformation in the proximal BCS is best interpreted in connection with an outward flow of the lower crust in a trend compatible with the margin shaping inherited from the major GOC-normal to GOC-oblique Miocene crustal deformation. This ductile flow would be enhanced by the heat input from the slab rupture. To conclude, the GOC evolution as an oblique rift system adjusts with no existing analogical or numerical model of strain field evolution of oblique rifting, a probable consequence to the complex dynamics in back-arc settings.

Ce travail est consacré à l'étude expérimentale des effets de la pression de radiation d'une onde laser continue sur une interface liquide. Les propriétés particulières du ménisque séparant deux phases liquides en coexistence au voisinage de leur point critique de démixtion nous ont permis de visualiser directement des déformations stationnaires d'interface de taille micrométrique. Au stade linéaire, i.e pour des intensités laser modérées, une loi d'échelle exprimant la hauteur des déformations a été validée, ceci pour les deux sens de propagation du faisceau relativement à l'interface. Pour des excitations laser plus élevées, une brisure de symétrie vis à vis du sens de propagation a été mise en évidence. On observe en effet la formation de doigts de grand rapport d'aspect, ou la brisure de l'interface suite à une instabilité optohydrodynamique, suivant que le faisceau se propage du milieu le moins réfringent au milieu le plus réfringent, ou inversement. Ces caractéristiques ont été exploitées pour créer et stabiliser sous champ laser des ponts liquides de rapports d'aspect bien supérieurs à la limite de l'instabilité de Rayleigh-Plateau des colonnes liquides. Les déformations thermocapillaires, engendrées par la faible élévation de température induite par le faisceau laser, ont également été caractérisées, afin de les distinguer sans ambiguité des déformations engendrées par la pression de radiation seule.

La sismicité induite par le gisement de gaz de Lacq a débuté en 1969, dans une région tectoniquement considérée comme asismique auparavant (M<3). L'exploitation est située au Sud du bassin Aquitain, à 30 km ad Nord de la chaîne pyrénéenne. Un réseau local a été installé en 1974 afin de suivre l'évolution spatio-temporelle des séismes *Plus de 1000 événements (1.0 <=MS <=4.2) ont été enregistrés entre 1974 et 1992. La très bonne connaissance de la structure géomécanique qui piège le réservoir (forages, profils sismiques, carottages ... ) a permis de construire un modèle de vitesse 3D permettant de relocaliser de manière fine la sismicité sur la période 1976- 1989. La position des hypocentres est étroitement corrélée avec la structure en dôme du réservoir ainsi qu'avec les discontinuités préexistantes. La construction des mécanismes au foyer (133 événements sélectionnés sur la période 1976-1989) montre une grande diversité dans les mécanismes de rupture. Cependant, la méthode des dièdres droits a permis de voir une évolution temporelle globale à l'échelle du gisement. Dans un premier temps (1976-1986) nous avons eu un régime en compression maximale verticale, suivi dans la seconde période (1987-1989) par un régime en compression maximale horizontale. L'inversion du tenseur des contraintes a confirmé ces résultats. La sismicité induite de Lacq étant très localisée (du fait de l'extension limitée du gisement), cela nous a permis de considérer le gisement comme un objet géologique isolé, à une échelle intermédiaire entre les processus tectoniques et les expériences de laboratoire. Nous n'avons pas retrouvé les résultats expérimentaux obtenus sur le coefficient b (relation fréquence-magnitude des séismes) car la présence d'événements asismiques (non décelables) sur Lacq semble jouer un rôle majeur dans la libération des contraintes. La géométrie fractale nous a permis de distinguer deux classes de séismes: (1) des séismes diffus répartis sur tout le gisement et (2) des nids de sismicité corrélés avec la structure en dôme du réservoir. Bien qu'il n'y ait pas de corrélation temporelle entre D et b, nous proposons cependant qu'il existe un lien géomécanique entre ces deux paramètres.

Cette thèse porte sur les membranes biologiques, étudiées du point de vue de la physique théorique. Nous nous intéressons aux effets génériques de la présence d'une ou deux inclusions membranaires, par exemple des protéines, et à ceux d'une modification chimique locale de l'environnement de la membrane.<br /> Tout d'abord, nous étudions une interaction entre deux inclusions membranaires, qui est analogue à la force de Casimir : elle provient des contraintes que les inclusions imposent aux fluctuations thermiques de la forme de la membrane. Nous calculons les fluctuations de cette force entre deux inclusions ponctuelles. Nous clarifions la définition de la force exercée par un fluide corrélé sur une inclusion, dans un micro-état du fluide. Cette définition joue un rôle clé dans l'étude de la force de Casimir au-delà de sa valeur moyenne à l'équilibre. Nous étudions également les interactions de Casimir entre des inclusions membranaires de forme allongée.<br /> Ensuite, nous nous intéressons à l'élasticité membranaire à l'échelle nanométrique, qui est mise en jeu dans les déformations locales de l'épaisseur de la membrane à proximité de certaines protéines. Nous soulignons l'importance d'un terme énergétique qui a été négligé jusqu'à présent.<br /> Enfin, nous présentons une description théorique, développée à partir de principes fondamentaux, de la dynamique d'une membrane soumise à une perturbation chimique locale de son environnement. Nous comparons nos prévisions théoriques à de nouveaux résultats expérimentaux portant sur la déformation dynamique d'une membrane biomimétique soumise à une augmentation locale de pH.

Dans cette thèse nous développons des estimations d'erreur a posteriori par équilibrage de flux pour la poro-mécanique et la poro-plasticité.En se basant sur ces estimations, nous proposons des algorithmes adaptatifs pour la résolution numérique de problèmes en mécanique des sols.Le premier chapitre traite des problèmes en poro-élasticité linéaire.Nous obtenons une borne garantie sur l'erreur en utilisant des reconstructions équilibrées et $H({rm div})$-conformes de la vitesse de Darcy et du tenseur de contraintes mécaniques.Nous appliquons cette estimation dans un algorithme adaptif pour équilibrer les composantes de l'erreur provenant de la discrétisation en espace et en temps pour des simulations en deux dimensions.La contribution principale du chapitre porte sur la reconstruction symétrique du tenseur de contraintes.Dans le deuxième chapitre nous proposons une deuxième technique de reconstruction du tenseur de contraintes dans le cadre de l'élasticité nonlinéaire.En imposant la symétrie faiblement, cette technique améliore les temps de calcul et facilite l'implémentation.Nous démontrons l'éfficacité locale et globale des estimateurs obtenus avec cette reconstruction pour une grande classe de lois en hyperélasticité.En ajoutant un estimateur de l'erreur de linéarisation, nous introduisons des critères d'arrêt adaptatifs pour le solveur de linéarisation.Le troisième chapitre est consacré à l'application industrielle des résultats obtenus. Nous appliquons un algorithme adaptatif à des problèmes poro-mécaniques en trois dimensions avec des lois de comportement mécanique élasto-plastiques. / In this Ph.D. thesis we develop equilibrated flux a posteriori error estimates for poro-mechanical and poro-plasticity problems.Based on these estimations we propose adaptive algorithms for the numerical solution of problems in soil mechanics.The first chapter deals with linear poro-elasticity problems.Using equilibrated $H({rm div})$-conforming flux reconstructions of the Darcy velocity and the mechanical stress tensor, we obtain a guaranteed upper bound on the error.We apply this estimate in an adaptive algorithm balancing the space and time discretisation error components in simulations in two space dimensions.The main contribution of this chapter is the symmetric reconstruction of the stress tensor.In the second chapter we propose another reconstruction technique for the stress tensor, while considering nonlinear elasticity problems.By imposing the symmetry of the tensor only weakly, we reduce computation time and simplify the implementation.We prove that the estimate obtained using this stress reconstuction is locally and globally efficient for a wide range of hyperelasticity problems.We add a linearization error estimator, enabling us to introduce adaptive stopping criteria for the linearization solver.The third chapter adresses the industrial application of the obtained results.We apply an adaptive algorithm to three-dimensional poro-mechanical problems involving elasto-plastic mechanical behavior laws.