Ce travail a pour objet l'analyse limite des structures par la méthode des éléments finis. Lorsqu'une structure atteint sa charge limite, certaines de ses composantes sont dans la phase inélastique de leur comportement, alors que dans les parties les plus critiques, du fait de la localisation des déformations inélastiques, se produit la rupture du matériau. Les effets de localisation sont, dans les matériaux fragiles liés à l'apparition et au développement de macro fissures alors qu'ils sont, dans les matériaux ductiles, gouvernés par les bandes de cisaillement localisées. L'étude de la charge limite est ainsi reliée à la modélisation du comportement inélastique standard du matériau mais également à la modélisation des effets localisés correspondant au comportement adoucissant des matériaux. Le comportement inélastique standard du matériau est, dans ce travail, décrit par des modèles élastoplastiques, élastoviscoplastiques ou élastiques non linéaires. Tous les modèles de comportement sont définis en termes d'efforts généralisés. Un certain nombre d'approches mathématiques et d'algorithmes numériques sont disponibles mais sont bien souvent inefficaces et manquent de précision. Ainsi, nous utilisons une approche développée plus récemment s'appuyant sur une méthode d'éléments finis enrichis de discontinuités. Nous avons développé de nouvelles formulations d'éléments standards prenant en compte des cinématiques et des descriptions des champs de déplacements discontinus complexes. Plusieurs formulations d'éléments finis ont été développées pour l'analyse de différents composants structurels. Nous présentons, dans un premier temps, un élément fini dédié à l'analyse limite des plaques en béton armé. La formulation d'un élément de plaque élastoplastique et élastoviscoplastique écrite en efforts généralisés associée à une procédure commune d'intégration sont présentées ensuite. Un élément de coque non linéaire, faisant intervenir une fonction seuil à deux surfaces incluant à la fois un écrouissage isotrope et un écrouissage cinématique est ensuite présenté. Les deux derniers éléments finis développés dans ce travail sont dédiés à la modélisation de la rupture localisée dans les poutres planes et les solides bidimensionnels. L'élément de poutre d'Euler-Bernouilli est enrichi par une discontinuité en rotation. Une stratégie s'appuyant sur l'analyse préalable, par un modèle raffiné, d'une partie de la structure est proposée afin d'obtenir les paramètres du modèle constitutif de la poutre. Enfin, nous présentons la formulation d'un élément quadrangulaire à discontinuité forte dont la cinématique permet de prendre en compte des sauts de déplacements linéaires dans les deux directions normale et tangentielle le long de la surface de discontinuité. Des résultats numériques montrent que les éléments développés ainsi que les algorithmes associés constituent un outil efficace et robuste d'analyse de la charge limite et de la rupture des structures. Parmi les exemples, nous présentons la simulation de la propagation d'une fissure dans un matériau fragile ainsi que le développement d'une bande de cisaillement dans un matériau ductile. Les codes numériques associés aux formulations présentées dans ce travail ont été générés par l'outil de programmation symbolique et d'optimisation de code AceGen. Les performances des éléments sont présentés à travers un grand nombre d'exemples numériques réalisés à partir du code AceFem.

Dans cette thèse, nous nous intéressons à la simulation numérique d'écoulements instationnaires de deux fluides visqueux non miscibles, séparés par une interface mobile. Plus particulièrement des écoulements sans choc constitués d'une phase gazeuse et d'une phase liquide sont considérés. Pour modéliser de tels écoulements, une approche dans laquelle le gaz est décrit par les équations de Navier-Stokes compressible et le liquide par les équations de Navier-Stokes incompressible est proposée. C'est le couplage de ces deux modèles qui constitue l'originalité et l'enjeu principal de de cette thèse. Pour traiter cette difficulté majeure, une méthode globale (i.e. la même dans chaque phase) et simple à mettre en \oe uvre est élaborée. L'utilisation des équations de Navier-Stokes formulées de façon unifiée pour les inconnues primitives (pression, vitesse et température) constitue le point de départ pour la construction de notre méthode qui repose sur les composants suivants: - une méthode d'éléments finis stabilisés pour la discrétisation spatiale des équations de Navier-Stokes; - une approche Level Set pour représenter précisément l'interface dont l'équation de transport a été résolue par une méthode de type Galerkin Discontinu; \item et des grandeurs moyennées pour traiter les discontinuités à l'interface. Le bon comportement de notre approche est illustré sur différents tests mono et bi-dimensionnels.

Un nouveau code de simulation numérique pour la propagation des ondes de choc acoustiques dans des géométries complexes a été développé. Le point de départ a été la méthode de Galerkin discontinu qui utilise des maillages non structurés (ici des éléments triangulaires), particulièrement adaptés aux géométries complexes. Cependant, cette discrétisation conduit à l'apparition d'oscillation de Gibbs. Pour pallier ce problème, nous avons choisi d'introduire de la viscosité artificielle au voisinage des chocs. Cela a nécessité le développement de trois outils originaux : (i) un nouveau détecteur de choc sensible aux ondes de chocs acoustiques sur des maillages non structurés, (ii) un nouveau terme de viscosité artificielle dans les équations de l'acoustique non linéaire défini élément par élément et (iii) un nouveau terme permettant de régler le niveau de viscosité locale à partir du raidissement des fronts d'onde. Le code de calcul a été utilisé pour étudier deux configurations différentes. La première concerne la réflexion d'ondes de choc acoustiques sur des surfaces rigides. Différents régimes de réflexion ont alors été observés allant, de la réflexion classique de Snell Descartes jusqu'à celui dit de réflexion faible de Von Neumann. La deuxième configuration était consacrée à la focalisation d'ondes de choc acoustiques produites par un transducteur à haute intensité (comme ceux utilisés en HIFU). Un soin particulier a été pris pour étudier le calcul de l'intensité et pour étudier l'interaction entre les ondes de choc et des obstacles placés dans la région du foyer. / A new numerical solver for the propagation of acoustical shock waves in complex geometry has been developed. This is done starting from the discontinuous Galerkin method. This method is based on unstructured mesh (triangular elements here), and so, naturally it is well-adapted for complex geometries. Nevertheless, the discretization induces Gibbs oscillations. To manage this problem, we choose to introduce some artificial viscosity only in the vicinity of the shocks. This necessitates the development of three original tools. First of all, a new shock sensor for unstructured mesh sensitive to acoustical shock waves has been designed. It senses where the local artificial viscosity has to be introduced thanks to a reformulation of a new element centred smooth artificial viscosity term in the equations. Finally, the amount of viscosity is computed by the introduction of an original notion of gradient factor linked to the steepening of the waveform. The numerical solver has been used to investigate two different physical situations. The first one is the nonlinear reflection of acoustical shock waves on rigid surfaces. Different regimes of reflection have been observed ranging from the linear Snell Descartes reflection to the weak von Neumann case. The second configuration deals with the focusing of shock waves produced by high intensity transducers (like in HIFU). Special attention has been given to the careful computation of intensity and to the interaction between the shock waves and obstacles in the region of the focus.

L'objet de cette thèse est l'étude et la mise en œuvre d'une méthode d’éléments finis multi-échelles pour la simulation d'écoulements miscibles en milieux poreux. La définition des fonctions de base multi-échelles suit l'idée introduite par F. Ouaki. La nouveauté de ce travail consiste à combiner cette approche multi-échelle avec des éléments finis de type Galerkine Discontinus (DG) de façon à pouvoir utiliser ces nouveaux éléments sur des maillages non-conformes composés de mailles de formes diverses. Nous rappelons, dans un premier temps, le principe des méthodes DG et montrons comment ces méthodes peuvent être utilisées pour discrétiser une équation de convection-diffusion instationnaire identique à celle rencontrée dans le problème d'écoulement considéré dans ce travail. Après avoir vérifié l'existence et l'unicité d'une solution à ce problème, nous redémontrons la convergence des méthodes DG vers cette solution en établissant une estimation d'erreur a priori. Nous introduisons, ensuite, les éléments finis multi-échelles non conformes et détaillons leur mise en œuvre sur ce problème de convection-diffusion. En supposant les conditions aux limites et les paramètres du problème périodiques, nous montrons une nouvelle estimation d'erreur a priori pour cette méthode. Dans une seconde partie, nous considérons le problème d'écoulement complet où l'équation considérée dans la première partie est résolue de manière couplée avec l'équation de Darcy. Nous introduisons différents cas tests inspirés de modèles d'écoulements rencontrés en géosciences et comparons les solutions obtenues avec les deux méthodes DG, à savoir la méthode classique utilisant un seul maillage et la méthode étudiée ici. Nous proposons de nouvelles conditions aux limites pour la résolution des problèmes de cellule qui permettent, par rapport à des conditions aux limites linéaires plus classiquement utilisées, de mieux reproduire les variations des solutions le long des interfaces du maillage grossier. Les résultats de ces tests montrent que la méthode multi-échelle proposée permet de calculer des solutions proches de celles obtenues avec la méthode DG sur un seul maillage et de réduire, de façon significative, la taille du système linéaire à résoudre à chaque pas de temps. / This work deals with the study and the implementation of a multiscale finite element method for the simulation of miscible flows in porous media. The definition of the multiscale basis functions is based on the idea introduced by F. Ouaki. The novelty of this work lies in the combination of this multiscale approach with Discontinuous Galerkin methods (DG) so that these new finite elements can be used on nonconforming meshes composed of cells with various shapes. We first recall the basics of DG methods and their application to the discretisation of a convection-diffusion equation that arises in the flow problem considered in this work. After establishing the existence and uniqueness of a solution to the continuous problem, we prove again the convergence of DG methods towards this solution by establishing an a priori error estimate. We then introduce the nonconforming multiscale finite element method and explain how it can be implemented for this convection-diffusion problem. Assuming that the boundary conditions and the parameters of the problem are periodic, we prove a new a priori error estimate for this method. In a second part, we consider the whole flow problem where the equation, studied in the first part of that work, is coupled and simultaneously solved with Darcy equation. We introduce various synthetic test cases which are close to flow problems encountered in geosciences and compare the solutions obtained with both DG methods, namely the classical method based on the use of a single mesh and the one studied here. For the resolution of the cell problems, we propose new boundary conditions which, compared to classical linear conditions, allow us to better reproduce the variations of the solutions on the interfaces of the coarse mesh. The results of these tests show that the multiscale method enables us to calculate solutions which are close to the ones obtained withDG methods on a single mesh and also enables us to reduce significantly the size of the linear system that has to be solved at each time step.

L'étude hydro géologique du massif du Rissiou et de sa couverture sédimentaire a été menée à partir d'observations réalisées sur les émergences de surface et sur les venues d'eau rencontrées pendant le percement des diverses galeries de l'aménagement EDF de Grand'Maison. Après une étude détaillée de la fracturation du massif cristallin du Rissi ou, il est montré qu 'il n 'existe pas une nappe unique mais une juxtaposition de réseaux aquifères, en accord avec l'état de la fissuration et dont les sources constituent les trop pleins: L'hétérogénéité des discontinuités détermine les vitesses d'écoulements , le temps de séjour dans l'encaissant et les caractères physicochimiques et isotopiques des eaux. La juxtaposition de la couverture sédimentaire beaucoup moins permléable, au massif cristallin, se fait par l'intermédiaire d'un accident tectonique majeur localement drainant qui détermine l'indépendance des systèmes hydrogéologiques . Dans la couverture sédimentaire , seul le Trias est apparu aquifère, jouant le rôle de collecteur des eaux superficielles . Les caractères physicochimiques des eaux y sont essentiellement cond itionnées par la nature des terrains traversés et les variations hydroclimatiques de surface.

Le travail présenté ici avait pour but le développement d'un schéma de type Galerkin discontinu (GD) d'ordre élevé pour la résolution des équations de la dynamique des gaz écrites dans un formalisme Lagrangien total, sur des maillages bi-dimensionnels totalement déstructurés. À cette fin, une méthode progressive a été utilisée afin d'étudier étape par étape les difficultés numériques inhérentes à la discrétisation Galerkin discontinue ainsi qu'aux équations de la dynamique des gaz Lagrangienne. Par conséquent, nous avons développé dans un premier temps des schémas de type Galerkin discontinu jusqu'à l'ordre trois pour la résolution des lois de conservation scalaires mono-dimensionnelles et bi-dimensionnelles sur des maillages déstructurés. La particularité principale de la discrétisation GD présentée est l'utilisation des bases polynomiales de Taylor. Ces dernières permettent, dans le cadre de maillages bi-dimensionnels déstructurés, une prise en compte globale et unifiée des différentes géométries. Une procédure de limitation hiérarchique, basée aux noeuds et préservant les extrema réguliers a été mise en place, ainsi qu'une forme générale des flux numériques assurant une stabilité globale L_2 de la solution. Ensuite, nous avons tâché d'appliquer la discrétisation Galerkin discontinue développée aux systèmes mono-dimensionnels de lois de conservation comme celui de l'acoustique, de Saint-Venant et de la dynamique des gaz Lagrangienne. Nous avons noté au cours de cette étude que l'application directe de la limitation mise en place dans le cadre des lois de conservation scalaires, aux variables physiques des systèmes mono-dimensionnels étudiés provoquait l'apparition d'oscillations parasites. En conséquence, une procédure de limitation basée sur les variables caractéristiques a été développée. Dans le cas de la dynamique des gaz, les flux numériques ont été construits afin que le système satisfasse une inégalité entropique globale. Fort de l'expérience acquise, nous avons appliqué la discrétisation GD mise en place aux équations bi-dimensionnelles de la dynamique des gaz, écrites dans un formalisme Lagrangien total. Dans ce cadre, le domaine de référence est fixe. Cependant, il est nécessaire de suivre l'évolution temporelle de la matrice jacobienne associée à la transformation Lagrange-Euler de l'écoulement, à savoir le tenseur gradient de déformation. Dans le travail présent, la transformation résultant de l'écoulement est discrétisée de manière continue à l'aide d'une base Éléments Finis. Cela permet une approximation du tenseur gradient de déformation vérifiant l'identité essentielle de Piola. La discrétisation des lois de conservation physiques sur le volume spécifique, le moment et l'énergie totale repose sur une méthode Galerkin discontinu. Le schéma est construit de sorte à satisfaire de manière exacte la loi de conservation géométrique (GCL). Dans le cas du schéma d'ordre trois, le champ de vitesse étant quadratique, la géométrie doit pouvoir se courber. Pour ce faire, des courbes de Bézier sont utilisées pour la paramétrisation des bords des cellules du maillage. Nous illustrons la robustesse et la précision des schémas mis en place à l'aide d'un grand nombre de cas tests pertinents, ainsi que par une étude de taux de convergence. / The intent of the present work was the development of a high-order discontinuous Galerkin scheme for solving the gas dynamics equations written under total Lagrangian form on two-dimensional unstructured grids. To achieve this goal, a progressive approach has been used to study the inherent numerical difficulties step by step. Thus, discontinuous Galerkin schemes up to the third order of accuracy have firstly been implemented for the one-dimensional and two-dimensional scalar conservation laws on unstructured grids. The main feature of the presented DG scheme lies on the use of a polynomial Taylor basis. This particular choice allows in the two-dimensional case to take into general unstructured grids account in a unified framework. In this frame, a vertex-based hierarchical limitation which preserves smooth extrema has been implemented. A generic form of numerical fluxes ensuring the global stability of our semi-discrete discretization in the $L_2$ norm has also been designed. Then, this DG discretization has been applied to the one-dimensional system ofconservation laws such as the acoustic system, the shallow-water one and the gas dynamics equations system written in the Lagrangian form. Noticing that the application of the limiting procedure, developed for scalar equations, to the physical variables leads to spurious oscillations, we have described a limiting procedure based on the characteristic variables. In the case of the one-dimensional gas dynamics case, numerical fluxes have been designed so that our semi-discrete DG scheme satisfies a global entropy inequality. Finally, we have applied all the knowledge gathered to the case of the two-dimensional gas dynamics equation written under total Lagrangian form. In this framework, the computational grid is fixed, however one has to follow the time evolution of the Jacobian matrix associated to the Lagrange-Euler flow map, namely the gradient deformation tensor. In the present work, the flow map is discretized by means of continuous mapping, using a finite element basis. This provides an approximation of the deformation gradient tensor which satisfies the important Piola identity. The discretization of the physical conservation laws for specific volume, momentum and total energy relies on a discontinuous Galerkin method. The scheme is built to satisfying exactly the Geometric Conservation Law (GCL). In the case of the third-order scheme, the velocity field being quadratic we allow the geometry to curve. To do so, a Bezier representation is employed to define the mesh edges. We illustrate the robustness and the accuracy of the implemented schemes using several relevant test cases and performing rate convergences analysis.

Ce travail est consacré au développement de schémas numériques pour approcher les solutions de modèles d'écoulement type “depth averaged”. Dans un premier temps, nous détaillons la construction d'approches Volumes Finis pour le système Shallow Water avec termes sources sur maillages non structurés. En se basant sur une reformulation appropriée des équations, nous mettons en place un schéma équilibré et préservant la positivité de la hauteur d'eau, et suggérons des extensions MUSCL adaptées. La méthode est capable de gérer des topographies irrégulières et exhibe de fortes propriétés de stabilité. L'inclusion des termes de friction fait l'objet d'une analyse poussée, aboutissant à l'établissement d'une propriété type “Asymptotic Preserving” à travers l'amélioration d'un autre récent schéma Volumes Finis. La seconde composante de cette étude concerne les méthodes Elements Finis type Galerkin discontinu. Certaines des idées avancées dans le contexte Volumes Finis sont employées pour aborder le système Shallow Water surmaillages triangulaires. Des résultats numériques sont exposés et la méthode se révèle bien adaptée à la description d'une large variété d'écoulements. Partant de ces observations nous proposons finalement d'exploiter ces caractéristiques pour étendre l'approche à une nouvelle famille d'équations type Green-Nadghi. Des validations numériques sont également proposées pour valider le modèle numérique. / This work is devoted to the development of numerical schemes to approximatesolutions of depth averaged flow models. We first detail the construction of Finite Volume approaches for the Shallow Water system with source terms on unstructured meshes. Based on a suitable reformulation of the equations, we implement a well-balanced and positive preserving approach, and suggest adapted MUSCL extensions. The method is shown to handle irregular topography variations and demonstrates strong stabilities properties. The inclusion of friction terms is subject to a thorough analysis, leading to the establishment of some Asymptotic Preserving property through the enhancement of another recent Finite Volume scheme.The second aspect of this study concerns discontinuous Galerkin Finite Elementmethods. Some of the ideas advanced in the Finite Volume context areemployed to broach the Shallow Water system on triangular meshes. Numericalresults are exposed and the method turns out to be well suited to describe a large variety of flows. On these observations we finally propose to exploit its features to extend the approach to a new family of Green-Nadghi equations. Numerical experiments are also proposed to validate this numerical model.

Ce travail vise à améliorer la description à l'échelle du nanomètre des sols argileux expansifs en utilisant la théorie de la fonctionnelle de densité (DFT). L’eau n’est plus considérée comme un solvant continu mais comme un fluide de molécules polaires individuelles. L'objectif est de reproduire les résultats issus de l'expérience ou de la modélisation numérique tels que la présence de couches d'eau discrètes dans l'espace interfolaire ou la variation de la pression de disjonction avec la distance interfolaire dans le régime de gonflement cristallin. Différents phénomènes physiques de complexité croissante sont successivement étudiés. La taille finie des molécules d'eau est tout d'abord prise en compte en modélisant l'eau comme un fluide de sphères dures traité par la théorie fondamentale de la mesure. La nature polaire du solvant est ensuite implicitement considérée en utilisant un potentiel intermoléculaire de Lennard-Jones pour reproduire les différents types d'interactions de Van der Waals. La nature dipolaire de l'eau est ensuite explicitement modélisée par un fluide dipolaire de sphères dures. Ces deux derniers modèles utilisent une approche perturbative de la théorie de la fonctionnelle de densité dans laquelle les effets de corrélation entre les molécules du fluide sont incorporés. Les ions sont finalement ajoutés afin de compléter la description de la double couche électrique. En vue d'une application au génie civil, l'expression améliorée de la pression de disjonction à l'échelle nanométrique est incluse dans une forme modifiée du principe de Terzaghi appliqué aux argiles expansives non-saturées récemment développée dans notre groupe afin de simuler numériquement le comportement hydro-mécanique des argiles gonflantes lors d’essais d'infiltration d’eau / This work aims at improving the nanoscale description of expansive clayey soils using the Density Functional Theory (DFT). Water is no longer considered as a continuous solvent but as a fluid of individual polar molecules in order to recover existing experimental and modeling results such as the presence of discrete water layers in the interplatelet space or the variation of the disjoining pressure with the interplatelet distance at low hydration level. Different physical phenomena of increasing complexity are successively considered. The finite size of the water molecules is firstly taken into account by modeling water as a Hard Sphere fluid using the Fundamental Measure Theory. The polar nature of the water solvent is then implicitly taken into account through a Lennard-Jones potential averaging the different types of Van der Waals interactions. Next the polar nature of the solvent is explicitly modelized by considering water as a Dipolar Hard Sphere fluid. These two fluid models are studied in the framework of the Density Functional Perturbation Theory in which correlation effects between the fluid molecules are incorporated. Ions are finally added in order to complete the Electrical Double Layer description at the nanoscale. With the objective of an application to civil engineering, the improved expression of the disjoining pressure at the nanoscale is included in a modified form of Terzaghi's effective stress principle for unsaturated expansive clays recently developed by our group in order to numerically simulate the hydro-mechanical behavior of expansive clays during water uptake

Mes travaux de thèse portent sur la construction de conditions aux limites absorbantes (CLAs) pour des problèmes de propagation d'ondes posés dans des milieux limités par des surfaces régulières. Ces conditions sont nouvelles car elles prennent en compte non seulement les ondes proagatives (comme la plupart des CLAs existantes) mais aussi les ondes évanescentes et rampantes. Elles sont donc plus performantes que les conditions existantes. De plus, elles sont facilement implémentables dans un schéma d'éléments finis de type Galerkine Discontinu (DG) et ne modifie pas la condition de stabilité de Courant-Friedrichs-Lewy (CFL). Ces CLAs ont été implémentées dans un code simulant la propagation des ondes acoustiques ainsi que dans un code simulant la propagation des ondes en régime harmonique. Les comparaisons réalisées entre les nouvelles conditions et celles qui sont les plus utilisées dans la littérature montrent que prendre en compte les ondes évanescentes et les ondes rampantes permet de diminuer les réflexions issues de la frontière artificielle et donc de rapprocher la frontière artificielle du bord de l'obstacle. On limite ainsi les coûts de calcul, ce qui est un des avantages de mes travaux. De plus, compte tenu du fait que les nouvelles CLAs sont écrites pour des frontières quelconques, elles permettent de mieux adapter le domaine de calcul à la forme de l'obstacle et permettent ainsi de diminuer encore plus les coûts de calcul numérique. / In my PhD, I have worked on the construction of absorbing boundary conditions (ABCs) designed for wave propagation problems set in domains bounded by regular surfaces. These conditions are new since they take into account not only propagating waves (as most of the existing ABCs) but also evanescent and creeping waves. Therefore, they outperform the existing ABCs. Moreover, they can be easily implemented in a discontinuous Galerkin finite element scheme and they do not change the Courant-Friedrichs-Lewy stability condition. These ABCs have been implemented in two codes that respectively simulate the propagation of acoustic waves and harmonic waves. The comparisons performed between these ABCs and the ABCs mostly used in the litterature show that when we take into account evanescent and creeping waves, we reduce the reflections coming from the artificial boundary. Therefore, thanks to these new ABCs, the artificial boundary can get closer to the obstacle. Consequently, we reduce the computational costs which is one of the advantages of my work. Moreover, since these new ABCs are written for any kind of boundary, we can adapt the shape of the computational domain and thus we can reduce again the computational costs.

Les structures mécaniques utilisées de nos jours ne cessent d’évoluer en utilisant des matériaux composites ou à gradient fonctionnel afin de répondre aux enjeux de résistance accrue, allégement de la structure et amélioration des performances. Ceux-ci nécessitent un contrôle adéquat de leur état de santé afin de s’assurer de l’intégrité de la structure. L’utilisation des ondes guidées ultrasonores fournit un moyen efficace et rapide d’inspection sur de longues distances. Néanmoins, ces ondes présentent certaines caractéristiques complexes qui rendent la tâche très difficile. L’utilisation d’outils d’analyse tels que les modèles numériques constitue un grand atout pour ce type d’application. Dans ce contexte, l’objectif de cette de thèse est le développement d’un outil de modélisation performant, permettant d’étudier la propagation des ondes guidées ultrasonores avec une grande précision et une faible consommation de ressources et de temps de calculs. De ce fait, l’intérêt est porté sur des méthodes numériques d’ordres élevés dont les propriétés de convergence sont beaucoup améliorées que les méthodes classiques. En particulier, la méthode semi-analytique éléments finis de Galerkin discontinue pour la détermination des courbes de dispersion des ondes guidées est développée. La méthode est applicable aux structures planes et cylindriques fabriquées de matériaux isotropes, anisotropes et hétérogènes (à gradient fonctionnel de propriétés). Une étude comparative sur l’analyse des performances de ces méthodes est effectuée. Celle-ci a démontré la capacité de la méthode à modéliser la propagation des ondes guidées ultrasonores dans des guides d’ondes à section arbitraire avec des performances prometteuses par rapport à la méthode des éléments finis classique. / The mechanical structures used today are constantly evolving using composite or functionally gradient materials to meet the challenges of increased strength, lightening the structure and improving performance. These require adequate control of their state of health to ensure the integrity of the structure. The use of Ultrasound Guided Waves (UGW) provides an efficient and fast way of inspection over long distances. Nevertheless, these waves have some complex features that make the task very difficult. The use of analysis tools such as numerical models is a great asset for this type of application. In this context, the objective of this thesis is the development of a powerful modeling tool, allowing to study the propagation of UGWs with a great precision, less computational time and consumption of resources. Accordingly, we are interested in higher order numerical methods whose convergence properties are much improved than the classical methods. In particular, a semi-analytical discontinuous Galerkin finite element method (SADG-FE) is developped for the determination of the dispersion properties of guided waves in arbitrary cross-section waveguides. The method is applicable to plates and cylindrical structures made of isotropic, anisotropic heterogeneous (functionally graded) materials. The performance analysis of these methods and their comparisons are performed with respect to the models based on the classical finite element method. The results demonstrated the ability of the proposed method to model the propagation of ultrasounic guided waves in arbitrary section waveguides with promising performance over the conventional finite element method.