Cette thèse porte sur l'étude d'une méthode de type Galerkin discontinu en domaine temporel (GDDT), afin de résoudre numériquement les équations de Maxwell instationnaires sur des maillages hybrides tétraédriques/hexaédriques en 3D (triangulaires/quadrangulaires en 2D) et non-conformes, que l'on note méthode GDDT-PpQk. Comme dans différents travaux déjà réalisés sur plusieurs méthodes hybrides (par exemple des combinaisons entre des méthodes Volumes Finis et Différences Finies, Éléments Finis et Différences Finies, etc.), notre objectif principal est de mailler des objets ayant une géométrie complexe à l'aide de tétraèdres, pour obtenir une précision optimale, et de mailler le reste du domaine (le vide environnant) à l'aide d'hexaèdres impliquant un gain en terme de mémoire et de temps de calcul. Dans la méthode GDDT considérée, nous utilisons des schémas de discrétisation spatiale basés sur une interpolation polynomiale nodale, d'ordre arbitraire, pour approximer le champ électromagnétique. Nous utilisons un flux centré pour approcher les intégrales de surface et un schéma d'intégration en temps de type saute-mouton d'ordre deux ou d'ordre quatre. Après avoir introduit le contexte historique et physique des équations de Maxwell, nous présentons les étapes détaillées de la méthode GDDT-PpQk. Nous réalisons ensuite une analyse de stabilité L2 théorique, en montrant que cette méthode conserve une énergie discrète et en exhibant une condition suffisante de stabilité de type CFL sur le pas de temps, ainsi que l'analyse de convergence en h (théorique également), conduisant à un estimateur d'erreur a-priori. Ensuite, nous menons une étude numérique complète en 2D (ondes TMz), pour différents cas tests, des maillages hybrides et non-conformes, et pour des milieux de propagation homogènes ou hétérogènes. Nous faisons enfin de même pour la mise en oeuvre en 3D, avec des simulations réalistes, comme par exemple la propagation d'une onde électromagnétique dans un modèle hétérogène de tête humaine. Nous montrons alors la cohérence entre les résultats mathématiques et numériques de cette méthode GDDT-PpQk, ainsi que ses apports en termes de précision et de temps de calcul.

La Fabrication Directe par Projection Laser (FDPL) est un procédé novateur qui permet la réalisation de pièces de géométrie complexe à partir d'un jet de poudre. Les pièces sont obtenues en superposant différentes couches de matière à l'aide d'une buse coaxiale avec un faisceau laser. Cette technique de fabrication, très rapide, reste cependant peu répandue, en particulier du fait de l'irrégularité des surfaces. L'objet de ce travail est de développer des modèles numériques capables de prédire le paramètre d'ondulations caractérisant l'état de surface à partir uniquement des paramètres opératoires et des propriétés thermophysiques du substrat et de la poudre, dans le but de mieux comprendre les mécanismes responsables de l'état de surface dégradé. Une première étude a porté sur la fusion locale d'un barreau vertical sous l'action d'un laser afin de valider la mise en données du modèle numérique. Ce modèle thermohydraulique inclut les effets de tension superficielle et utilise la méthode ALE (Arbitrary Lagrangian-Eulerian) pour suivre la déformation de la surface libre. Ce modèle prédictif a été transposé au procédé FDPL en intégrant l'apport de matière dû au jet de poudre. Les caractéristiques de ce jet de poudre sont obtenues grâce à un modèle 3D prédisant la trajectoire des particules et leur échauffement sous l'action du faisceau laser. Ces données servent à définir les conditions aux limites de modèles 2D thermohydrauliques du bain fondu. Nous avons ainsi retrouvé, pour un substrat mince, la corrélation établie expérimentalement entre l'amplitude des ondulations et le taux de dilution. Il a été montré que, dans le cas de l'alliage de titane Ti-6Al-4V, un meilleur état de surface est obtenu pour une forte puissance laser, une vitesse de défilement élevée et un faible débit massique. La simulation 2D de dépôts multicouches a permis d'étudier l'influence de la stratégie de balayage et du temps de pause entre couches. Une étude de sensibilité a également été menée afin d'analyser le rôle de paramètres tels que le coefficient thermocapillaire ou la viscosité. L'étude a par la suite été étendue à un cas 3D pour étudier l'effet de la distribution énergétique du faisceau laser ainsi que les propriétés thermophysiques. La prédiction des différents modèles numériques développés à l'aide du code de calcul COMSOL Multiphysics® est validée grâce à des données expérimentales.

La simulation numérique de comportement des futurs produits est largement utilisée sur les modèles virtuels de produits avant leur fabrication physique. Cependant, le processus pourrait encore être optimisé en particulier pendant la phase d'optimisation du comportement de produit. Ce processus implique la répétition de quatre étapes principales de traitement : conception de CAO, création de maillage, enrichissement de sémantique par la modélisation du comportement physique et enfin calcul par éléments finis (EF). L'analyse de comportement de produit est effectuée à partir de la première solution de conception puis sur les nombreuses boucles successives d'optimisation de produit. Chaque évaluation de solution nécessite le même volume de temps que celui nécessaire pour la première conception de produit, cela est particulièrement crucial dans le contexte de maintenance de produit et d'évaluation de cycle de vie de produit. Cette thèse propose un nouveau cadre de travail pour l'optimisation de produit à partir de simulation par EF menées successivement sans retour à la CAO initiale du produit, ce qui réduit les activités de préparation de maillages et d'enrichissement sémantique E.F. Plus concrètement, l'idée est d'opérer directement le maillage enrichi par la sémantique E.F pour optimiser le produit. Dans cette thèse, les concepts sous-jacents et les composants conçus pour le développement de ces opérateurs de modification sont présentés et analysés. Une spécification d'opérateur de haut niveau est proposée selon une structure modulaire qui permet ensuite une réalisation facile des différents opérateurs de modification de maillage. Enfin, quatre déclinaisons de cet opérateur de maillage de haut niveau sont présentées: la fusion, la fissuration, le perçage et le congé d'arête. Ces opérateurs ont été prototypés et validés sur des modèles E.F. académiques et industriels, permettant de démontrer leur efficacité et la pertinence de l'approche proposée.

L'animation de personnages virtuels guidée par des données fait l'objet de nombreuses études dans le domaine de l'informatique graphique. Dans ce contexte, le mouvement est classiquement défini par une suite de squelettes au cours du temps, chacun étant décrit par un vecteur de positions et de rotations. Le maillage 3D est ensuite guidé par les squelettes au moyen d'un couplage défini explicitement. L'enchaînement des différentes étapes de ce processus est difficile à mettre en oeuvre, et conduit à des approximations et des erreurs de modélisation, à la fois au niveau de l'animation du squelette et de la déformation du maillage. Dans ce manuscrit nous proposons d'étudier d'autres représentations du mouvement par le biais d'un ensemble de relations spatiales. Cette approche nous permet de tenir compte implicitement de contraintes de distance entre les points de la structure articulée, du maillage et de l'environnement, et d'exprimer en particulier la notion de contact. Deux axes d'étude sont principalement abordés~: le premier considère le mouvement dans l'espace métrique, et le second caractérise chaque posture par son information différentielle dans l'espace Laplacien. Dans un premier temps, nous proposons de représenter les squelettes associés aux postures du mouvement par un ensemble de distances. Caractériser une telle structure dans l'espace métrique se ramène à un formalisme mathématique connu sous le nom de problème de la géométrie des distances. Ainsi, nous nous inspirons des différentes techniques existantes et les appliquons au contrôle du mouvement. L'objectif est de produire de nouveaux mouvements à partir de processus d'édition ou d'inversion cinématique. Nous montrons que cette représentation permet un contrôle simple et intuitif de l'animation d'un personnage. Elle possède également plusieurs propriétés exploitables dans le cadre de l'analyse du mouvement. Ce dernier point est illustré par une application originale de recherche de mouvements dans des grandes bases de données. Dans un second temps, nous définissons le mouvement par un ensemble de graphes dont les sommets sont caractérisés par une information différentielle. A travers cette représentation, nous proposons une nouvelle méthode d'édition du mouvement couplant des contraintes de distance avec l'opérateur Laplacien discret. Cet opérateur permet de préserver les relations spatiales lors de l'édition du mouvement alors que les contraintes de distance préservent certaines propriétés inhérentes au squelette. Ce concept donne lieu à plusieurs applications dédiées à la reconstruction et l'édition de mouvement : (i) l'édition interactive d'animation de squelette, où nous proposons d'éditer tout type de mouvement avec de fortes déformations tout en préservant l'information spatio-temporelle ; (ii) la reconstruction de trajectoires de marqueurs~: en faisant l'hypothèse qu'il existe un lien entre la trajectoire d'un marqueur et celles de ses voisins, nous proposons de reconstruire les trajectoires incomplètes ; (iii) l'animation de maillage où nous proposons un nouveau processus d'animation directement guidé par les trajectoires des marqueurs.

Avec le développement de l’internet des objets et l’augmentation des applications sans fil, les antennes sont de plus en plus présentes au quotidien. Cependant, l’implantation de ces antennes est un challenge tant d’un point de vue technologique (intégration des antennes dans les dispositifs), que psychologique (acceptabilité des antennes par le grand public). Dans ce contexte, le développement d’antennes optiquement transparentes permet non seulement leur implantation sur de nouvelles surfaces (vitrages d’immeubles, écrans de smartphones ...), mais promeut aussi leur acceptabilité par le grand public grâce à leur faible impact visuel. Ce travail présente la conception, la fabrication et la caractérisation d’antennes actives optiquement transparentes. Le matériau transparent et conducteur utilisé est un maillage métallique à pas micrométrique développé spécifiquement, alliant conductivité électrique et transparence optique élevées. Dans ce cadre, un premier prototype d’antenne transparente et miniature en bande FM utilisant un transistor MESFET de dimensions sub-millimétriques a été réalisé. Des antennes agiles en fréquence en bande X (~10 GHz) couplées, soit à une diode varicap localisée (agilité ~10%), soit à un matériau ferroélectrique (agilité ~2%), ont été développées et étudiées. Une antenne passive transparente a été conçue en bande V (~60 GHz). Enfin, une transition optique (1540 nm) / hyperfréquence (1,4 GHz) a été réalisée et caractérisée, basée sur la transmission optique d’un faisceau laser au travers du matériau constitutif de l’antenne. Pour l’ensemble des prototypes réalisés, une transparence optique supérieure à 80% dans le domaine du visible associée à une résistance par carré inférieure à 0,1 ohm/sq ont été utilisées. / Within the development of the Internet of Things (IoT) and the increase of the wireless communications, antennas are even more present on everyday life. However, antenna implementation is a real challenge, from a technological point of view (antenna integration into the devices) and from a psychological point of view (acceptability by the general public). Within this framework, the development of optically transparent antennas on new surfaces (glass windows, smartphone screens . . . ) is of great interest to improve the network coverage and to assist the general public in acceptability thanks to the low visual impact of such printed antennas. The present work deals with the design, the fabrication and the characterization of optically transparent and active antennas. The transparent and conducting material used is a micrometric mesh metal film specifically developed, associating high electrical conductivity and high optical transparency. A first optically transparent and miniature FM antenna based on a MESFET transistor with micrometric size has been designed and fabricated. Frequency agile antennas operating in X-band (~10 GHz), based on a beam-lead varactor (agility ~10%) and on a ferroelectric material agility ~2%), have been developed and characterized. An optically transparent and passive antenna has been studied in V-band (~60 GHz). At last, optics (1540 nm) / microwave (1.4 GHz) transition has been performed based on the transmission of a laser beam through the transparent antenna. For all prototypes, an optical transparency level higher than 80% coupled with a sheet resistance value lower than 0.1 ohm/sq have been used.

Les machines électriques sont modélisées pour prédire leurs performances et optimiser leur rendement. Cette modélisation peut être faite par des simulations avec la méthode des éléments finis. En particulier, les machines asynchrones nécessitent des simulations 3D pour prendre en compte les courants de Foucault et les têtes de bobines. Dans le logiciel Flux®, des formulations 3D basées sur le potentiel scalaire magnétique sont utilisées avec succès depuis de nombreuses années. Néanmoins, des coupures mathématiques artificielles sont nécessaires, lorsque le domaine n'est pas simplement connexe.Afin de se libérer de ces contraintes de connexité, des formulations en potentiel vecteur magnétique ont été étudiées et développées. En 3D, leur mise en œuvre nécessite l'utilisation d'éléments finis d'arêtes afin de respecter la nature des champs. Avec les éléments d'arêtes, les formulations sont généralement résolues avec une condition de jauge pour les solveurs directs comme pour les solveurs itératifs. De nouvelles formulations en potentiel vecteur magnétique auto-jaugées ont été développées permettant la prise en compte des bobines maillées et des bobines non maillées. La prise en compte du mouvement est relativement simple à mettre en œuvre pour les formulations en potentiel scalaire magnétique avec l'interpolation nodale.Avec les éléments d'arête, l'interpolation est plus délicate. C'est pourquoi la méthode des éléments avec joints a été développée pour prendre en compte le mouvement dans un cas général. / Electric machines are modeled in order to predict their performance and to optimize their output. This modeling can be done by simulation with the finite element method. In particular, induction machines require 3D simulation to take into account eddy currents and coils overhangs. In the Flux® software, 3D formulations based on magnetic scalar potential has been used with success for many years. Nevertheless, artificial mathematical cuts are necessary, when the domain is not simply connected.In order to avoid connection constraints, magnetic vector potential formulations have been studied and developed. In 3D, their implementation requires the use of edge elements to respect the nature of fields. With edge elements, formulations are generally solved with a gauge condition for direct solvers as well as for iterative solvers. New auto-gauged magnetic vector potential formulations have been developed to take into account meshed coils and non-meshed coils. Consideration of movement is relatively simple to implement for magnetic scalar potential formulations with nodal interpolation. With edge elements, the interpolation is more delicate. For this reason, the mortar method has been developed to take into account movement in a general case.

Ce travail est une contribution à la simulation numérique d'écoulements turbulents autour de corps non profilés. Après avoir précisé les ingrédients numériques et les modèles de turbulence utilisés dans nos simulations, nous présentons une étude sur l'évaluation des effets de la procédure dynamique des modèles de sous-maille dans un modèle VMS-LES et une approche hybride RANS/VMS-LES. Des problèmes d'écoulements autour d'un cylindre seul et en tandem sont considérés. Nous étudions ensuite le comportement de modèles de turbulence hybrides pour la simulation d'écoulements en régime sous-critique autour d'un cylindre circulaire. Le calcul de l'écoulement autour d'un cylindre de section rectangulaire par l'approche VMS-LES est aussi présenté. Enfin, dans une dernière partie, après avoir fait une revue des travaux importants sur les schémas d'avancement en temps multirate, nous proposons une nouvelle approche explicite multirate par agglomération de volumes finis que nous appliquons à des calculs d'écoulements turbulents complexes en utilisant un modèle de turbulence hybride. / This work is a contribution to the numerical simulation of turbulent flows around bluff bodies. After specifying the numerical ingredients and the turbulence models used in our simulations, we present a study on the impact of the dynamic sub-grid scale modeling in VMS-LES model and a RANS/VMS-LES hybrid turbulence approach. Simulations of flows around a cylinder and a tandem are performed. Next, we assess the behaviour of some hybrid turbulence models for the simulation of flows around a circular cylinder in the subcritical regime. The computation of the flow around a rectangular cylinder with the VMS-LES approach is also presented. At last, after a review of some important works on multirate time advancing schemes, we propose a new volume-agglomeration explicit multirate approach that is applied to the computation of complex turbulent flows by a hybrid turbulence model.

Les méthodes d'éléments finis de type Galerkine discontinu (DG FEM) ont démontré précision et efficacité pour résoudre des problèmes d'ondes dans des milieux complexes. Cependant, elles nécessitent un très grand nombre de degrés de liberté, ce qui augmente leur coût de calcul en comparaison du coût des méthodes d'éléments finis continus.Parmi les différentes approches variationnelles pour résoudre les problèmes aux limites, se distingue une famille particulière, basée sur l'utilisation de fonctions tests qui sont des solutions locales exactes des équations à résoudre. L'idée vient de E.Trefftz en 1926 et a depuis été largement développée et généralisée. Les méthodes variationnelles de type Trefftz-DG appliquées aux problèmes d'ondes se réduisent à des intégrales de surface, ce qui devrait contribuer à réduire les coûts de calcul.Les approches de type Trefftz ont été largement développées pour les problèmes harmoniques, mais leur utilisation pour des simulations en domaine transitoire est encore limitée. Quand elles sont appliquées dans le domaine temporel, les méthodes de Trefftz utilisent des maillages qui recouvrent le domaine espace-temps. C'est une des paraticularités de ces méthodes. En effet, les méthodes DG standards conduisent à la construction d'un système semi-discret d'équations différentielles ordinaires en temps qu'on intègre avec un schéma en temps explicite. Mais les méthodes de Trefftz-DG appliquées aux problèmes d'ondes conduisent à résoudre une matrice globale, contenant la discrétisation en espace et en temps, qui est de grande taille et creuse. Cette particularité gêne considérablement le déploiement de cette technologie pour résoudre des problèmes industriels.Dans ce travail, nous développons un environnement Tre#tz-DG pour résoudre des problèmes d'ondes mécaniques, y compris les équations couplées de l'élasto-acoustique. Nous prouvons que les formulations obtenues sont bien posées et nous considérons la difficulté d'inverser la matrice globale en construisant un inverse approché obtenu à partir de la décomposition de la matrice globale en une matrice diagonale par blocs. Cette idée permet de réduire les coûts de calcul mais sa précision est limitée à de petits domaines de calcul. Etant données les limitations de la méthode, nous nous sommes intéressés au potentiel du "Tent Pitcher", en suivant les travaux récents de Gopalakrishnan et al. Il s'agit de construire un maillage espace-temps composé de macro-éléments qui peuvent être traités indépendamment en faisant une hypothèse de causalité. Nous avons obtenu des résultats préliminaires très encourageants qui illustrent bien l'intérêt du Tent Pitcher, en particulier quand il est couplé à une méthode de Trefftz-DG formulée à partir d'intégrales de surface seulement. Dans ce cas, le maillage espace-temps est composé d'éléments qui sont au plus de dimension 3. Il est aussi important de noter que ce cadre se prête à l'utilisation de pas de temps locaux ce qui est un plus pour gagner en précision avec des coûts de calcul réduits. / Discontinuous Finite Element Methods (DG FEM) have proven flexibility and accuracy for solving wave problems in complex media. However, they require a large number of degrees of freedom, which increases the corresponding computational cost compared with that of continuous finite element methods. Among the different variational approaches to solve boundary value problems, there exists a particular family of methods, based on the use of trial functions in the form of exact local solutions of the governing equations. The idea was first proposed by Trefftz in 1926, and since then it has been further developed and generalized. A Trefftz-DG variational formulation applied to wave problems reduces to surface integrals that should contribute to decreasing the computational costs.Trefftz-type approaches have been widely used for time-harmonic problems, while their implementation for time-dependent simulations is still limited. The feature of Trefftz-DG methods applied to time-dependent problems is in the use of space-time meshes. Indeed, standard DG methods lead to the construction of a semi-discrete system of ordinary differential equations in time which are integrated by using an appropriate scheme. But Trefftz-DG methods applied to wave problems lead to a global matrix including time and space discretizations which is huge and sparse. This significantly hampers the deployment of this technology for solving industrial problems.In this work, we develop a Trefftz-DG framework for solving mechanical wave problems including elasto-acoustic equations. We prove that the corresponding formulations are well-posed and we address the issue of solving the global matrix by constructing an approximate inverse obtained from the decomposition of the global matrix into a block-diagonal one. The inversion is then justified under a CFL-type condition. This idea allows for reducing the computational costs but its accuracy is limited to small computational domains. According to the limitations of the method, we have investigated the potential of Tent Pitcher algorithms following the recent works of Gopalakrishnan et al. It consists in constructing a space-time mesh made of patches that can be solved independently under a causality constraint. We have obtained very promising numerical results illustrating the potential of Tent Pitcher in particular when coupled with a Trefftz-DG method involving only surface terms. In this way, the space-time mesh is composed of elements which are 3D objects at most. It is also worth noting that this framework naturally allows for local time-stepping which is a plus to increase the accuracy while decreasing the computational burden.

Les travaux de thèse présentés dans ce mémoire s’inscrivent dans le cadre du projet ANR BARPPRO. Ce programme de recherche vise à étudier l’influence d’une barrière de protection face à une explosion en régime de détonation. L’objectif est d’établir des méthodes de calcul rapides de classement des zones d’effets pour aider les industriels au dimensionnement des barrières de protection. L’une à partir d’abaques, valable pour des configurations en géométrie 2D, sur des plages spécifiées de paramètres importants retenus, avec une précision de +/- 5%. L’autre à partir d’une méthode d’estimation rapide basée notamment sur les chemins déployés, valable en géométrie 2D et en géométrie 3D, mais dont la précision estimée est de +/- 30%. Afin d’y parvenir, l’étude s’appuie sur trois volets : expérimental, simulation numérique et analytique. La partie expérimentale étudie plusieurs géométries de barrière de protection à petites échelles pour la détonation d’une charge gazeuse (propane-oxygène à la stoechiométrie). Les configurations expérimentées servent à la validation de l’outil de simulation numérique constitué du solveur HERA et de la plateforme de calcul TERA 100. Des abaques d’aide au dimensionnement ont pu être réalisés à partir de résultats fournis par l’outil de simulation (3125 configurations de barrière de protection, TNT). L’étude des différents phénomènes physiques présents a également permis de mettre en place une méthode d’estimation rapide basée sur des relations géométriques, analytiques et empiriques. L’analyse de ces résultats a permis d’établir quelques recommandations dans le dimensionnement d’une barrière de protection. Les abaques et le programme d’estimation rapide permettent à un ingénieur de dimensionner rapidement une barrière de protection en fonction de la configuration du terrain et de la position de la zone à protéger en aval du merlon. / This thesis is a part of the ANR BARPPRO project. This research program studies this influence of the protection barrier during an explosion detonation. The goal of this project is to establish fast-computation methods of area classification effects to help the industrial to design the protection barrier on the SEVESO sites. One from abacus, for configurations in 2D geometry on specified parameters used, with an accuracy of +/- 5%. The other from a fast-running method based on broken lines for configurations in 2D and 3D geometries, but the accuracy is +/- 30%. This study includes three approaches: experimental, numerical simulation and analytical approaches. The experimental part studies several geometries of the protection barrier for a gaseous explosion (stoichiometric propane-oxygen mixture) at small scales. The experimental configurations used to validate the numerical simulation tool constituted of the HERA software and the TERA 100 supercomputer. The overpressure charts were able to generate from the numerical results (3125 configurations of the barrier for a TNT charge). The analysis of these results allows to establish different recommendations in the design of the protection barrier. The study of the different physical phenomena present has also helped to set up a fast-running method based on the geometrical, empirical and analytical relations. All these tools will enable an engineer to analyze and estimate the evolution of overpressure around the barrier as a function of the site’s dimensions.

Pour comprendre et diagnostiquer des pathologies telles que l'ostéoporose, qui est un problème majeur de santé publique et, est un facteur de risque important de fractures notamment la Fracture de l'Extrémité Supérieure du Fémur (FESF), il est essentiel d’aborder ces problématiques en trois dimensions (3D) pour fournir au praticien un outil de diagnostic et de dépistage du risque fracturaire. De plus, la visualisation en 3D de l'anatomie joue un rôle important dans le domaine de la chirurgie orthopédique guidée (assistée par ordinateur). En préopératoire, pour la planification chirurgicale ou la conception de prothèses sur mesure, en per opératoire, pour assister le chirurgien durant l'acte chirurgical et enfin en postopératoire pour le suivi. La reconstruction de modèles anatomiques en 3D peut être réalisée par l'utilisation de techniques d'imagerie 3D directes telles que la tomodensitométrie. Cependant, l'utilisation d'une telle imagerie est limitée à des procédures complexes en raison des contraintes imposées par le coût, la disponibilité et les risques de radiation pour le patient. Ainsi, l'alternative à ce type d'imagerie 3D est de développer des méthodes de reconstruction 3D qui s'appuient uniquement sur quelques radiographies 2D. L'objectif principal de cette thèse est de proposer une technique permettant de reconstruire de façon automatique une surface 3D du fémur proximal à partir d'un nombre restreint de radiographies. Les études précédentes sur la reconstruction de surface ont généralement besoin de connaissances supplémentaires comme l'utilisation d'un modèle générique ou statistique 3D de la forme à reconstruire. La méthode décrite dans cette thèse nécessite seulement les coordonnées 3D de points calculées à partir de quelques paires de clichés radiographiques. Deux approches sont proposées. La première méthode repose sur la mise en correspondance de contours extraits de paires de radiographies et sur un modèle mathématique basé sur le principe de la stéréovision pour le calcul d'un nuage de points 3D. La deuxième technique utilise les résultats de l’approche précédente ainsi que des points extraits d’un autre fémur pour améliorer la précision au niveau de certaines régions sensibles choisies par l’opérateur. La reconstruction du modèle surfacique à partir des nuages de points obtenus par les deux techniques est obtenue par un maillage basé sur l'équation de Poisson. Un recalage 3D/3D est effectué entre le nuage de points calculé et le nuage de points extrait du modèle générique connu ("Gold Standard" obtenu avec des coupes CT-Scan) du même fémur proximal afin de pouvoir comparer la surface reconstruite à un modèle "vérité terrain" et ainsi estimer la précision de la méthode. / To understand and diagnose pathologies such as osteoporosis, which is considered as a major public health issue and, an important risk factor of fractures in particular the proximal femur fracture. The classical tools of diagnosis are mainly based on the analysis of X-rays photographs. These techniques have shown many limitations to carry out the key information for the physician. Through the last decade the 3D visualization of anatomy demonstrated the effectiveness for analysis and diagnosis, particularly for the guided orthopedic surgery (computer aided). In preoperative, for the surgical planning or the design of prostheses, in per operative, to assist the surgeon during the surgical act and finally in postoperative for the monitoring. It is also essential to provide the practitioner with a 3D tool for the diagnosis and analysis of the osteoporosis and the fracture risk. Reconstruction of 3D anatomical models can be achieved by the use of direct 3D imaging modalities such as Computed Tomography. However, such technique is limited to complex procedures because of the constraints imposed by cost, availability and risk of radiation to the patient. Thus, the alternative to this kind of 3D imaging is to develop methods for 3D reconstruction which are based only on few 2D radiographs. The main objective of this work is to propose a tool able to reconstruct automatically a 3D surface of the proximal femur from a limited number of X-ray images. Previous studies on the reconstruction of surfaces usually need additional knowledge such as the use of a 3D generic or statistical model of the shape to be reconstructed. The method described in this thesis requires only the 3D coordinates of points calculated from a few pairs of radiographs. Two approaches are proposed. The first method is based on the matching of extracted contours from pairs of radiographs and on a mathematical model based on the principle of the stereovision for the calculation of a 3D point cloud. The second technique uses the results of the previous method as well as new points, chosen by an operator from another proximal femur to improve accuracy at sensitive areas. The reconstruction of the surface model from this cloud of points is obtained by a meshing based on the Poisson's equation. A 3D/3D registration is made between the cloud of the calculated points and the cloud of the extracted points from the generic model ("Gold Standard" obtained with CT-Scan) of the same proximal femur in order to compare the reconstructed surface with a model "ground truth" and thus estimate the accuracy of the method.