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

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.

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 présentons une étude numérique de la localisation des déformations en bande de cisaillement dans les géomatériaux.<br />Nous utilisons les modèles de second gradient locaux: sur un problème de cisaillement dans un milieu unidimensionnel, et de simulation de creusement de puits bidimensionnels.<br /><br /> Dans un premier temps, nous avons réalisé une analyse de problème de cisaillement unidimensionnel avec la mécanique des milieux continus classique et pour un comportement quasi-fragile élasto-plastique. Afin de jouer des calculs non restreintes, nous avons employé la méthode de longeur d'arc dans un code aux éléments finis unidimensionnel pour un milieu de second gradient. Nous avons choisi de réaliser notre expérience numérique en utilisant une loi de comportement plus réaliste.<br /><br />Le modèle de second gradient bidimensionnel a été utilisé pour réaliser des simulations numérique d'excavation avec une loi de comportement quasi-fragile. Nous avons réalisé la simulation sur un puits complet en utilisant le code aux éléments finis en grandes déformation Lagamine.<br /><br />Les résultats obtenus sont significatifs. Ils ont illustré la non unicité des solutions, dont différents modes localisés ont été obtenus, qualitativement semblable à la rupture de puits observés dans les laboratoires.

Ce mémoire concerne l'étude et la réalisation d'un code d'éléments finis pour les circulations océaniques engendrées par les vents dans un bassin ferme aux latitudes moyennes. Deux modèles physiques sont examines: 1) la modélisation barotrope océan de densité constante proposée en 1982 par Dumas-Le Provost-Poncet; 2) une extension de cette modélisation a un océan stratifie en densité. Dans les deux cas, une étude théorique démontre l'existence et l'unicité de la solution du système d'équations traite. La methode des éléments finis autorise des domaines a géométrie complexe, ce qui a permis d'expérimenter le modèle barotrope sur un bassin réaliste: l'atlantique nord. Le code de modélisation des circulations barotrope et barocline, écrit pour un ordinateur vectoriel, est teste sur un bassin carre découpé en trois couches de densité variable selon la verticale. Les résultats des simulations effectuées sont compares a ceux obtenus avec un code aux différences finies. Cette modélisation multicouche fournit un système d'équations se prêtant naturellement a l'utilisation d'algorithmes parallèles, un code de calcul a été étudié et implante sur l'hypercube t20 de la société FPS

Dans le domaine de la simulation numérique en électromagnétisme, l'évaluation d'erreur en régime magnétodynamique n'a pas trouvée de réponse pleinement satisfaisante jusqu'à maintenant. Après avoir introduit la méthode des éléments finis à partir d'une approche thermodynamique de l'électromagnétisme, nous proposons un nouveau critère d'erreur énergétique, qui permet d'évaluer l'erreur dans les parties conductrices sièges de courants de Foucault. Les caractéristiques de ce critère en fonction de la formulation et de l'ordre des éléments sont détaillées. Nous présentons ensuite une procédure de maillage adaptatif reposant sur ce critère appliquée à un cas en régime magnéto-harmonique 2D. Une extension aux régimes transitoires est introduite. La convergence des grandeurs énergétiques globales valide la pertinence du critère d'erreur proposé.

L'objectif de ce travail consiste à étudier numériquement la propagation des ondes mécaniques dans un milieu poreux continubiphasique en régime impulsionnel. Un code de calcul permettant la simulation et l'analyse du comportement mécanique est mis au point à cet effet. Une analyse préliminaire présente des généralités sur la notion de milieu poroviscoélastique, sur la théorie de Biot ainsi que sur la propagation des ondes mécaniques. Un bilan des résultats rencontrés dans la bibliographie montre la nécessité de mener une étude complémentaire sur le modèle complet du milieu poreux continu biphasique dans le domaine temporel.<br /><br />Une approche par éléments finis est proposée dans le cadre de la théorie générale de Biot. Les caractéristiques de l'outil numérique développé sont précisées. En particulier, la structure orientée objet donne un code compact et souple. Un travail<br />semi-analytique, préalablement effectué, s'intéresse aux phénomènes de dispersion, d'atténuation et à la détermination des vitesses de propagation des différentes ondes.<br /><br />Une modélisation bidimensionnelle permet d'obtenir les déplacements temporels des phases solide et fluide en surface et en profondeur d'un sol poreux semi-infini. Une étude paramétrique des couplages mécaniques est effectuée. La seconde onde de compression est mise en évidence. Une première approche de sols hétérogènes ou partiellement saturés est en outre proposée.<br /><br />L'étude de problèmes tridimensionnels est ensuite envisagée. La taille numérique importante de cette problématique nécessite alors la parallélisation du code de calcul. Des essais sur différents supercalculateurs sont réalisés pour mesurer la performance du<br />calcul parallèle et conduisent à des résultats tridimensionnels.

Les paliers soumis à des charges dynamiques, comme par exemple les paliers de moteurs, ont une géométrie simple mais un fonctionnement assez complexe. La modélisation de tous les phénomènes physiques qui interviennent dans le fonctionnement d'un palier dynamique est très difficile à réaliser et nécessite parfois des temps de calcul prohibitifs. Le but de cette étude est d'aider les concepteurs des paliers en proposant des outils prédictifs rapides et performants. <br />Avec les lois de variation de la viscosité adaptées, les influences respectives des effets non-newtonien et piezo-visqueux sont étudiées dans le cas complexe des paliers de tête de bielle. Pour les cas étudiés l'effet piezo-visqueux s'avère être plus significatif sur le comportement du palier, que l'effet non-newtonien. Deux modèles thermiques sont développés : un premier modèle, qui donne une information 2D de la température dans le contact et un second, plus évolué, qui permet d'obtenir une information 3D de la température. L'efficacité des modèles est prouvée par le traitement d'un palier de tête de bielle. Une comparaison des résultats en utilisant un modèle EHD, deux modèles TEHD 3D avec et sans les effets non-newtoniens et un modèle TEHD dit "global" est présentée et analysée. <br />Afin de valider les modèles théoriques et numériques une machine d'essais (MEGAPASCALE) pour étudier la lubrification des paliers de bielle dans des conditions de fonctionnement réelles et sévères est étudiée et réalisée. Les diagrammes de charge obtenus (90kN en compression et 60 kN en traction pour une vitesse de rotation de 20000 tr/min) sont des diagrammes essentiellement "inertiels" comme ceux rencontrés dans des moteurs à hauts régimes. Une description détaillée des nombreuses solutions originales adoptées ainsi que le fonctionnement du banc est présentée.

Dans les accidents de la route, les lésions cervicales représentent 15% des lésions traumatiques graves et plus de 50% des demandes de remboursements. Les mécanismes lésionnels associés sont encore mal connus, et médecins et assureurs sont demandeurs de critères quantifiant les lésions de type whiplash. Le but de cette étude est de contribuer à l'amélioration de modèles Eléments-Finis permettant d'explorer les mécanismes lésionnels survenant lors de ces sollicitations. A cet effet des mécanismes lésionnels probables ont été identifiés et comparés aux résultats de simulations numériques. L'influence de paramètres extrinsèques et intrinsèques sur le risque lésionnel associé a été évaluée, permettant de mettre en relief l'intérêt d'une étude de critères locaux dans ce contexte. En particulier, l'étude de l'influence de la morphologie a été accompagnée d'une évolution des modèles vers la personnalisation, basée sur la reconstruction géométrique d'un sujet à partir de stéréo-radiographies.

Les traumatismes subis à la face lors d'accidents sur la voie publique, d'activités sportives ou bien de rixes ne nuisent généralement pas à la sur-vie de l'individu mais peuvent requérir des soins coûteux et laisser des séquelles esthétiques importantes. Par ailleurs, la face joue un rôle protecteur du contenu intracrânien et sa réponse au choc conditionne les sollicitations appliquées à celui-ci.<br />En vue de prédire précisément les risques de blessures de la face et du contenu intracrânien, l'objectif de cette thèse est de développer un modèle en éléments finis de la structure osseuse de la face et du crâne pour la simulation de sa réponse au choc.<br />Une première partie du travail a consisté à construire le maillage d'une structure osseuse cranio-faciale à partir de coupes scanner de tête. Le choix d'un maillage en éléments de type plaque, de densité suffisante pour représenter fidèlement la géométrie complexe de cette structure osseuse et d'épaisseur variable, a été fait.<br />L'hypothèse d'un matériau osseux homogène et isotrope a été choisie pour l'ensemble de la structure cranio-faciale. Des essais de flexion statique sur des échantillons d'os crânien, associés à leur simulation numé-rique et une méthode d'identification, ont permis de définir les propriétés élasto-plastiques de ce matériau. Les résultats se situent correctement par rapport aux intervalles de valeurs de la littérature.<br />La réponse du modèle cranio-facial a été validée sous sollicita-tions statiques. Pour cela des essais spécifiques de compression de la face ont été réalisés sur pièce anatomique. Les courbes globales effort - déplacement expérimentale et numérique ont été comparées pour « calibrer » les propriétés du matériau. La réponse du modèle a ensuite été validée par comparaison du champ de déplacement mesuré expérimentalement par une méthode de corrélation d'images et obtenu par simulation.<br />Le modèle ainsi validé en statique a été évalué sous sollicitations dynamiques. Sa réponse au choc a été comparée aux résultats d'un impact sur la face réalisé spécifiquement. La réponse au choc globale du modèle est similaire à celle enregistrée expérimentalement.<br />Plusieurs perspectives d'exploitation de ce modèle sont envisa-geables dans le domaine du choc ou celui de la chirurgie. En particulier, il aidera à définir des critères de blessures de la tête en cas de choc sur la face.