L'objectif de cette thèse est de comprendre et modéliser, sur des bases physiques, les phénomènes métallurgiques et mécaniques induisant des changements de textures cristallographiques. Des caractérisations microscopiques et texturales sont menées sur un acier sans interstitiel (IF pour Interstitial Free) au cours de son procédé d'élaboration : à l'état de tôle à chaud, puis de tôle à froid et enfin après recristallisation. Parallèlement à cette analyse expérimentale, un couplage de modèles de prévision des textures de déformation et de recristallisation est mené. Cela concerne principalement la prévision des paramètres de la déformation déterminants pour la genèse des textures de recristallisation. La caractérisation des hétérogénéités de déformation tout comme la modélisation a été réalisée à différentes échelles. Les paramètres extraits de la caractérisation expérimentale, notamment grâce aux analyses EBSD (au MEB et au MET), ont permis d'estimer la fragmentation qui a lieu pour certaines orientations cristallographiques. Ces résultats s'accordent avec l'estimation de l'énergie stockée par DRX avec, en moyenne deux fois plus d'énergie dans les grains de la fibre γ après de fortes réductions mais peu de différences avant. Une hiérarchie de la fragmentation F est proposée telle que F{100}<110> < F{112}<110> < F{111}<110> < F{111}<112> < F{554}<225> . Des comparaisons expérience-modèles ont complété ces analyses. La fragmentation de la fibre γ s'explique par un nombre de murs créés au cours de la déformation plus important que pour α. D'autre part les axes de désorientations de ces deux fibres sont différents, une rotation autour d'un axe DL pour la fibre α, et DN pour la fibre γ.

Les pièces multimatériaux à gradient fonctionnel (Functionally Graded Materials - FGM) sont des structures dont la composition et la microstructure du matériau changent graduellement à l'intérieur de la pièce. Cette distribution des matériaux permet de réaliser des gradients de propriétés au niveau mécanique, physique, chimique, etc. Les domaines d'application sont nombreux pour ces pièces, en particulier l'aérospatial et le biomédical mais également l'électronique, l'énergie nucléaire, la production d'outillage, le design, etc. L'utilisation des procédés innovants tels que les procédés de fabrication additive est indispensable pour la réalisation de pièces multimatériaux complexes. Bien que ces procédés aient les caractéristiques attendues pour la réalisation de pièces multimatériaux, on constate qu'aucune pièce fonctionnelle n'a encore été fabriquée à ce jour. Pour permettre la fabrication de pièces fonctionnelles, il est indispensable de proposer une méthodologie de fabrication complète permettant de passer de l'objet imaginé par le concepteur à la fabrication. Cette méthodologie doit comporter les étapes suivantes : description de la pièce à fabriquer, détermination d'une stratégie de fabrication adaptée et génération des instructions de fabrication. Parmi les étapes du processus de fabrication, celle de choix d'une stratégie de fabrication occupe une place importante. En effet, les caractéristiques de pièces - géométrie et répartition des matériaux - sont fortement dépendantes de la stratégie de fabrication choisie. Les travaux de thèse portent principalement sur les méthodes mises en place pour la détermination de trajectoires appropriées à la fabrication des pièces multimatériaux. Ces méthodes reposent sur la modélisation du procédé nécessaire à l'évaluation des stratégies et une optimisation du procédé permettant de diminuer les différences entre la répartition des matériaux souhaitée et celle fabriquée. Une des méthodes proposées permet d'obtenir automatiquement des trajectoires parfaitement adaptées aux pièces multimatériaux et repose sur la modélisation et l'optimisation du procédé. Ces travaux sont intégrés dans une méthodologie de fabrication de pièces multimatériaux. De plus, une maquette informatique a été développée pour mettre en avant les possibilités d'utilisation de cette méthodologie.

Ce travail s'inscrit dans une démarche visant à quantifier l'impact de la diagenèse sur les ratios élémentaires Sr/Ca, Mg/Ca et U/Ca caractéristiques de l'exosquelette (aragonite biogénique) de coraux Porites sp. Ces ratios élémentaires sont très utilisés comme paléo thermomètres en paléo climatologie tropicale pour reconstruire les paléo variables environnementales comme la température de surface de l'océan. On considère dans un premier temps un corail moderne prélevé in vivo et un corail fossile daté au 14C de l'Holocène Moyen (5445 ans BP). Les deux coraux ne présentent aucune trace d'altération diagénétique. On mesure Sr/Ca, Mg/Ca, U/Ca dans l'exosquelette par spectrométrie de masse couplée à un plasma inductif et à la microsonde de Castaing pour Sr/Ca. Les températures de l'océan du temps du vivant des coraux sont alors reconstruites à partir de la géochimie des échantillons et validées. L'approche est ensuite étendue à des coraux pléistocènes (~125000 ans BP) altérés par la diagénèse e.g., apparition de calcite au détriment de l'aragonite initiale sous l'action de l'eau douce percolant dans les récifs. Cette calcitisation affecte les ratios Sr/Ca et Mg/Ca originels et est susceptible d'entacher les reconstructions d'un biais qu'il est important de bien connaître. En s'appuyant sur la spectrométrie Raman et la microanalyse X, on confirme les tendances qualitatives rapportées dans la littérature sur l'impact de la calcite sur les proxies du climat : la présence de calcite diminue (resp. rehausse) le ratio Sr/Ca (resp. Mg/Ca) créant donc des artéfacts chauds dans les températures reconstruites.

L'éolienne Darrieus connaît un intérêt accru ces dernières années parce qu'elle représente une solution alternative potentielle de production d'électricité dans les milieux urbains. En particulier, une éolienne de forte solidité peut être choisie car certaines de ses propriétés peuvent être avantageuses pour son implantation proche de zones habitées. A l'inverse, certaines difficultés aérodynamiques émergent. Ce type d'éolienne fonctionne à de faibles vitesses réduites pour lesquelles le décrochage dynamique a un rôle très significatif. L'objectif de ce travail de thèse consiste à compléter la connaissance du phénomène de décrochage dynamique sur une éolienne à axe vertical afin d'améliorer les modèles numériques de prédiction existants. Cette étude s'appuie sur une analyse combinée de résultats numériques et expérimentaux. Les simulations numériques sont produites avec une méthode des panneaux bidimensionnelle instationnaire. Les effets de la viscosité sont introduits par des corrections utilisant notamment un modèle semi-empirique de décrochage dynamique. Le travail expérimental s'est concentré sur la dynamique tourbillonnaire à proximité immédiate du rotor résultante du décrochage dynamique. Le montage se compose d'une éolienne à pale droite placée dans une soufflerie. Des mesures instationnaires de la répartition de pression pariétale le long de la corde et des mesures de champ de vitesse par vélocimétrie par images de particules ont été accomplies. Les résultats révèlent la manière dont les caractéristiques du décrochage dynamique sont conditionnées par la vitesse réduite. Le retard au décrochage, l'intensité de l'effet du tourbillon de décrochage dynamique et sa convection ont été quantifiés. Enfin, un examen critique de l'applicabilité du modèle de Leishman-Beddoes pour simuler efficacement les effets du décrochage dynamique a été réalisé.

La séparation aveugle de sources est une technique prometteuse pour l'identification, la localisation, et la classification des sources sonores. L'objectif de cette thèse est de proposer des méthodes pour séparer des sources sonores incohérentes qui peuvent se chevaucher à la fois dans les domaines spatial et fréquentiel par l'exploitation de l'information spatiale. De telles méthodes sont d'intérêt dans les applications acoustiques nécessitant l'identification et la classification des sources sonores ayant des origines physiques différentes. Le principe fondamental de toutes les méthodes proposées se décrit en deux étapes, la première étant relative à la reconstruction du champ source (comme par exemple à l'aide de l'holographie acoustique de champ proche) et la seconde à la séparation aveugle de sources. Spécifiquement, l'ensemble complexe des sources est d'abord décomposé en une combinaison linéaire de fonctions de base spatiales dont les coefficients sont définis en rétropropageant les pressions mesurées par un réseau de microphones sur le domaine source. Cela conduit à une formulation similaire, mais pas identique, à la séparation aveugle de sources. Dans la seconde étape, ces coefficients sont séparés en variables latentes décorrélées, affectées à des "sources virtuelles" incohérentes. Il est montré que ces dernières sont définies par une rotation arbitraire. Un ensemble unique de sources sonores est finalement résolu par la recherche de la rotation (par gradient conjugué dans la variété Stiefel des matrices unitaires) qui minimise certains critères spatiaux, tels que la variance spatiale, l'entropie spatiale, ou l'orthogonalité spatiale. Il en résulte la proposition de trois critères de séparation à savoir la "moindre variance spatiale", la "moindre entropie spatiale", et la "décorrélation spatiale", respectivement. De plus, la condition sous laquelle la décorrélation classique (analyse en composantes principales) peut résoudre le problème est établit de une manière rigoureuse. Le même concept d'entropie spatiale, qui est au coeur de cette thèse, est également iv exploité dans la définition d'un nouveau critère, la courbe en L entropique, qui permet de déterminer le nombre de sources sonores actives sur le domaine source d'intérêt. L'idée consiste à considérer le nombre de sources qui réalise le meilleur compromis entre une faible entropie spatiale (comme prévu à partir de sources compactes) et une faible entropie statistique (comme prévu à partir d'une faible erreur résiduelle). La méthode proposée est validée à la fois sur des expériences de laboratoire et des données numériques et illustrée par un exemple industriel concernant la classification des sources sonores sur la face supérieure d'un moteur Diesel. La méthodologie peut également séparer, de façon très précise, des sources dont les amplitudes sont de 40 dB inférieur aux sources les plus fortes. Aussi, la robustesse vis-à-vis de l'estimation du nombre de sources actives, de la distance entre le domaine source d'intérêt et le réseau de microphones, ainsi que de la taille de la fonction d'ouverture est démontrée avec succès.

Le développement de techniques de mesure optiques de surfaces libres permettant la mesure des ondes générées par un modèle de navire en bassin des carènes est proposé dans ce mémoire. Trois méthodes, toutes basées sur un système de prises de vues stéréoscopiques, sont présentées. La première repose sur l'emploi de procédures d'inter-corrélation des images des caméras composant le système de stéréovision. La seconde s'appuie sur la théorie de la réfraction de la lumière au travers de l'interface air/eau. La troisième méthode est basée sur l'identification des projections de faisceaux laser sur la surface libre vue par un système stéréoscopique. Ces méthodes sont appliquées à la mesure du champ de vague autour d'un modèle de navire tracté dans le bassin des carène de l'Institut Pprime et dans un bassin des carènes de grande taille (B600 de la DGA Techniques hydrodynamique) et permettent la reconstruction et l'analyse du sillage. Les principales caractéristiques des champs de vagues et leurs dépendances au nombre de Froude sont étudiées en détail. Afin de compléter l'étude de l'écoulement généré par le modèle de navire, des mesures de PIV stéréoscopique autour de la carène, des mesures spécifiques de la vague d'étrave et des simulations numériques sont réalisées. L'ensemble de ces données permettent l'analyse de champs de vitesse tridimensionnelle proche de la coque et l'identification des caractéristiques du champ de vague de l'étrave jusqu'au champ lointain.

Ce travail est dédié au traitement numérique de la fissuration au sein d'une structure dont la rupture résulte d'un endommagement ductile. Trois points importants sont ainsi à considérer : la modélisation du comportement du matériau, la modélisation de la fissure au sein de la structure, la méthode de propagation couplant les deux modélisations précédentes. L'étude est conduite en 2D en déformations planes (DP), dans le cadre des petites perturbations. Une procédure dynamique a été choisie avec un schéma de résolution implicite, en utilisant le code de calculs industriel Abaqus. Le comportement du matériau est modélisé via un modèle de type GTN modifié (GTNm), voir Longère et al, Int. J. Dam. Mech. (2012), prenant en compte les effets combinés de l'écrouissage, de l'adoucissement thermique, de la viscoplasticité et de l'endommagement ductile. Le modèle GTNm a été implanté en 2D DP dans Abaqus via une routine utilisateur (umat). Pour décrire les conséquences cinématiques de la présence d'une fissure au sein de la structure ainsi que de sa propagation, nous avons retenu la méthode des éléments finis étendus (X-FEM), voir Moes et al, Int. J. Numer Meth. Engng (1999). La méthode X-FEM a été implantée en 2D DP dans Abaqus via une routine utilisateur (uel). Afin de permettre le couplage entre la X-FEM et le comportement fortement non linéaire du matériau, certains choix ont dû être faits, notamment en ce qui concerne la forme du champ de déplacement enrichi ainsi que l'intégration numérique. Comme dans le travail de Haboussa et co-auteurs, voir Haboussa et al, Int. J. Numer Meth. Engng (2011), certaines grandeurs sont moyennées dans une pastille située en pointe de fissure dans le but d'atténuer la dépendance des résultats numériques au maillage. L'originalité de ce travail réside sans doute dans la méthode de propagation, dont le but est de reproduire une rupture induite par endommagement ductile. En premier lieu, la direction de propagation potentielle de la fissure est calculée via une analyse de bifurcation dans le cas d'un état matériau figé. Le principe est en effet de considérer que la direction de propagation de la fissure est la direction dans laquelle la déformation s'est localisée, comme c'est le cas dans la phase de coalescence Dans un second temps, un critère d'amorçage basé sur une énergie stockée en pointe de fissure a été défini. Si cette énergie dépasse une valeur critique, nous considérons que la fissure se propage. Nous calculons alors la longueur du nouvel incrément de fissure via une des deux méthodes retenues dans notre étude : une méthode dite par épuisement (MEPUI) et une méthode par contrôle de la vitesse de fissuration (MCVF). La MEPUI est basée sur le principe que la fissure se propage tant que le critère d'amorçage de la fissure est vérifié et est valable quelle que soit la vitesse de chargement de la structure. A contrario, la MCVF est basée sur la loi de Kanninen qui permet de contrôler la vitesse de fissuration et n'est valable que pour des vitesses de chargement élevées. Ces deux méthodes de propagation (MEPUI et MCVF) ont été implantées dans Abaqus via une routine utilisateur (uel). Afin de tester les deux méthodes mises en place dans notre étude (MEPUI et MCVF), des simulations numériques sur des plaques fissurées et entaillées ont été conduites. Considérant quelques simplifications, ce travail permet de reproduire la propagation d'une fissure dans une structure 2D dont la rupture résulte d'un endommagement par germination et croissance de cavités.

Developing numerical tools for polymer lifetime prediction constitutes a promising opportunity for shortening the duration of material certification procedures in the automotive industry without decreasing their reliability. This PhD thesis aims at modeling the photothermal oxidation of isotactic polypropylene (iPP), which is responsible for the alteration of both its mechanical and aspect properties. The adopted approach consists in coupling the kinetics of photo- and thermo-oxidation reactions with physical phenomena, such as oxygen transport and UV-light attenuation in the material thickness, in order to describe all the physico-chemical changes. Upper-scale properties, from which will be defined the end-of-life criteria, will be calculated afterwards by applying the suitable structure-property relationships. The main challenge was to extend the pre-existing kinetic model of thermal ageing to photothermal ageing by taking into account initiation reactions of photolysis. Heavy campaigns of ageing and characterization tests made on a reference iPP, as well as an exhaustive capitalization of literature data of other iPPs, have allowed elaborating a kinetic model of photothermal oxidation and to generalize it to the whole iPP family in large domains of oxygen partial pressure (from 0.2 to 50 bars), temperature (from 40 to 230°C) and UV-light exposure (variable intensities and light sources) describing both natural and accelerated ageing conditions. The experimental validation of the model has allowed substantiating the kinetic approach and showing its limitations, as well as highlighting some numerical issues. The model has been designed in order to be an upgradable numerical tool which will allow, at term, optimizing the representativeness of the ageing testing devices and the performance of commercial iPP formulations. All these theoretical and numerical developments are prone to be applied to the photothermal degradation of other types of polymer substrates, but also in other application fields of the macromolecular photochemistry such as UV-photopolymerization.Keywords: Polypropylene, photothermal oxidation, oxygen diffusion control, screen effect, kinetic modeling, lifetime prediction.

La dispersion dans des milieux poreux homogènes (empilements de grains) a été étudiée par des mesures par résonance magnétique nucléaire (RMN) et des simulations de marches aléatoires dans un réseau de pores. La RMN permet de mesurer l'ensemble des déplacements des molécules d'eau durant un temps tΔ, et d'obtenir propagateurs et moments caractéristiques. L'évolution temporelle du second moment σ (σ2 ∝ taΔ) permet de caractériser de manière précise le régime de dispersion des molécules (Gaussien ou anormal). Des mesures pour des écoulements de 15 < Pe < 45 dans un empilement de grains de 30μm ont permis d'observer une dispersion anormale faiblement super-dispersive (a = 1.17) en écoulement saturé et une augmentation progressive du caractère super-dispersif avec la diminution de la saturation en eau (jusqu'à a = 1.5 pour 42 %)lors d'une co-injection stationnaire eau-huile. En écoulement saturé, les propagateurs et courbes de percée sont quasi-gaussiennes, tandis qu'en écoulement insaturé, les propagateurs sont asymétriques et les courbes de percée présentent des trainées aux grands temps. Dans ces conditions, on montre que la dispersion anormale observée est mieux décrite par des lois stables de Lévy que par des lois gaussiennes. Des simulations de marche aléatoire ont été réalisées dans un réseau de pores extrait d'un milieu poreux réel par imagerie microscanner.Elles permettent d'obtenir les mêmes informations que la RMN, les marcheurs se déplaçant par advection et diffusion. Ces simulations montrent l'existence d'une stagnation non observée dans les expériences, montrant que la simplification du réseau poreux est trop importante et empêche de reproduire certains aspects du champ de vitesses détecté par la RMN. Toutefois, l'évolution temporelle du second moment a également un caractère super-dispersif à temps long à 100 % de saturation

Les matériaux composites stratifiés sont de plus en plus utilisés dans l'aéronautique, mais ils peuvent être sujets à large délaminage si soumis à impact. La nécessité d'effectuer des simulations numériques pour prédire l'endommagement devient essentielle pour l'ingénieur. Dans ce contexte, l'utilisation d'une modélisation fine semble préférable. En revanche, le coût de calcul associé serait prohibitif pour larges structures. Le but de ce travail consiste à réduire ce coût de calcul, en couplant le modèle fin, restreint à la zone active de délaminage, avec un modèle grossier appliqué au reste de la structure. En raison du comportement transitoire des problèmes d'impact, l'adaptabilité dynamique des modèles pour suivre les phénomènes évolutifs représente un point crucial de la stratégie de couplage. Des méthodes avancées sont utilisées pour coupler différents modèles. Par exemple, la méthode de Décomposition de Domaines, appliquée à l'adaptabilité dynamique, doit être combinée avec une stratégie de remaillage, considérée comme intrusive pour la mise en œuvre d'un logiciel pour Analyse à Eléments Finis. Dans ce travail, les bases d'une approche faiblement intrusive, la méthode de Substitution, sont présentés dans le domaine de la dynamique explicite. Il s'agît d'une formulation globale-locale, conçue pour appliquer un modèle grossier sur tout le domaine pour obtenir une réponse globale: ce pré-calcul est ensuite corrigé itérativement par l'application du modèle raffiné appliqué seulement où nécessaire. La vérification de la méthode de Substitution en comparaison avec la méthode de Décomposition de Domaines est présentée.