Ce travail de thèse est consacré à l'étude des phénomènes d'initiation et de propagation d'une fissure dans un matériau bitumineux. À l'aide d'un essai innovant développé au Laboratoire Central des Ponts et Chaussées, on reproduit les conditions d'initiation et de propagation d'une fissure confinée dans un film de bitume. Dans la première partie de ce travail de thèse, on aborde l'approche utilisée pour identifier les propriétés rhéologiques à l'aide du modèle de Maxwell généralisé. La détermination des séries de Prony passe par les méthodes d'optimisation linéaire et non-linéaire pour représenter précisément le module complexe du liant bitumineux dans le domaine fréquentiel. Ce comportement viscoélastique identifié est ensuite validé par l'essai de

Dans ce travail de thèse, on s'intéresse à la réponse vibratoire d'un rotor affecté de fissures afin d'exhiber les paramètres vibratoires susceptibles de favoriser leur détection dans le cadre d'une procédure de surveillance. Dans un premier temps, les travaux expérimentaux, numériques et analytiques antérieurs effectués visant la modélisation de telles structures sont résumés et commentés. Puis la démarche de la prise en compte des effets de cisaillement dans la loi de comportement en flexion d'une section de rotor fissurée est présentée et est validée par des calculs tridimensionnels tenant compte du contact unilatéral entre les lèvres de fissure. La démarche est basée sur une formulation énergétique du problème et peut être appliquée à toute géométrie de fissure. L'exploration par voie numérique de la réponse d'un rotor fissuré est présentée dans le troisième chapitre de cette thèse. Différents systèmes dynamiques sont considérés. Les parties saines du rotor sont représentées part des éléments de type barre ou poutre, et la section fissurée par un ressort non-linéaire pour tenir compte du mécanisme de respiration de la fissure. Vers la fin de cette partie, on présente une méthode originale de construction d'un élément fini de poutre comportant des fissures. Le dernier chapitre est consacré à l'étude analytique du système à deux ddl. Ainsi, plusieurs fonctions de respiration de la fissure sont examinées. La stabilité linéaire des solutions périodiques est étudiée vers la fin de ce chapitre par la théorie de Floquet. Finalement quelques paramètres vibratoires sont proposés comme indicateurs de la présence de fissures.

Cette thèse comporte deux parties. L'objectif de la première partie était l'implémentation d'une méthode de résolution de l'équation dynamique de Landau-Lifshitz-Gilbert fondée sur l'approximation des éléments finis. Deux approches ont été présentées: la première reposant sur une formulation avec contrainte et la seconde mettant en œuvre des fonctions tests dans le plan tangent à l'aimantation. Seule la seconde approche reproduit la dynamique obtenue en différences finies sur des cas tests. Dans la seconde partie, le but était d'étudier le déplacement de parois de domaines magnétiques sous l'effet d'un champ magnétique ou d'un courant polarisé en spin dans des systèmes à anisotropie perpendiculaire. Il a été nécessaire d'introduire dans l'équation dynamique les termes dus au transfert de spin. Des systèmes idéaux et des systèmes avec différents types de défauts ont été étudiés. Les résultats numériques ont été comparés avec les données expérimentales disponibles.

Ce travail concerne l'étude des machines asynchrones à cage dans les applications de traction automobile, en particulier pour les véhicules hybrides. Nous avons développé des modèles et des méthodes utiles pour analyser et dimensionner de telles machines électriques. Nous avons tout d'abord mis au point un modèle électromagnétique non linéaire de la machine, déterminé à partir d'un nombre restreint de calculs "éléments finis". Ce modèle a ensuite été adapté pour réaliser des études de sensibilités sur quelques dimensions géométriques importantes. Il permet d'adapter rapidement une machine à un nouveau cahier des charges. Nous avons finalement étendu cette méthode de modélisation pour prendre en compte un plus grand nombre de paramètres géométriques. Ce dernier modèle paramétrique a été utilisé pour l'optimisation sous contrainte des dimensions d'une machine. Pour cela, nous avons proposé une nouvelle méthode "d'optimisation à modèle recalé" qui concilie précision, rapidité et simplicité. Cette démarche a été appliquée au cas concret d'un dimensionnement de machine asynchrone avec un cahier des charges typique d'un véhicule hybride.

SALOME est une plate-forme générique Open Source de Pré-Post traitement destinée à être spécialisée pour y intégrer des codes de calcul existants. Elle intègre plusieurs modules dans une architecture de composants distribués. Notre mission dans le cadre du projet SALOME a été la conception et la mise en œuvre technologique du module DATA, dédié à la description des propriétés physiques.En ce sens, nous avons réalisé un nouveau langage dédié à la description des modèles de données physiques : le SPML (SALOME Physics Modelling Language). Un métamodèle dédié à la description des propriétés physiques offre au langage SPML sa base sémantique. Pour réaliser une partie commune de communication entre des modèles représentant différents domaines de la physique, il a été développé un Modèle de Données Commun, matérialise par une librairie SPML réutilisable. La réalisation de l'IHM, notamment l'adaptation automatique de l'interface graphique aux modèles physiques décrits en SPML, font du module DATA un outil performant, qui permet une adaptation facile de la plate-forme, pour tout domaine de la physique. Nous avons réalisé une première connexion du solveur Flux, dans le cadre de la plate-forme pour des analyses magnétostatiques sur des problèmes décrits en SALOME.

Le chargement dynamique des paliers hydrodynamiques et des amortisseurs à film fluide soumet le film lubrifiant à une action d'écrasement périodique. Cette sollicitation conditionne la portance du film lubrifiant et peut engendrer sa rupture, créant ainsi des conditions sévères de fonctionnement. L'étude approfondie, tant expérimentale que numérique du comportement du film soumis à une action d'écrasement périodique, s'avère donc indispensable pour pouvoir maîtriser les performances de ces dispositifs. <br /><br />Pour ce faire, trois codes de calcul numérique sont développés dans le cadre de cette thèse. Ces codes sont basés sur la théorie de JFO. Les équations établies sont discrétisées en différences finies et en éléments finis. Pour un film d'huile écrasé périodiquement, à faibles fréquences, entre deux surfaces conformes, les prédictions numériques de ces trois codes sont très proches. Par ailleurs, les visualisations numériques de la rupture dans les différents types de contacts ont montré l'impact de la fréquence des oscillations sur l'étendue et la durée de la rupture de film. <br /><br />L'étude expérimentale est conduite via un banc d'essai conçu puis réalisé. Deux modes de rupture de film sont identifiés : la rupture par séparation des gaz saturés et la rupture par séparation due à l'introduction de l'air ambiant. La visualisation de chaque mode est conditionnée par la fréquence des oscillations. La forme de la zone de rupture due à l'introduction de l'air ambiant dépend de la fréquence, de la hauteur du bain d'huile et de la géométrie du contact. <br /><br />Pour les surfaces non conformes, les résultats des modèles numériques sont en bon accord avec les résultats expérimentaux. Les différences constatées dans le cas des surfaces conformes sont dues aux imperfections géométriques et à la souplesse des structures qui composent le dispositif expérimental.

Nous nous intéressons au cas de structures poutres vibrantes soumises à des précharges. Ces précharges influent sur les fréquences propres. Elles induisent un comportement non linéaire de la structure. L'attention est portée en particulier sur les poutres préchargées en rotation uniforme. Celles-ci subissent alors les effets des forces d'inertie centrifuge et aérodynamiques. Les effets rigidifiants augmentent la raideur en flexion des poutres et décalent les modes propres. Il apparaît alors une matrice de raideur géométrique, et une matrice de raideur cinétique dans le cas des poutres en rotation. Nous utilisons une méthode développée par D. GAY permettant une modélisation unidimensionnelle d'une structure réelle composite tridimensionnelle. Cette méthode s'appuie sur le concept d'homogénéisation avec prise en compte des effets de gauchissement. Des paramètres intégraux de déplacements et de rotation «moyens» de la section permettent alors la constitution d'un élément fini de poutre. Nous avons rédigé et mis au point le code de calculs correspondant, qui permet d'effectuer des calculs de fréquences et modes de poutres composites à ligne moyenne rectiligne avec des conditions aux limites variées en déplacement et en efforts. Afin de valider les outils informatiques développés, nous avons réalisé des vérifications expérimentales sur des poutres soumises à des précontraintes de traction et de flexion pure. Ces essais, menés conjointement aux calculs numériques ont montré une bonne corrélation

Ces travaux portent sur le développement d'un capteur capacitif de pluie à électrodes protégées fonctionnant en mode différentiel. Dans un premier temps, les propriétés électriques des différents matériaux constitutifs du capteur sont déterminées expérimentalement. A partir de considérations théoriques, nous avons établi des expressions semi-empiriques traduisant la dépendance en fréquence [10^2 à 10^5Hz] et en température [−40 à +80°C] des permittivités diélectriques, nécessaires aux calculs des capacités. Le capteur est ensuite dimensionné par simulation numérique (méthode des éléments finis) puis fabriqué selon différents procédés de dépôt pour les électrodes. Deux motifs d'électrodes sont étudiés : électrodes linéaires et électrodes interdigitées. La sensibilité à l'eau est évaluée par simulation numérique (film puis gouttes d'eau de volume variable) avant validation expérimentale. L'influence de la température sur les performances du capteur est finalement abordée : les effets de la dilatation thermique des matériaux et de la variation des permittivités diélectriques sont successivement traités. A l'issue de cette étude, l'efficacité du mode différentiel contre l'influence de la température est évaluée.

Cette thèse porte sur le développement de méthodes numériques pour le calcul d'écoulements incompressibles multi-fluides sur grille de calcul, s'inscrivant ainsi dans le cadre du projet MecaGrid.Une étude de la grille MecaGrid met en évidence son caractère hétérogène et ses conséquences. Plusieurs techniques d'optimisation sont présentées afin d'améliorer son utilisation : répartir la masse de calcul de façon adaptée, et privilégier le travail des processeurs vis-à-vis des communications réseaux pénalisantes.Des méthodes numériques sur maillage non structuré composé de tétraèdres sont choisies pour réaliser la simulation directe d'écoulements multi-fluides avec capture d'interface. Nous adoptons une approche unique qui rappelle celle du Multiscale, dans laquelle la condensation d'une fonction Bulle pyramidale est utilisée comme technique universelle de stabilisation.Les équations de Navier-Stokes incompressible sont résolues par une méthode éléments finis mixtes à interpolation P1+/P1. Un schéma temporel d'Euler implicite est appliqué, en association avec un algorithme de Newton pour linéariser le problème.Les interfaces sont capturées par une technique Level Set à interpolation continue P1 qui consiste à résoudre une équation de transport stabilisée par la condensation d'une Bulle (Residual-Free Bubbles). Le couplage avec une équation d'Hamilton-Jacobi permet de réinitialiser la fonction Level Set au court de son transport. La comparaison avec une méthode Galerkin discontinue proche du Volume of Fluid montre que le Level Set se distingue par sa simplicité et l'absence de diffusion numérique.Enfin, les simulations numériques sont validées par plusieurs cas test reconnus.

Dans le cadre de la conception des structures aéronautiques et de la modélisation de leur tenue mécanique, la représentation des comportements non linéaires locaux conduisant à la rupture est nécessaire et notamment ceux des liaisons. Ainsi en situation de crash survivable, les liaisons sont soumises à des charges dynamiques qui peuvent conduire à la rupture et à la désarticulation. Pour caractériser la tenue au crash des assemblages rivetés, une méthodologie dite de 'plans d'expérience numériques' est proposée ; son intérêt est de limiter le coût des procédures expérimentales. Fondée sur les techniques de modélisation par éléments finis et d'optimisation par méthode inverse, la méthode repose sur la détermination des paramètres du comportement non linéaire et à la rupture des matériaux constitutifs de la liaison rivetée. Le développement de l'outil informatique d'optimisation est présenté dans la première partie du mémoire. Les modèles de comportement de deux alliages d'aluminium composant l'assemblage sont identifiés grâce à l'optimiseur. Ils seront utilisés pour mener les simulations de différents essais effectués sur les assemblages. L'évaluation des contraintes et des déformations résiduelles dues au procédé de rivetage est abordée dans la deuxième partie. Les mécanismes de mise en forme des rivets sont analysés expérimentalement. Des simulations sont d'autre part entreprises à l'aide du code de calcul explicite PAM-SOLID(TM), afin d'introduire, une fois validées, les contraintes et les déformations résiduelles post-rivetage dans le calcul de la tenue mécanique de l'assemblage. La faisabilité des plans d'expérience numériques est démontrée dans la troisième partie. Les paramètres d'endommagement du modèle de Gurson sont identifiés par méthode inverse pour les deux matériaux composant l'assemblage. La qualité de l'ensemble des résultats obtenus permet notamment d'envisager l'emploi de l'outil de simulation dans le cadre de la caractérisation de modèles équivalents macroscopiques du rivet. Ces techniques de modélisation simplifiée des liaisons sont abordées dans la dernière partie. Un critère de rupture de la liaison est caractérisé par l'expérience à l'aide du dispositif ARCAN et par plan d'expérience numérique. Un nouveau modèle non linéaire de la liaison est développé dans le code de calcul PAM-SOLID(TM) et ses paramètres identifiés. L'ensemble des résultats obtenus est finalement appliqué à la simulation de l'écrasement d'une structure comportant 700 rivets.