L'évolution actuelle des exigences en termes de performance énergétique des bâtiments a fait apparaître de nouveaux enjeux et problématiques scientifiques, dont ceux liés à l'humidité. Cette étude s'appuie sur une cellule expérimentale construite sur la technologie des maisons à ossature bois et soumise aux conditions climatiques réelles de Grenoble. L'instrumentation de ce bâtiment et le suivi de l'évolution en température et en humidité dans les différentes couches de l'enveloppe permettent de définir des séquences nécessaires à la validation de modèles numériques. Dans cet objectif, un modèle existant nommé HAM-Tools a été utilisé pour simuler les transferts couplés de chaleur, d'air et d'humidité à l'échelle du bâtiment. La démarche de validation a été décomposée en plusieurs étapes, de manière à cibler des transferts spécifiques et d'en améliorer la modélisation. Ces études localisées concernent les transferts couplés de chaleur et de masse à travers les parois solides, la modélisation des transferts de chaleur à travers une lame d'air ventilée et enfin la modélisation du renouvellement de l'air intérieur en conditions naturelles. Pour estimer la précision globale du modèle, c'est-à-dire à l'échelle du bâtiment, une séquence expérimentale a été simulée en prenant en compte l'ensemble des transferts couplés simultanément. Les performances du modèle sont discutées à partir des mesures locales, c'est-à-dire dans les parois, puis globales. La bonne concordance entre mesures et résultats de simulation permet de conclure sur la validité et la généricité de la démarche mise en œuvre et les hypothèses de simulation. Plus particulièrement, il est apparu que l'outil de modélisation permet de prédire correctement le comportement moyen des parois en humidité et en température. Il est donc envisageable de l'utiliser pour simuler et estimer l'impact des constituants des parois en termes de durabilité, de performances énergétiques et de confort de l'occupant.

This study is devoted to the thermo-mechanical constitutive modeling for saturated stiff clays and the development of a corresponding efficient stress integration algorithm. The mechanical behavior of natural Boom Clay in isothermal conditions was first characterized. The Modified Cam Clay model (MCC) was then applied to simulate the natural Boom Clay behavior. It has been found that the MCC gives poor-quality predictions of the natural Boom Clay behavior. Thereby, an adapted Cam Clay model (ACC-2) was developed by introducing a new yield surface and a new plastic potential as well as a Two-surface plastic mechanism. This model allows satisfactory prediction of the main features of the mechanical behavior of natural Boom Clay. Moreover, the constitutive equations of this model can be formulated mathematically as in a classic elasto-plastic model. Thus, the classic stress integration algorithm can be applied. The thermal effects were considered by assessing the performance of some advanced thermo-mechanical models (Cui et al., 2000; Abuel-Naga et al., 2007; Laloui and François, 2008; 2009). It appears that all the three models can capture the main features of the thermo-mechanical behavior of saturated clays. However, each constitutive model has its own limitations or unclear points from the theoretical point of view. The stress integration algorithm of the thermo-mechanical model proposed by Cui et al. (2000) at the stress point level was also developed using a specifically designed adaptive time-stepping scheme. The computation time required to achieve a given accuracy is largely reduced with the adaptive sub-stepping considered for both mechanical and thermal loadings. A Two-surface thermo-mechanical model (TEAM model) was developed based on the Two-surface plastic mechanism. The proposed model extends the model of Cui et al. (2000) to a Two-surface formulation, considering the plastic strain coupling between the thermal and the mechanical loading paths. The simulation of drained tests shows that this model can capture the main thermo-mechanical features of natural Boom Clay along different loading paths. The TEAM model was finally extended to undrained conditions. After setting up an appropriate effective stress principle and defining a volumetric strain condition, the undrained heating process was analyzed. The validity of the thermo-hydro-mechanical constitutive equations was examined based on the data from typical tests

Certains besoins économiques, sociétaux ou environnementaux nécessitent une caractérisation approfondie des structures géologiques du sous-sol. Parmi les différentes méthodes géophysiques, la sismique est sans doute celle qui est le plus souvent privilégiée. Imager et caractériser le sous-sol à partir de données réelles acquises sur le terrain est cependant un problème inverse généralement di fficile à résoudre. Il l'est d'autant plus si le problème direct associé n'a pas été au préalable "maîtrisé" a minima, c'est-à-dire si l'on ne dispose pas d'une modélisation réaliste de la propagation des ondes. Pour cela, il semble essentiel, pour modéliser les données sismiques synthétiques a fin de les rendre comparables aux données réelles, d'identi fier et de comprendre en premier lieu les facteurs géologiques et/ou physiques de premier ordre qui in fluencent la propagation des ondes dans une con figuration donnée. Cela suppose qu'on prenne en compte à la fois les lois de comportement des couches géologiques traversées et les lois de contact entre ces couches, mais aussi les limitations introduites par le système d'observation constitué par les ondes sismiques. Ce mémoire décrit les travaux que j'ai menés essentiellement sur la compréhension et la modélisation de l'interaction ondes/interfaces en me basant sur des approches pluridisciplinaires. Outre l'aspect fondamental d'une meilleure connaissance des mécanismes physiques sous-jacents à la propagation des ondes dans des milieux complexes (e.g., les résolutions sismiques verticale et latérale, la dé finition du ré flecteur sismique 3D, la diffusion d'interface et l'in fluence de l'anisotropie sur les ondes de surface), des aspects plus appliqués, dédiés notamment à l'exploration de gisements énergétiques ou à la surveillance de réservoirs géologiques de stockage du CO2, ont été également développés.

Les propriétés mécaniques des tissus biologiques sont des paramètres importants en médecine : ce sont des biomarqueurs du fonctionnement normal ou pathologique d'un tissu. En effet, ces propriétés peuvent être affectées par certaines conditions mécaniques telles que l'application d'une contrainte externe, comme l'hypertension ou un traumatisme, mais également par la présence de certaines maladies, telles que le cancer, la fibrose, l'inflammation, la maladie d'Alzheimer, ou bien tout simplement avec l'âge. La palpation réalisée par le médecin permet de discerner ces changements mais ce geste est qualitatif et ne peut accéder à des organes profonds. L'élastographie-IRM reste une méthode quantitative, robuste, d'une grande précision, qui permet de sonder l'élasticité et la viscosité des tissus. Elle consiste à mesurer le champ de déplacement d'une onde de cisaillement induite dans l'organe ciblé par une technique IRM en contraste de phase. Les modules viscoélastiques sont alors déduits après inversion de l'équation d'onde. Malgré cela, la justesse de cette technique n'a pas encore été pleinement établie. L'élastographie-IRM est en cours d'implémentation en routine clinique sur des patients atteints de maladies hépatiques chroniques ou bien pour caractériser des tumeurs dans le cas de cancer du sein. L'application aux autres organes protégés, tels que le cerveau ou les poumons, reste encore du domaine de la recherche à cause de la difficulté d'y induire des ondes mécaniques (protection naturelle de la boîte crânienne ou de la cage thoracique). C'est dans ce contexte qu'intervient un volet de mon travail de thèse : la mise en place, la caractérisation et l'optimisation d'un système induisant des ondes mécaniques dans les organes profonds. L'approche originale suivie a été d'utiliser les voies naturelles permettant d'amener l'onde de pression aux poumons ou bien à l'encéphale, différente des approches classiques consistant à traverser les barrières protectrices. Ce générateur d'onde de pression nous a permis d'obtenir des amplitudes d'onde allant de 6 µm à 30 µm dans l'ensemble du cerveau, amplitudes suffisantes afin d'en déduire les modules viscoélastiques du cerveau entier. D'autre part, un travail important s'est attaché à la réalisation d'un schéma original de correction des mouvements du patient en élastographie-IRM. Nous avons mis en évidence comment ces mouvements peuvent entraîner une discordance des composantes du champ de déplacement, nécessitant alors d'être corrigées. La correction proposée est composée d'une première étape dont la finalité est de recaler spatialement l'ensemble des volumes acquis, puis d'une seconde étape permettant de rétablir les composantes du champ de déplacement dans la même base orthonormée. Nous avons évalué numériquement et expérimentalement le biais induit quand aucunes corrections n'étaient appliquées sur ces données ainsi que l'apport de ces deux étapes de correction. Un travail préliminaire sur l'étude de la reproductibilité des acquisitions (phase en particulier) a été nécessaire. Enfin, l'ensemble des résultats de ces deux volets nous ont permis de réaliser des acquisitions d'élastographie du cerveau complet et d'obtenir des cartes du champ de déplacement de qualité. Ainsi, nous avons pu montrer la tendance des ondes mécaniques à suivre les directions privilégiées des fibres du cerveau, résultats que nous avons commencé à confronter aux observations faites en DTI.

Au moyen d'une expérience d'excitation paramétrique d'onde de surface, nous mettons en évidence l'existence de nouveaux types d'ondes solitaires et stationnaires à la surface de l'eau. Ces ondes de grande amplitude sont très non-linéaires et l'étude théorique réalisée ne permet pas de rendre compte de la forme des vagues mais permet de comprendre l'origine du phénomène d'hystérésis observé qui est nécessaire à la compréhension des phénomènes observés. En effet, l'existence de ces ondes (dans notre configuration expérimentale) est conditionnée par la présence d'un domaine de bistabilité dans le plan amplitude d'excitation - amplitude des vagues au coeur duquel nous avons montré qu'il était possible d'avoir coexistence de deux solutions, une d'amplitude nulle et une d'amplitude non nulle. Ces expériences en géométrie Hele-Shaw ont aussi permis de mettre en évidence des ondes enveloppes qui ne sont encore décrit par aucun modèle existant. Il s'agit à notre connaissance de la première onde enveloppe stationnaire observé à la surface de l'eau. Nous mettons aussi en évidence des ondes de gravité de très grande amplitude, qui sont formées alternativement d'étoiles et de polygones. Nous montrons que la symétrie du motif (nombre de branche de l'étoile) est indépendante de la taille et de la forme du récipient vibré. Nous montrons qu'un mécanisme de couplage non-linéaire résonant à trois ondes peut expliquer cette géométrie, bien que cette possibilité fut rejetée pour des ondes purement gravitaire.

Cette thèse contribue à la simulation numérique des écoulements d'un fluide compressible modélisé par les équations de Euler et Navier-Stokes: étude d'un schéma d'ordre élévé basé sur une matrice de masse, modélisation des écoulement turbulents compressibles à très haut Reynolds, développement des estimateurs et correcteurs d'erreurs a posteriori et a priori, et adaptation de maillage anisotrope pour les fonctionnelles d'observation. Concernant la prédiction des écoulements turbulents, on s'est intéressé aux modèles hybrides de type RANS/LES comportant les nouveautés suivantes: traitement des tourbillons de grande échelle utilisant la formulation VMS (Variational Multi-Scale) et du RANS employé sur la paroi sur une distance imposée via une zone de protection conçue afin d'éviter le phénomène assez commun apelé "grid induced model depletion". Le niveau de viscosité du modèle VMS-LES est de plus controllé par un procédé de double filtre dynamique. La seconde partie concerne l'adaptation de maillage anisotrope pour mieux observer une fontionnelle d'observation. Les estimations a priori sont réalises pour le modèle des équations d'Euler et Navier-Stokes en instationnaire en 2D et 3D. A partir de ces estimations on sait définir les maillages optimaux au cours du calcul instationnaire, en fonction de l'état et de l'état adjoint. Le système d'optimalité est discrétisé et résolu à l'aide d'une méthode de point fixe instationnaire global, comportant une stratégie de stockage/recalcul pour le couplage état/ état adjoint. Des applications à la propagation d'ondes de choc et d'ondes acoustiques sont présentées.

Les robots poly-articulés industriels sont un moyen de production moins couteux que les machines outils. De part leur structure, ils sont moins rigides, mais ils disposent d'une agilité et d'une zone de travail plus importante. L'exploitation de ces avantages pour la réalisation de certaines opérations continues, comme l'usinage par exemple, fait l'objet d'une demande croissante de l'industrie manufacturière. Ces nouvelles applications des robots poly-articulés pour l'usinage nécessitent de progresser sur le front de l'amélioration de la précision statique et dynamique de ces structures. Ainsi, afin d'améliorer la précision des robots, nous avons développé dans ce travail de thèse une méthode de planification de trajectoire basée sur l'interpolation paramétrique des courbes géométriques. Cette méthode permet de maîtriser le positionnement et la cinématique de l'outil pour les applications nécessitant un suivi de profil continu et notamment pour l'usinage. Nous proposons ainsi de qualifier les différentes souplesses des robots industriels 6 axes afin de déduire une cartographie de rigidité dans l'espace de travail cartésien. Une méthode exploitant cette cartographie permettant l'optimisation de la configuration géométrique du robot pour l'usinage est présentée. Les souplesses axiales des articulations sont intégrées dans un modèle élasto-statique utilisé pour la commande. Ce modèle permet d'anticiper les déviations statiques induites par ces souplesses articulaires. Enfin, nous mettons en évidence les défauts de transmission associés aux chaînes cinématiques des axes du robot. Nous montrons que ces défauts sont à l'origine d'une erreur de position au niveau de l'organe terminal de l'ordre de quelques dixièmes de millimètre. Un protocole d'identification de ces défauts est proposé. Ces défauts sont modélisés et intégrés dans une stratégie de correction hors ligne de position.

La Décomposition Modale Empirique (EMD) est un outil de traitement de signal piloté par les données et dédié aux signaux non-stationnaires issus ou non de systèmes linéaires. L'idée de base de l'EMD est l'interpolation des extrema par des splines pour extraire de composantes oscillantes appelées modes empiriques intrinsèques (IMFs) et un résidu. Dans cette thèse, un nouvel algorithme de l'EMD est introduit où au lieu d'une interpolation rigide, un lissage est utilisé pour la construction des enveloppes supérieures et inférieures du signal à décomposer. Ce nouvel algorithme est plus robuste au bruit que l'EMD conventionnelle et réduit le nombre d'IMFs "artificielles" (sur-décomposition). En combinant le nouvel algorithme et la méthode de séparation d'énergie (ESA) basée sur l'Opérateur d'Energie de Teager-Kaiser (OETK), un nouveau schéma de démodulation des signaux AM-FM multi-composante appelé EMD-ESA est introduit. Différentes versions de l'EMD-ESA sont analysées en terme de performance. Pour l'analyse Temps-Fréquence (TF), une nouvelle formulation de la carte TF de l'EMD-ESA appelée Transformation de Teager-Huang (THT) est présentée. Cette nouvelle Représentation TF (RTF) ne présentant pas de termes d'interférences est comparée aux RTF classiques telles que le spectrogramme, le scalogramme, la distribution de Wigner-Ville Distribution (WVD), la Pseudo-WVD et la réallocation de la Pseudo-WVD. En combinant la nouvelle formulation de la THT et la transformée de Hough, une nouvelle méthode de détection des signaux multi-composante à modulation linéaire de fréquence dans le plan TF est présentée. Cette méthode de détection est appelée transformation de Teager-Huang-Hough (THHT). Les résultats de la THHT sont comparés à ceux de la transformée WVD-Hough. Finalement, l'analyse TF par THT et par des RTF classiques (WVD, spectrogramme, etc.) de signaux réels de rétrodiffusion par des coques cylindriques de dimensions et de caractéristiques physiques différentes est présentée. Les résultats obtenus montrent l'apport de la THT comme un outil TF.

L'évaluation des risques naturels et du risque sismique en particulier a pris une importance croissante pour les autorités publiques ces dernières années, entraînant une adaptation importante des outils utilisés jusqu'à présent dans l'étude de la vulnérabilité du bâti. Il est nécessaire de mieux évaluer ces risques afin de mieux y faire face. De nouvelles approches doivent être proposées pour répondre à ces objectifs dont les approches numériques font partie. Parmi les méthodes numériques développées pour tenir compte de la spécificité du bâti maçonné figurent les approches par éléments discrets comme la méthode NSCD (" Non Smooth Contact Dynamics "). Les travaux présentés dans cette thèse visent à mieux appréhender le comportement mécanique et dynamique des ouvrages maçonnés en utilisant ce type de méthode implémenté dans le code de calcul LMGC90, afin d'affiner la compréhension et la caractérisation de la vulnérabilité sismique des édifices en maçonnerie. La méthodologie que nous cherchons à mettre en place vise à traiter la vulnérabilité sismique de bâtiments spécifiques ou au contraire d'ensembles de bâtiments de même typologie.

Cette thèse s'inscrit dans le sillage des nombreuses recherches menées autour du problème d'estimation du flot optique. Elle se distingue par le fait d'avoir été menée dans un cadre pluridisciplinaire. Elle traite en effet d'un problème fondamental en traitement d'images (flot optique) auquel ont été appliqué des techniques de calcul scientifique (parallélisation) en visant une application expérimentale de la mécanique des fluides (vélocimétrie par images de particules). Le flot optique désigne le champ de vitesses apparent observé entre deux images successives d'une même scène. Son estimation est un intermédiaire essentiel dans de très nombreuses applications et nous proposons d'appliquer des techniques venant de la communauté de la vision par ordinateur à la vélocimétrie par images de particules (PIV). À la suite d'une évaluation préliminaire des approches classiques d'estimation du flot optique employant une régularisation soit locale soit globale du champ de vecteurs, nous décidons de considérer l'approche CLG de Bruhn et al. [2005b] qui propose de combiner les régularisations locale et globale. Nous étudions dans un premier temps la mise en œuvre du "warping" dans CLG puis nous le justifions dans notre approche baptisée "modified CLG" (MCLG). Nous décrivons ensuite deux méthodes de résolution. La première repose sur une approche par calcul variationnel et permet le développement d'un solveur de type Picard-SOR par bloc employant uniquement des linéarisations locales. Ce solveur est parallélisé par une technique de "color labelling". La deuxième méthode est une approche par "variable splitting" qui repose sur le fait d'introduire une seconde variable dans le modèle ce qui permet de le séparer en sous problèmes plus faciles à résoudre. Les solveurs intermédiaires employés dans cette techniques se composent d'opérations décorrélées, se qui facilite leur parallélisation. Une évaluation complète des différentes méthodes mises en œuvre au cours de ces travaux est présentée. On s'intéresse particulièrement au parallélisme des algorithmes développés, et cela sur différentes architectures (CPU multi-cœurs et GPU). On mesure également leur précison à l'aide de séquences classiques de traitement d'images mais aussi à l'aide d'images de PIV. On montre ainsi l'efficacité de la parallélisation de ces algorithmes (aussi bien sur CPU que sur GPU) et on prouve également l'interrés de considérer une approche CLG en PIV.