Le disque intervertébral est principalement composé de deux types de tissu : le nucleus pulposus au centre contenu par l'annulus fibrosus en périphérie.Le nucleus est relativement bien décrit dans la littérature contrairement à l'annulus dont la constitution bi-phasique couplée à une structure à fibres orientées reste délicate à modéliser. La caractérisation hydro-chimio-mécanique représente un point clé dans la compréhension des phénomènes qui régissent le comportement de l'organe.Pour répondre au manque d'information, deux axes de travail sont mis en place : expérimental et numérique.La réalisation de tests de traction cyclique associés à des techniques de corrélation d'images permet de fournir une base de données complète du comportement.Elle est ensuite utilisée pour identifier un modèle hydro-chimio-mécanique qui traduit le comportement non-linéaire, anisotrope, dissipatif ainsi que la sensibilité aux variations d'environnement chimique du tissu observés expérimentalement.L'intérêt majeur de cette méthode est l'identification simultanée de l'ensemble des paramètres du modèle au cours de sollicitations proches des conditions physiologiques sur chacune des éprouvettes.Ce protocole réduit ainsi les incertitudes liées aux variabilités inter-individus et inter-localisations pour obtenir des valeurs de paramètres plus fiables.La formulation identifiée permettra de proposer un modèle compréhensif des transferts de fluide impliqués dans les processus de nutrition cellulaire.Ce sujet s'intègre dans une réflexion plus large relative à la compréhension des couplages hydro-chimio-mécaniques dans les tissus biologiques.L'enjeu est de quantifier l'impact des stimuli mécaniques sur l'évolution des tissus humains et leur adaptation aux conditions environnementales. / Intervertebral disc is mainly divided into two parts : nucleus pulposus surrounding by annulus fibrosus.Nucleus is relatively well described in literature, on the other hand annulus remains difficult to model due to its bi-phasic constitution coupled with an oriented fibres structure.Hydro-chemo-mechanical characterization represent a key point to reach a better understanding of this organ.To overcome this lack of knowledge, two axes are developed in this work : experimental and numerical.Cyclic tensile tests associated to digital image correlation techniques allow to obtain a large amount of data on global tissue behaviour.It permit to identify a hydro-chemo-mechanical model which takes in account the non-linear, anisotropic and dissipative behaviours.It also highlights sensibility to chemical variations in tissue environment observed experimentally.The most important point of this work is the simultaneous parameter identification performed on each sample during solicitations similar to in vivo conditions.This protocol reduce inter-individuals and inter-localisations discrepancies and leads to an accurate bundle of parameters.The formulation aims to accurately model fluid motion inside intervertebral disc tissues which is related to cell nutrition, growth and homoeostasis.This work is involved in a larger problematic about hydro-chemo-mechanical coupling within biological tissues.It will allow to quantify mechanical stimuli influence on human tissues evolution and adaptation to environment conditions.

Dans le domaine de l'ingénierie, connaitre le chargement appliqué sur une structure permet de résoudre des problèmes directs dont le résultat est le champ de déplacement, de déformation dans une structure. Il est alors possible d'effectuer un dimensionnement. Cependant, parfois ce chargement doit être identifie a posteriori. Malheureusement, il n'est pas toujours possible de mesurer ce chargement : ainsi, par exemple, on ne sait pas a priori où aura lieu le chargement, ou bien il n'est pas possible de placer un capteur sans l'endommager ou encore il peut nécessiter un encombrement trop important. On a alors recours à des mesures indirectes de déplacement, de déformation, d'accélération et on est alors amené à résoudre des problèmes inverses, qui sont en général mal posés. Il est alors nécessaire d'ajouter une (des) conditions supplémentaire(s) pour obtenir une solution unique et stable : c'est la régularisation du problème. Ces techniques sont bien connues et leur essor est dû à l'utilisation des décompositions en valeurs singulières des matrices de transfert. Toutefois, elles nécessitent l'utilisation d'un paramètre additionnel qui pondère cette condition supplémentaire : la détermination de ce paramètre est délicate. Peu de travaux ayant été réalisé pour tester de façon intensive les méthodes usuelles de régularisation (Tikhonov et troncature de la (G)SVD), en association avec les différents critères de détermination du paramètre de régularisation et les différentes réponses possibles, on a effectué un tel travail pour tirer des conclusions sur la méthodologie optimale. On a pu mettre en évidence que la mesure de l'accélération associée à un critère faisant intervenir les dérivées du signal à reconstruire donne en général les meilleurs résultats sous réserve d'utiliser le critère GCV pour déterminer le paramètre de régularisation. Ces méthodes supposent que la localisation de la zone de chargement est connue. Aussi on s'est intéressé à déduire cette zone de chargement en tentant de reconstruire des chargements identiquement nuls. Cette identification a été effectuée aisément sous réserve qu'on ait peu de forces à identifier par rapport au nombre de mesures disponibles. En revanche une telle identification est délicate lorsqu'on n'a pas plus de mesures que de forces à identifier. Finalement on s'est tourné vers l'identification de chargement ayant plastifié la structure étudiée. On a alors essayé de reconstruire le chargement en supposant que la structure reste linéaire élastique, alors qu'elle a été plastifiée : on a utilisé la méthode du double chargement et effectue des simulations à l'aide du logiciel de simulation Ls-dyna.La force reconstruite fait alors apparaitre une composante statique traduisant la déformation résiduelle dans la structure. Dans ce cas, la réponse à utiliser pour identifier le chargement est une déformation dans une zone non plastifiée / In the field of the engineering, knowing the load applied on the structure which allows to solve the direct problem of which the results are given the field of displacement and strain in a structure. It is possible to perform a dimensioning. However, sometimes this load must be identified a posteriori. Unfortunately, it is not always possible to measure this load. Thus, for example, we do not know a priori where it will be loaded, either it is not possible to place a sensor without damaging it or needs too much space. We then have to use indirect measures of displacement, strain, acceleration and then we are lead to solve the inverse problems which are generally an ill-posed. It is then necessary to add one (or more) conditions to obtain a unique and stable solution: it is the regularization of the problem. These techniques are well known and their development is due to the use of the singular value decomposition of the transfer matrix. However, they require the use of an additional parameter that weights this additional condition: the determination of this parameter is difficult. Few studies having been realized in way the usual regularization methods of (Tikhonov and truncation of the (G)SVD), in association with the various criteria for determining the regularization parameter and the various possible responses, we conducted a such work, to draw conclusions on the optimal methodology. It has been highlighted that the measurement of the acceleration associated with a criterion involving the derived signal to reconstruct generally gives the best results via the GCV criterion to determine the regularization parameter. These methods suppose that the location of the loading area is known. We also were interested to deduct this loading area while trying to reconstruct load that is identically zero. This identification was performed easily that has little load to identify compared to the number of measurements available. However such identification is difficult when there are no more measures than loads to identify. Finally we turned to the identification of loading with the plastic structure. We then tried to reconstruct the load assuming that the structure remains linear-elastic, while it was plasticized: we used the method of the double load and performed simulations using the software ls-dyna. The reconstructed load then shows a static component reflecting the residual strain in the structure. In this case, the response used to identify the load is a strain in a non-plasticized zone

La compression élastique de la jambe est le traitement de référence des pathologies liées à l’insuffisance veineuse. Bien que l’efficacité ne soit plus à prouver, les objectifs thérapeutiques restent non atteints pour certains patients. Un objectif de la compression élastique est la réduction de la pression pariétale des veines afin de rétablir ou d’augmenter le retour du sang vers le cœur par une transmission de pression au travers des tissus mous. Ce mécanisme est complexe et peut être prédit par des modèles éléments finis personnalisés. Pour être personnalisés, ces modèles doivent prendre en compte la géométrie et la carte des propriétés mécaniques du sujet.Dans cette étude, on développe deux méthodologies permettant d’identifier les propriétés mécaniques des tissus mous. Dans un premier temps, on mesure par élastographie ultrasonore par onde de cisaillement la distribution du module élastique au sein des tissus mous superficiels. Dans un deuxième temps, on identifie par une méthode inverse les propriétés mécaniques des tissus mous profonds. Cette méthode associe l’acquisition de données d’un essai expérimental de compression localisée de la jambe à un modèle éléments finis bidimensionnel. Ces deux méthodologies nous permettent d’évaluer l’hétérogénéité des propriétés mécaniques de la peau au fascia cruris et de caractériser le comportement non-linéaire des tissus mous profonds. Enfin, les résultats de ces deux méthodologies sont couplés afin de générer un modèle biomécanique de la jambe sous compression élastique pour prédire la distribution de pression au sein des tissus mous pour quatre sujets sains. / Elastic compression of the leg is a widely used treatment in case of pathologies related to venous insufficiency. Its benefits are not to be proven, but still, for some patients, the therapeutic goal is not reached. One goal of this treatment is to reduce transmural pressure applied to veins in order to restore or increase blood return to the heart by the transmission of the external pressure through soft tissues. This is a complex mechanism that can be predicted by patient-specific finite element models. To be patient-specific, these models must take into account the geometry and the distribution of mechanical properties of each subject.In this study, two methodologies are developed to identify the mechanical properties of soft tissues. First, the elastic modulus distribution inside the superficial soft tissues is measured by shear wave ultrasound elastography. Second, the mechanical properties of deep soft tissues are identified through an inverse method combining the data acquired from an experimental localized compression of the leg to a bi-dimensional finite element model.These two methodologies allow to evaluate the mechanical properties heterogeneity from the skin to the fascia cruris and to characterize the non-linear behaviour of deep soft tissues. Finally, the results from both methodologies are brought together to generate a biomechanical model of the leg under elastic compression to predict pressure distribution inside soft tissues for four healthy subjects.

La compréhension du comportement musculaire est un champ d’investigation primordial dans divers domaines tels que le sport, le crash automobile et la médecine. Peu de techniques in vivo permettent aujourd’hui de caractériser les propriétés mécaniques du muscle. Aussi, l’objectif de ce travail est de développer une méthode d’identification in vivo à partir demesures de champs de déplacement. La méthode d’identification est composée de 3 étapes interdépendantes. La première consiste à réaliser une IRM de la cuisse pour segmenter manuellement différents tissus musculaires (quadriceps, ischios, gracilis et sartorius) et le tissu adipeux. Un comportement Néo-Hookéen est choisi pour modéliser le comportement hyperélastique (C10, D). Dans un second temps, un dispositif expérimental de compression est développé pour mesurer le champ de déplacement in vivo à partir des techniques d’imagerie échographique et de Corrélation d’Images Numériques. Finalement, une méthode inverse est mise en œuvre pour identifier les paramètres C10 et D de chaque tissu. Par le biais d’un exemple numérique, l’erreur sur les paramètres identifiés est évaluée. Les cartographies des champs de déplacement expérimentaux confirment les observations qualitatives obtenues sur les images échographiques et sont validées par les champs de déformation obtenus par approximation diffuse. Une faible erreur d’identification (C10<3%, D<7%) est obtenue à partir de l’exemple numérique, et les paramètres mécaniques identifiés sont en accord avec la littérature. Ces résultats valident la démarche inverse mise en œuvre qui permettra, à terme, de suivre l’évolution des pathologies et de mener des simulations prédictives. / The understanding of the mechanical behavior of the muscle tissue is an important field of investigation with different applications in medicine, car crash and sport. Currently, few in vivo imaging techniques are able to characterize the mechanical properties of muscle. Thus, this study aims at developing an in vivo identification method based on displacement fieldmeasurements. The identification approach is composed of 3 dependent steps. The first step consists in performing a 2D MRI acquisition of the thigh in order to segment manually the muscle (quadriceps, ischio, gracilis and sartorius) and fat tissues. A Neo-Hookean model is chosen to characterize the hyperelastic behavior (C10, D). Secondly, an experimental compressiondevice is developed to measure the in vivo displacement field using ultrasound and Digital Image Correlation (DIC) techniques. Finally, an inverse method is implemented to identify the C10 and D parameters of each soft tissue. A numerical example is used to quantify the identification error on each parameter. Displacement field measurements confirm the ultrasound observations. They are also validated by the cartographies of the strain fields, which are obtained by the diffuse approximation method. Using the numerical example, the identification leads to low errors on the C10 (< 3%) and D (< 7%) parameters. Identified values of the mechanical parameters are in good agreement with the literature. All results validate the implemented identification method. In the long term, this protocol will allow to follow the evolution of pathologies and to conduct predictive simulations.

La trempe à eau par jet impactant est une étape clé des traitements thermiques subis par les tubes d'acier sans soudure. Elle permet, par le contrôle des évolutions microstructurales de l'acier, de conférer des propriétés mécaniques élevées, requises par l'utilisation des tubes dans des environnements de plus en plus extrêmes (nouvelles générations de centrales électriques, forages très profonds,…). En cours de trempe, des phénomènes complexes se produisent à la fois en surface des tubes et au sein du matériau. Les mécanismes au sein du tube (diffusion thermique – transformations de phases – contraintes - déformations) sont en effet étroitement liés aux phénomènes hydrauliques et thermiques survenant en surface de celui-ci. L'ensemble de ces phénomènes et leur couplage requièrent une parfaite compréhension et maîtrise pour anticiper l'apparition de certains défauts de trempe comme le cintrage des tubes ou les tapures et optimiser le traitement thermique. Nous avons donc mis en place au cours de ce travail de thèse, une démarche en deux temps. Dans un premier temps, à l'aide de deux dispositifs originaux de refroidissement, par un jet d'eau plan impactant, d'un cylindre tournant en Ni201 préchauffé à 600°C, nous avons étudié les effets de plusieurs paramètres opératoires sur les transferts thermiques: sous-refroidissement, vitesse du jet, vitesse de déplacement de la paroi, orientation et angle d'impact du jet. A l'issue de cette campagne, de nouvelles corrélations prédisant la densité de flux maximale échangée en paroi, qui prennent en compte la mobilité de la paroi, ont été proposées. En parallèle de ces expérimentations, nous avons entrepris la simulation numérique de la configuration expérimentale en utilisant le code de mécanique des fluides Fluent. Cependant la difficulté rencontrée pour simuler, à partir des modèles d'ébullition disponibles par défaut dans le code, l'ensemble des régimes de transfert thermique –et leur occurrence, n'a pas pu être surmontée. Dans un second temps, nous avons construit, au sein du centre de recherche de Vallourec, un dispositif expérimental de refroidissement hétérogène d'un tube en acier 42CrMo4 par un jet d'eau plan impactant la génératrice supérieure. La température du tube (en plusieurs points) et le déplacement vertical en son milieu étaient mesurés lors des essais. De plus, des visualisations à l'aide d'une caméra rapide ont permis de suivre le front de remouillage en tout début de trempe. En parallèle, nous avons réalisé des simulations de ces expériences pour calculer les évolutions de température, les cinétiques de transformations de phases et les évolutions des contraintes et des déformations au sein du tube d'acier. Pour cela nous avons utilisé le code de calcul par éléments finis, Sysweld, auquel est intégré le modèle de prédiction des cinétiques de transformations de phases, PhaseRC. Un jeu de paramètres thermophysiques, thermométallurgiques et thermomécaniques a été établi en se basant sur des travaux antérieurs et sur des caractérisations expérimentales menées au laboratoire.. L'analyse des résultats des simulations a permis de comprendre l'évolution de la flèche du tube en fonction des gradients de température et de la progression des transformations de phases au cours du refroidissement. La comparaison calcul - expérience a mis en évidence des écarts que nous avons tenté d'expliquer. Nous avons proposé des pistes pour poursuivre le travail entrepris au cours de cette thèse / Jet impingement quenching is one of the key steps among the heat treatments undergone by seamless steel tubes. By controlling the steel microstructural evolutions, the heat treatment leads to specific mechanical properties which are required by the use of the tubes in extreme environment (new generation of power plants, deep drilling,…) During the quenching, complex phenomena occur both at the pipe surface and within the material. The latter (thermal diffusion – phase transformations – stress – strain) are indeed closely related to the hydraulic and thermal mechanisms occuring on the surface. All these phenomena and their interconnection must be perfectly understood and controlled to prevent the appearance of quenching defects such as tube bending and quench cracks and to get the desired metallurgical phases. To achieve this goal we have adopted a two-stage approach. Firstly, using two innovative cooling experimental devices, including a rotating Ni201 cylinder preheated up to 600°C, we have studied the effect of several operating parameters on thermal transfer: subcooling, jet velocity, velocity of the displacement of the cooled surface, direction and impact angle of the jet. Following these sets of experiments, new correlations predicting the maximum of wall heat flux density have been proposed. These correlations take into account the wall motion. Besides this experimental work, we have undertaken numerical simulations of the experimental configuration using the CFD software Fluent. However, simulating all the thermal transfer regimes, particularly the boiling regime and their transitions is still challenging when using the default encoded boiling models; we were unable to successfully complete this work. In a second step, we have built an experimental device allowing heterogeneous quenching of a 42CrMo4 steel tube by an impinging water jet. During the quenching, the tube temperature at different locations and vertical displacement evolutions were recorded. In addition, data obtained from high-speed camera recordings allowed us to monitor the evolution of the rewetting front at the onset of cooling. Numerical simulations of these experiments were conducted in order to compute the time evolutions of temperature, phase transformations, stress and strain throughout the steel tube. To this end we used the finite element calculation software Sysweld, which includes a predictive model for the kinetics of phase transformations called PhaseRC. Using bibliographic data and some laboratory experimental characterizations, we were able to build a set of thermo-physical, metallurgical and thermo-mechanical data needful for these calculations. The simulation results have allowed to understand well the bending evolution of the tube during cooling, considering the thermal gradients and the progress of the phase transformation. Comparison of the simulation results with those obtained from measurements has shown discrepancies that we have tried to explain. Some suggestions have been given for the progress of the work accomplished during this thesis

Les matériaux composites présentent une forte rigidité pour une faible masse. Les méthodes courantes pour représenter le comportement vibratoire de ces matériaux ne sont souvent pas adaptées, car elles ne caractérisent pas le matériau mis en forme. Or, les propriétés dynamiques de celui-ci peuvent varier lors de sa fabrication et peuvent être dépendante de l’'espace. De nouvelles approches sont donc nécessaires pour mieux appréhender ces phénomènes.L’'approche proposée dans cette thèse utilise une méthode inverse locale, inspirée fortement de la méthode RIFF (Résolution Inverse Fenêtrée Filtrée) pour l’'identification des propriétés de matériaux. Ces travaux proposent d’'étendre le champ d’'application de cette méthode à des structures composites ayant une géométrie complexe, en remplaçant l’'opérateur analytique par un opérateur Éléments Finis. Les matériaux composites sont alors considérés comme homogènes et leurs propriétés sont recherchées. Dans le cas d’'une géométrie complexe deux couples de paramètres sont identifiés : le module d’'Young homogénéisé de traction complexe et le module d’'Young homogénéisé de flexion complexe, résultant du couplage existant entre les mouvements longitudinaux et transverses. Les méthodes inverses étant connues pour être sensibles aux incertitudes de mesure, une approche probabiliste est présentée pour régulariser le bruit de mesure. La régularisation est alors automatique et ne nécessite pas de paramètres à ajuster.L’'identification des paramètres structuraux, qui peut être globale ou locale, est présentée sur poutres droites, poutre courbes, plaques et coques. / Composite materials have high stiffness for a low mass. Common methods to represent the vibratory behavior of these materials are often not appropriate, as they do not characterize the material being shaped. However, the dynamic properties of the material can vary during its manufacture and can be dependent on space. New approaches are therefore needed to better understand these phenomena.The approach proposed in this thesis uses a local inverse method, strongly inspired by the FAT (Force Analysis Technique) for the identification of material properties. This work proposes to extend the scope of this method to composite structures with complex geometry, replacing the analytical operator with a Finite Element operator. Composite materials are then considered homogeneous and their properties are sought. In the case of a complex geometry two pairs of parameters are identified, the homogenized Young's modulus of complex traction and the homogenized Young's modulus of complex bending, resulting from the coupling between the longitudinal and transverse movements. As inverse methods are known to be sensitive to measurement uncertainties, a probabilistic approach is presented to regularize measurement noise. The regularization is then automatic and does not require any parameters to be adjusted.The identification of structural parameters, which can be global or local, is presented on straight beams, curved beams, plates and shells.

Résumé : La mesure des profils transitoires et de la vitesse de solidification sont deux données importantes pour le contrôle de procédés industriels impliquant un changement de phase. Dans le cas de l’électrolyse de l’aluminium, ce processus de solidification assure la protection du système et influe sur la performance énergétique du procédé de fabrication. Malheureusement, ces données se révèlent, dans la plupart des cas, difficilement accessibles. Ce travail de thèse porte sur le développement de nouveaux outils permettant l’étude et la caractérisation de la solidification de matériaux à changement de phase et à haute température. L’objectif est de développer un système de mesure du front de solidification de matériaux à changement de phase non destructif et ne perturbant pas le milieu de mesure, tout en assurant une précision et une réponse suffisamment rapide pour exploiter de nouvelles stratégies de contrôle dans les cuves d’électrolyse. Ce travail couple une étude expérimentale fondamentale de la solidification de la cryolithe avec une modélisation numérique de phénomène de changement de phase solide-liquide dans des conditions proches du fonctionnement de cuves d’électrolyse. // Abstract : Measurement of transient solidification fronts and of solidification rate are two important data for controlling industrial processes involving a solid-liquid phase change. In the case of aluminium electrolysis, this solidification process protects the system and affects the energy performance of the manufacturing process. Unfortunately, these data are not easy to obtain in most cases. This thesis focuses on the development of new tools for the study and on the solidification characterization of phase change materials at high temperature. The goal is to develop a nondestructive solidification front measurement system for phase change materials without disturbing the measurement medium, while ensuring accuracy and a fast enough response time to exploit new control strategies in electrolysis cells. This work couples a fundamental experimental study of the cryolite solidification with numerical modeling of solid-liquid phase change phenomenon under conditions close to those during normal operation of electrolytic cells.

La problématique de la sauvegarde de l’énergie pose un certain nombre de défis à la science, en particulier celui de son efficacité. La conception et la caractérisation de nouveaux matériaux isolants thermiques plus performants se révèlent donc fondamentales dans cette perspective. Les aérogels se présentent comme de sérieux candidats dans ce domaine, leur procédé de fabrication confère à certains d’entre eux des caractéristiques extrêmes telles qu’une grande porosité et une faible masse volumique. La caractérisation thermique de tels matériaux est délicate, leur faible sensibilité aux flux thermiques qui les traversent rend les méthodes connues difficiles à mettre en œuvre. A travers l’étude d’échantillons d’aérogels de faible poids moléculaire conçus au LCPM, une méthode de caractérisation adaptée a été développée. Cette méthode de type « tri-couche » offre les avantages d’être robuste et de s’affranchir de la connaissance de paramètres difficiles à atteindre dans de tels cas. La description et la validation de cette méthode sont l’objet principal de ce travail. Par ailleurs, les mesures de conductivité thermique sous vide ont été exploitées et ont permis une compréhension plus poussée de la structure de ces aérogels. Les résultats obtenus dans cette étude ouvrent donc des perspectives en vue de l’optimisation de nouvelles solutions pour l’isolation thermique / The issue of preserving energy raises a number of challenges to science, particularly its efficiency. The conception and characterization of new more efficient thermal insulating materials prove fundamental in this regard. Aerogels appear as serious candidates in this area, their manufacturing process provides extreme characteristics such as high porosity and low density for some of them. Thermal characterization of such materials is tricky, their low sensitivity to heat flux makes well-known methods difficult to implement. Through the study of low molecular weight aerogel samples designed by the LCPM a characterization method suitable to these samples has been developed by the LEMTA. This “three-layers” method offers the advantages of being robust and to overcome the knowledge of parameters that are difficult to reach in such cases. Describing and validating this method is the main object of this work. In addition, thermal conductivity measurements under vacuum have been processed which allowed a deeper understanding of the structure of aerogels. The results obtained this study open perspectives for the optimization of new solutions for thermal insulation

Les propriétés des argilites du Callovo-Oxfordien les désignent comme une barrière naturelle sûre pour le stockage de déchets radioactifs. Afin d’optimiser la prédiction de leur comportement, leur réponse à diverses sollicitations à court et long termes est étudiée. Notamment, lors de la phase d’exploitation, l’air ventilé dans les galeries n’est pas à l’équilibre hydrique avec la roche et peut perturber ses propriétés hydromécaniques (HM). Il semble alors essentiel de caractériser la réponse HM de la roche à des sollicitations hydriques.La démarche adoptée consiste à proposer un modèle physique adéquat, sur la base d’essais au laboratoire. L’essai de séchage est retenu puisqu’il permet d’étudier la réponse des échantillons en conditions non-saturées. A partir des résultats HM, un modèle élastique 2D isotrope est proposé, puis est élevé à la 3D avec un module de Young isotrope transverse. Ensuite, les données expérimentales servent à estimer certains paramètres poroélastiques et de transfert du modèle. Cette étape est accomplie par procédure inverse (minimisation de l’erreur mesures-calculs). Enfin, une modélisation 2D du comportement in situ est proposée, et compare les prédictions de modèles plastique isotrope et élastique isotrope transverse.Les simulations d’essais au laboratoire reproduisent assez bien les données expérimentales. Pour la modélisation in situ, une bonne corrélation est obtenue entre les prédictions et les mesures, et ce sans ajustement préalable des paramètres. Toutefois, le modèle souligne une influence limitée de la plasticité à l’échelle du laboratoire, alors que les phénomènes dissipatifs sont marqués in situ. Les modélisations 3D au laboratoire ne donnent pas de résultats plus fins qu’en 2D, mais reproduisent plus de données expérimentales (variations de masse, déformations axiales et latérales). De plus, injecter plusieurs types de données dans la formulation inverse permet d’améliorer la précision de l’algorithme. Par ailleurs, une meilleure stabilité de l’algorithme est obtenue en adoptant une convergence en deux étapes (simplex, puis méthode de type gradient). Les estimations numériques des paramètres corroborent les mesures expérimentales obtenues par ailleurs. / The Callovo-Oxfordian claystones’ properties make them reliable as a geological barrier for the confinement of radioactive wastes. In order to optimally predict their behavior, how they respond to various short and long terms loadings has to be studied. Particularly during the exploitation phase, air is continuously ventilated throughout the galleries. The climatic properties of this air are not balanced with those of the rock, and may perturb its hydromechanical (HM) attributes. Thus, assessing the HM response of the rock under hydric loading seems to be a priority.This dissertation begins with laboratory tests to propose an appropriate physical model. Drying tests were studied as they focus on the HM response of samples undergoing hydric loadings. A first 2D isotropic model is proposed, and then enhanced to 3D by considering a transversely isotropic Young modulus. Secondly, experimental results provide relevant data to estimate poroelastic and transport parameters involved in the model. Estimation is achieved according to an inverse procedure, which minimizes the error between measurements and model predictions. Finally, a real-size test is simulated using 2D models: an isotropic plastic one and a transversely isotropic elastic one.Model predictions reproduce well the laboratory tests data. When simulating the in situ behavior, a rather good agreement is obtained between the numerical and experimental results (although using the parameters estimated at the laboratory scale). However, the model highlights a limited influence of plasticity in the laboratory tests, while dissipative phenomena obviously occur in situ. 3D laboratory simulations do not improve the precision of 2D results, but reproduce more experimental data (mass variations, axial and lateral strains). Moreover, the inversion process is more efficient when ran over various kinds of data. Furthermore, stability of the algorithm is improved when adopting a two-phase convergence (simplex, followed by a gradient-like method). Numerical estimates of the parameters are in agreement with the direct experimental measurements obtained through other tests.

EDF exploite sur son parc de production de nombreuses machines électriques tournantes. Les contraintes thermiques subies par celles-ci engendrent des échauffements locaux qui nuisent à leur intégrité. Le présent travail contribue à fournir des méthodes de calcul adaptées à la détection et à la localisation des points chauds. Il participe à améliorer la compréhension des écoulements en rotation et leurs effets sur les transferts thermiques. Plusieurs dispositifs expérimentaux, de complexité ascendante, ont été utilisés pour comprendre et valider les simulations numériques. Une première étude sur une marche descendante (demi-pôle) parcourue par un jet de paroi non-confiné a mis en avant des différences par rapport à un jet confiné ; ces deux cas existent dans un alternateur. Une seconde étude menée sur une cavité tournante confinée a analysé l’impact d’un écoulement de Taylor-Couette-Poiseuille sur la température et la position des points chauds créés, en balayant l’ensemble des régimes d’écoulement. Ces études ont mis en exergue une première méthode de calcul fiable, fondée sur l’étude numérique CHT. Une autre méthode, basée sur la FEM couplée à une méthode inverse, a été testée sur une maquette d’alternateur hydraulique afin de pallier aux temps de calcul longs de la première. Cette méthodologie remonte aux coefficients d’échanges convectifs numériques à partir des mesures du champ thermique du rotor, mais n’est envisageable que lorsque l’on dispose de données expérimentales suffisantes. Ces travaux ont aussi mis en évidence de nouvelles techniques de mesures sans contact, comme l’utilisation d’un pyromètre à haute fréquence pour la mesure de température sur des machines tournantes. / EDF operates a large number of electrical rotating machines in its electricity generation capacity. Thermal stresses which affect them can cause local heating, sufficient to damage their integrity. The present work contributes to provide methodologies for detecting hot spots in these machines, better understanding the topology of rotating flows and identifying their effects on heat transfer. Several experimental scale model were used by increasing their complexity to understand and validate the numerical simulations. A first study on a turbulent wall jet over a non-confined backward-facing step (half-pole hydrogenerator) notes significant differences compared to results from confined case : both of them are present in an hydrogenerator. A second study was done on a small confined rotating scale model to determinate the effects of a Taylor-Couette-Poiseuille on temperature distribution and position of hot spots on the heated rotor, by studying the overall flow regimes flow. These studies have helped to obtain a reliable method based on conjugate heat transfer (CHT) simulations. Another method, based on FEM coupled with the use of an inverse method, has been studied on a large model of hydraulic generator so as to solve the computation time issue of the first methodology. It numerically calculates the convective heat transfer from temperature measurements, but depends on the availability of experimental data. This work has also developped new no-contact measurement techniques as the use of a high-frequency pyrometer which can be applied on rotating machines for monitoring temperature.