Calculation of high degree shape functions by a perturbation technique
Calcul de fonctions de forme de haut degré par une technique de perturbation
Résumé
Most problems of physics and mechanics lead to partial differential equations. The many methods that exist are relatively low degree. In this thesis, we propose a method of very high degree. Our idea is to increase the order of interpolation function via a perturbation technique to avoid or reduce the difficulties caused by the high cost operations such as integrations. In dimension 1, the proposed technique is close to the P-version finite elements. At a basic level, approximates the solution by a power series of order p. In the case of a linear equation of order 2, the local resolution can build an element of degree, with two degrees of freedom per element. For nonlinear problems, a linearization of the problem by Newton's method is needed. Tests involving linear and nonlinear equations were used to validate the method and show that the technique has a similar convergence in the p-version finite elements. In dimension 2, the problem is discretized through reorganizing polynomials in homogeneous polynomials of degree k. After a definition of variables called principal and secondary combined with a vertical scanning field, the problem becomes a series of 1D problem. A collocation technique allows to take into account the boundary conditions and coupling conditions and determine the solution of the problem. The collocation technique coupled with the least-squares enabled to improve the initial results and has made more robust the perturbation technique
La plupart des problèmes de la physique et de la mécanique conduisent à des équations aux dérivées partielles. Les nombreuses méthodes qui existent déjà sont de degré relativement bas. Dans cette thèse, nous proposons une méthode de très haut degré. Notre idée est d'augmenter l'ordre des fonctions d'interpolation via une technique de perturbation afin d'éviter ou de réduire les difficultés engendrées par les opérations très coûteuses comme les intégrations. En dimension 1, la technique proposée est proche de la P-version des éléments finis. Au niveau élémentaire, on approxime la solution par une série entière d'ordre p. Dans le cas d'une équation linéaire d'ordre 2, cette résolution locale permet de construire un élément de degré élevé, avec deux degrés de liberté par élément. Pour les problèmes nonlinéaires, une linéarisation du problème par la méthode de Newton s'impose. Des tests portant sur des équations linéaires et nonlinéaires ont permis de valider la méthode et de montrer que la technique a une convergence similaire à la p-version des éléments finis. En dimension 2, le problème se discrétise grâce à une réorganisation des polynômes en polynômes homogènes de degré k. Après une définition de variables dites principales et secondaires associé à un balayage vertical du domaine, le problème devient une suite de problème 1D. Une technique de collocation permet de prendre en compte les conditions aux limites et les conditions de raccord et de déterminer la solution du problème. La collocation couplée avec la technique des moindres carrés a permis de d'améliorer les premiers résultats et a ainsi rendu plus robuste la technique de perturbation