.

Cette étude présente des algorithmes de génération et d'améliorations de maillages en simplexes pour la methode des éléments finis (triangles en 2d et tetraedres en 3D). Dans une première partie, après avoir montre que tout domaine plan défini par un contour extérieur et un ensemble de points internes représentant sa discrétisation peut être triangule, on expose et on justifie un algorithme simple de maillages en triangles. On prouve ensuite que l'opération de changement de diagonales dans un quadrilatère convexe permet d'atteindre toute les topologies realisables. Dans la seconde partie, on étudie le cas tridimensionnel et on donne quelques propriétés sur l'existence d'une tetraedrisation quand on se fixe le maillage en facettes triangulaires de la surface extérieure du domaine et les points internes pour la discrétisation. On décrit plusieurs heuristiques d'améliorations topologiques de maillages tetraedriques et on compare leur efficacité pratique sur quelques exemples industriels tires de la mécanique des fluides: on examine ainsi un algorithme glouton généralisant la démarche bidimensionnelle, un algorithme base sur la methode statistique du recuit simule et des algorithmes de type Divide and Conquer couples avec des procédures non topologiques simples comme le barycentrage des Nuds d'une discrétisation. Ces derniers algorithmes provoquent des améliorations conséquentes même sur des maillages de topologies très défavorables

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é.

Dans la première partie de ce travail on présente des méthodes numériques pour résoudre l'équation de convection-diffusion et on propose une méthode qui ramène le problème à une formulation symétrique qui permet de traiter les cas de faible viscosité. Les problèmes unidimensionnel, bidimensionnel stationnaires et unidimensionnel non stationnaire sont abordes. La deuxième partie est consacrée à la résolution d'un modèle test de circulation océanique aux latitudes moyennes dont le mouvement est du au vent. On introduit une condition d'adhérence aux parois est et ouest qui représentent les cotes terrestres. On utilise une méthode d'éléments finis qui permet de résoudre le modèle comme une "cascade" des problèmes très simples. On compare avec des méthodes de différences finies du point de vue de la précision, du cout de calcul et de l'applicabilité.

.

.

On rappelle les notions concernant le transducteur général d'états finis du systeme PIAF (programmes interactifs de l'analyse du français) et on décrit certaines applications réalisées à partir du transducteur général d'états finis. On l'étend à un réseau d'automates appelé MULTIPIAF permettant d'envisager de nouvelles applications.

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.

Avec le vieillissement de la population, l'ostéoporose et les fractures de l'extrémité supérieure du fémur sont devenues un problème de santé publique majeur. L'intervention chirurgicale, à laquelle s'ajoutent la rééducation des patients et la gestion de leur dépendance, représentent un poids socio-économique considérable qui doit être réduit, sans pour autant proposer un traitement préventif systématique, lui aussi coûteux. La mesure de densité minérale osseuse par DXA fait aujourd'hui référence dans l'évaluation du risque de fracture mais cet outil reste limité par le chevauchement des valeurs mesurées entre les sujets fracturés et non fracturés. Cette étude propose de développer des modèles par éléments finis 3D du fémur proximal permettant d'améliorer la prédiction du risque de fracture. Un modèle anatomique personnalisé tenant compte des spécificités géométriques et mécaniques complexes de chaque fémur a conduit, sur 39 spécimens, à un bon compromis entre une prédiction individuelle de l'effort à rupture (Fexp = 1,006 FMEF avec r = 0,87, IC95% = +- 2600 N) et un coût numérique adapté à une future exploitation clinique de cet outil (simulation numérique 30 mn). Afin de consolider la personnalisation, une étude expérimentale menée sur 92 éprouvettes (46 traction et 46 compression) a permis d'améliorer la prédiction des propriétés mécaniques de l'os cortical à partir de l'imagerie médicale (r = 0,42 à 0,72). Des essais in-vitro macroscopiques sur 40 paires de fémurs ont également été réalisés afin de constituer une base de données solide pour l'évaluation du modèle. Parallèlement, un modèle paramétré a été utilisé dans le cadre d'une étude par plans d'expérience afin d'éclaircir le rôle des paramètres géométriques dans la résistance mécanique de l'extrémité supérieure du fémur. Ces travaux ont souligné le potentiel de la modélisation par éléments finis personnalisée à devenir un outil complémentaire de prédiction du risque, et ont mis en évidence le rôle prépondérant de l'enveloppe corticale et de sa limite d'endommagement.

Les polymères amorphes manifestent un comportement viscoélasto-viscoplastique non linéaire qui dépend fortement de la température et de la vitesse de déformation. Ce comportement est décrit au moyen du modèle des défauts quasi ponctuels présentant un fort contenu physique. En effet, ce modèle distingue les parts élastique, viscoélastique et viscoplastique de la déformation telles qu'elles peuvent être mises en évidence expérimentalement. L'expression des lois cinétiques relatives au développement des déformations non élastiques repose sur des concepts de mobilité moléculaire. Le modèle rend compte de manière cohérente de la réponse du matériau, en petites et grandes déformations, à des sollicitations mécaniques sur une large plage de vitesse de déformation et de température. Conçue sur la base de ce modèle de mobilité moléculaire, une loi de comportement a été établie sous forme tensorielle et implémentée dans le code de calcul de structures ABAQUS au moyen d'une routine utilisateur umat. <br />Un jeu de paramètres matériaux a été déterminé pour le polycarbonate Bisphénol A Lexan. Des essais de traction, de compression et de torsion ont été réalisés et simulés au moyen de l'outil développé. La validation a porté, d'un point de vue global, sur la confrontation des réponses calculée et mesurée de la structure en termes de courbes force-déplacement. A un niveau local, elle a consisté à vérifier l'aptitude de la simulation à rendre compte d'hétérogénéités de déformation associées à des concentrations de contrainte et au développement d'une déformation en tonneau d'un plot de compression en présence de frottement au contact. La simulation a également montré l'aptitude du modèle éléments finis à décrire la propagation de la striction au cours d'un essai de traction.