Contribution à la formulation symétrique du couplage équations intégrales - éléments finis : application à la géotechnique, Contributing to the symmetric formulation of the coupling integral equations - finite elements : application to the geotechnics

Sous la direction de Guy Bonnet
Thèse soutenue le 17 septembre 2010: Paris Est
Un des outils numériques les plus utilisés en ingénierie est la méthode des éléments finis, qui peut être mise en o euvre grâce à l'utilisation de nombreux codes de calcul. Toutefois, une difficulté apparaît lors de l'utilisation de la méthode des éléments finis, spécialement en géotechnique, lorsque la structure étudiée est en interaction avec un domaine de dimensions infinies. L'usage courant en ingénierie est alors de réaliser les calculs sur des domaines bornés, mais la définition de la frontière de tels domaines bornés pose de sérieux problèmes. Pour traiter convenablement les problèmes comportant des frontières à l'infini, l'utilisation d'éléments discrets infinis est maintenant souvent délaissée au profit de la méthode des équations intégrales ou méthode des éléments de frontière qui permet de résoudre un système d'équations aux dérivées partielles linéaire dans un domaine infini en ne maillant que la frontière du domaine à distance finie. La mise en oeuvre du couplage entre la méthode des éléments finis et la méthode des éléments de frontière apparaît donc comme particulièrement intéressante car elle permet de bénéficier de la flexibilité des codes de calcul par éléments finis tout en permettant de représenter les domaines infinis à l'aide de la méthode des éléments de frontière. La méthode est basée sur la construction de la matrice de raideur du domaine infini grâce à l'utilisation de la méthode des équations intégrales. Il suffit alors d'assembler la matrice de raideur du domaine infini avec la matrice de raideur du domaine fini représenté par éléments finis. L'utilisation de la méthode la plus simple de traitement des équations intégrales, dite méthode de « collocation » conduit à une matrice de raideur non-symétrique. Par ailleurs, la méthode dite «Singular Galerkin» conduit à une formulation symétrique, mais au prix du calcul d'intégrales hypersingulières. La thèse porte sur une nouvelle formulation permettant d'obtenir une matrice de raideur symétrique sans intégrales hypersingulières, dans le cas de problèmes plans. Quelques applications numériques sont abordées pour des problèmes courants rencontrés en géotechnique
-Équations intégrales
-Éléments finis
-Géotechnique
-Formulation symétrique
-Valeurs propres
-Matrice de raideur
One of the most used numerical tools in engineering is the finite element method, which can be implemented through the use of many computer codes. However, a difficulty arises when using the finite element method, especially in geotechnical engineering, where the structure is studied in interaction with a field of infinite dimensions. The commonly used in engineering is then performming the calculations on bounded domains, but the definition of the border of the domain also poses serious problems. To properly solve the problems which have the boundary at infinity, the use of discrete elements infinite is now often neglected in favor of the integral equations method or boundary element method, which allows to solve a linear partial differential equations system in an infinite domain by the discretization of the only boundary of the domain at finite distance. The implementation of coupling between the finite element method and boundary element method is therefore particularly interesting because it allows to benefit the flexibility of computer codes by the finite element method, while the infinite domains is represented by the help of the integral equations method. It is sufficient to assemble the stiffness matrix of infinite domain with the stiffness matrix of finite domain represented by finite elements. Using the simplest method of treatment of integral equations, known as method of collocation leads to a non-symmetric stiffness matrix. Furthermore, a method known “Galerkin Singular” leads to a symmetric formulation, but it is at the cost of computing hypersingular integrals. The thesis focuses on a new formulation to obtain a symmetric stiffness matrix without full hypersingular, in the case of plane problems. Some numerical applications are discussed for common problems encountered in geotechnical engineering
-Integral equations
-Finite elements
-Geotechnics
-Symmetric formulation
-Eigenvalues
-Stiffness matrix
Source: http://www.theses.fr/2010PEST1081/document