Les éboulements rocheux sont difficilement prévisibles en raison de leur soudaineté et du peu d'informations disponibles sur la structure interne du massif rocheux. L'analyse en retour d'éboulements dans des falaises calcaires a montré qu'ils étaient initiés par la fracturation de ponts rocheux et que leur persistance était le principal paramètre à prendre en compte pour des calculs de stabilité. En 1er lieu, nous avons utilisé l'information contenue dans le bruit sismique afin d'évaluer le degré d'instabilité d'une écaille rocheuse. Les données viennent de la falaise calcaire de Chamousset (Vercors) où un réseau sismique a fonctionné 4 mois jusqu'à 2 semaines avant la chute d'une écaille de 21000 m3. Le calcul des spectres de bruit sismique a permis d'identifier la 1ère fréquence propre de l'écaille qui diminue significativement avant l'éboulement. Ce même réseau a enregistré des milliers d'événements dont les plus énergétiques ont été classés. Certains de ces signaux ont des arrivées d'ondes P et S distinctes, permettant d'analyser leur mécanisme de rupture et de calculer leur magnitude et distance hypocentrale. Des simulations favorisent l'hypothèse de ruptures en traction pour expliquer les polarisations. D'autres simulations confirment que les événements situés à l'interface écaille-massif entraînent la résonance de l'écaille, expliquant l'observation de codas mono-fréquentielle avec une amplitude sismique plus forte sur l'écaille. Enfin, l'étude de la réponse thermomécanique d'une 2nde écaille sur un an montre l'influence des paramètres météorologiques sur ses fréquences propres. Les effets réversibles et irréversibles sur ces fréquences ont été identifiés.

L'augmentation des performances des moteurs conduisent à un accroissement des températures dans les collecteurs d'échappement. Les fontes ferritiques, GS50 ou SiMo, utilisées par RENAULT présentent des évolutions microstructurales à température élevée (800oC) qui peuvent compromettre leur résistance. Plutôt que d'adopter des matériaux plus résistants, mais plus chers, une meilleure connaissance des évolutions microstructurales des fontes GS et de leur comportement permettraient d'optimiser la géométrie des pièces dans la recherche d'un compromis légèreté-fiabilité. Une expertise de collecteurs vieillis sur banc, comparés aux mêmes pièces à l'état brut, a permis de recenser trois principaux modes d'endommagement : oxydation, décarburation et transformation de phases. Ils ont été étudiées par des essais de vieillissement isotherme et anisotherme. Les cinétiques d'oxydation et de décarburation ont été déterminées à différentes températures. La connaissance des mécanismes de vieillissement ont permis de reproduire sur des éprouvettes massives certaines des microstructures observées. Les matériaux décarburés et bruts de coulée ont été comparés lors d'essais de dilatométrie et de fatigue oligocyclique. Et des lois de comportement élasto-viscoplastiques des nuances brutes et vieillies ont été identifiées. Des essais de fatigue thermique ont été simulés avec une bonne précision à partir des lois de comportements déterminées par des essais isothermes. Les résultats de simulation d'essais thermomécaniques sur des matériaux en partie décarburés, effectués sur une géométrie simple, sont suffisamment satisfaisants pour envisager d'étendre ces calculs à des collecteurs d'échappement.

L'objectif de cette étude est de fournir un outil d'aide à la conception de gammes d'emboutissage à chaud. Ce procédé appliqué à des aciers trempants (dans notre cas à l'Usibor 1500) permet l'obtention de pièces avec une résistance mécanique très élevée. La simulation numérique représente un outil de prédiction très puissant permettant de minimiser les coûts de conception. Grâce à la simulation numérique et aux analyses nécessaires à la modélisation, nous sommes arrivés à élucider les points clés qui font qu'une pièce emboutie par ce procédé soit réussie. Deux logiciels, Forge 3 et Abaqus/Explicit, ont été choisis. Ces logiciels permettent d'effectuer un calcul thermomécanique indispensable à la simulation de notre procédé.<br/> Pour réaliser ces calculs, la connaissance du comportement du matériau ainsi que celles des échanges thermiques gérant l'évolution de température dans la tôle sont cruciales. Le comportement mécanique comme l'évolution de la température sont décrits par des lois mathématiques dépendant d'un certain nombre de paramètres propres à chaque matériau. Ces paramètres ont été identifiés pour notre matériau dans les conditions du procédé d'emboutissage à chaud.

Nous avons travaillé au cours de cette étude sur des plaques stratifiées en verre/époxy fabriquées à partir de préimprégnés pseudo-tissés. Ces composites sont sujets au délaminage et à la fissuration intralaminaire. L'étude se décompose en quatre grandes étapes. Tout d'abord, une campagne d'essais a été menée, au cours de laquelle les caractéristiques des endommagements intralaminaires transverses et leur cinétique d'évolution ont été déterminées. Cet endommagement présente beaucoup de points communs avec celui observé dans les matériaux non tissés. Toutefois, la structure tissée du matériau a des conséquences sur certaines de ces caractéristiques, en particulier sur sa cinétique d'évolution. Une deuxième campagne d'essais a été réalisée afin d'étudier le phénomène d'initiation du délaminage aux voisinages des bords. Dans ce but, une méthode d'observation optique innovante a été établie afin de pouvoir visualiser durant les essais les différents mécanismes de dégradation présents. Nous avons ainsi étudié l'influence de la fissuration intralaminaire et de la structure tissée sur le début du délaminage. En parallèle, ces endommagements ont été modélisés. Ainsi, un modèle de comportement à l'échelle mésoscopique, prenant en compte les dégradations intralaminaires, a été mis en place. Les effets dus à la refermeture incomplète des fissures a également été prise en compte dans le modèle. De plus, un critère d'initiation du délaminage a été établi. Le critère est construit à partir des composantes du vecteur contrainte de l'interface, moyennés sur une distance caractéristique. Enfin, une dernière phase a permis de valider le modèle et le critère sur des structures non trouées et trouées.

L'objet de ces travaux est de se doter d'outils de modélisation qui constituent aujourd'hui un passage incontournable pour maîtriser le procédé R.T.M. Ces outils reposent sur la prédiction de grandeurs macroscopiques caractérisant les renforts. La représentation fidèle du renfort à l'échelle microscopique dans une cellule élémentaire est la clé de voûte de notre travail. Une méthode d'initialisation de fonction caractéristique a ainsi été mise au point permettant de déterminer fidèlement les interfaces entre la résine et les renforts. Cette méthode possède la particularité d'être applicable à des renforts à structure périodique mais aussi aléatoire. La simulation de l'écoulement 3D de la résine dans les pores entre les fibres permet de déterminer la perméabilité de la structure représentée. Le passage de l'échelle microscopique à l'échelle macroscopique est réalisé en homogénéisant les champs de vitesse et de pression obtenus numériquement dans la cellule élémentaire.Une démarche similaire a été développée pour étudier les phénomènes thermiques dans les milieux poreux. Cette approche permet de déterminer des valeurs de conductivité thermique équivalente en considérant dans un premier temps, le fluide au repos et ensuite en couplant les phénomènes thermiques à l'écoulement engendrant des fluctuations thermiques à l'échelle microscopique. Ce phénomène est important et modifie grandement les résultats. Sa prise en compte introduit des termes non linéaires couplant la conductivité thermique équivalente à la vitesse d'écoulement du fluide au travers du renfort.La connaissance précise de ces paramètres permet la réalisation de simulations 3D réalistes du procédé R.T.M.

L'instabilité de Vishniac est invoquée pour expliquer la fragmentation et la filamentation de la coquille fine de matière choquée présente dans les restes de supernova radiatifs. Toutefois l'implication et les conséquences de ce processus spécifique sur la morphologie complexe des restes de supernova ne sont pas entièrement démontrées. Nous avons réalisé des simulations numériques 2D d'une onde de souffle perturbée qui se propage dans un milieu ambiant homogène afin de mieux comprendre le mécanisme de l'instabilité de Vishniac. Le code hydrodynamique HYDRO-MUSCL 2D a été utilisé afin de déclencher la perturbation de la coquille fine dans les géométries plan-parallèle et sphérique. Dans ce travail, je montrerai que nous avons obtenu la surstabilité de Vishniac (" Vishniac overstability ") comme prédit par l'analyse théorique. Ce processus repose sur une oscillation d'amplitude et de période temporelle croissante des grandeurs fluides et spatiales. Nous avons testé l'effet de paramètres spécifiques sur la variation de masse d'une région donnée dans le cadre d'une étude paramétrique que nous avons réalisée sur le supercalculateur Titane. Nous avons trouvé que la perturbation est atténuée après quelques oscillations de la variation de masse pour tous les jeux de paramètres. Ainsi nous concluons que dans notre modèle, l'instabilité de Vishniac ne permet pas la fragmentation de la coquille fine à cause d'effets qui ne sont pas pris en compte par l'analyse théorique.

Les travaux présentes s'inscrivent dans le cadre de la réduction de variance pour les méthodes de Monte-Carlo et plus généralement dans l'optimisation de méthodes numériques à l'aide de couplage entre des méthodes déterministes et des méthodes probabilistes. Trois thèmes principaux seront abordés à l'aide de ces techniques: l'intégration numérique sur un hypercube, la résolution d' équations aux dérivées partielles linéaires et le calcul des éléments propres principaux (valeur propre et vecteur propre) de certains opérateurs linéaires.

Au cours des dernières décennies, les retours d'expérience de tunnels ont mis en évidence de nombreux dommages causés par le phénomène de gonflement. Ce phénomène, qui concerne des formations variées très répandues dans les ères géologiques (terrains argileux ou anhydrite), est susceptible d'affecter un certain nombre de tunnels futurs. Cette thèse envisage de mettre au point une méthodologie de calcul des revêtements de tunnels en terrain gonflant.<br /><br />En s'inspirant des retours d'expérience et des méthodes existantes, une loi de comportement élastoplastique a été établie pour simuler le développement de gonflement apparaissant dans le massif du fait d'un changement d'état de contraintes et d'une arrivée d'eau. La loi de gonflement supposée élastique et tridimensionnelle relie le logarithme de la contrainte à la déformation, en tenant compte d'une éventuelle anisotropie de gonflement associée à la stratification du terrain.<br /><br />Après avoir été implanté dans le code par éléments finis CESAR-LCPC, ce modèle de gonflement a été validé à deux niveaux : d'une part, sur un échantillon de matériau gonflant soumis à un essai de gonflement et, d'autre part, sur un ouvrage instrumenté directement concerné par le phénomène de gonflement.<br />Des essais de gonflement spécifiques permettent d'identifier les différents paramètres de gonflement. Les tentatives de calage de la loi de gonflement ont donné des résultats satisfaisants sur plusieurs matériaux gonflants différents.<br />La modélisation numérique du tube Nord du tunnel de Chamoise pour lequel les mesures in situ ont mis en évidence le développement du gonflement des marnes sous le radier, a montré une meilleure concordance du modèle de gonflement avec les mesures in situ que les calculs classiques ; en particulier, on a pu retrouver une asymétrie marquée du comportement du radier et obtenir des déplacements dans le massif plus proches de la réalité.

Dans cette thèse, nous modélisons le passage d'une onde acoustique à travers un matériau poreux supposé périodique (mousse ou laine de verre). L'objectif est d'établir un système d'équations macroscopiques, prenant en compte la microstructure du domaine, sur un matériau homogène équivalent. Nous commençons par définir la modélisation du problème et nous obtenons un système non coercif écrit en fréquence afin de répondre aux problématiques industrielles. Nous établissons ensuite le caractère bien posé du système avant de le faire converger double échelle. Nous obtenons ainsi deux systèmes de cellule : un système de Stokes pour le fluide et un système de type élasticité linéaire pour décrire le déplacement du matériau. Au niveau macroscopique, nous obtenons un système couplé non intuitif. Nous finissons la partie théorique par une comparaison entre le modèle que nous avons obtenu et le modèle de Biot-Allard utilisé par les physiciens et industriels pour traiter ce type de problème. Enfin, nous illustrons le travail précédent grâce à des résultats numériques.

Cette étude est hébergée par un programme hospitalier de recherche clinique. L'objectif était de développer un outil d'évaluation de la motricité de la main pour caractériser la réorganisation motrice associée à la restauration d'une paralysie radiale (technique de Tsugé). Nous avons montré que la minimisation des moments secondaires était un principe biomécanique robuste pour expliquer les interactions biomécaniques entre les doigts. Ce principe reste toutefois à mieux explorer en extension. L'utilisation du modèle biomécanique a montré la nécessité d'inclure l'équilibre du poignet dans la procédure. L'analyse de la réorganisation motrice montre la redistribution dynamique des tensions des tendons et la mise en place de co-contraction. L'ensemble de cette démarche ouvre des perspectives d'analyse de la motricité de la main intéressantes.