Les fichiers présents sur ce dépôt permettent d'évaluer la solution d'un problème mis sous forme d'équation différentielle sur un domaine de dimension 2D définit par sa frontière de telle sorte que la solution est nulle sur celle ci. L'utilisateur sera requis d'entrer les valeur de sa choix pour la configuration initiale du script, qui serviront à déterminer le type d'éléments finis à utiliser (NQ1/NQ2), le nombre de points de Gauss pour l'intégration numérique et le nombre d'éléments suivant chaque axe du domaine 2D.
- Maillage du domaine 2D selon le type d'éléments et leurs nombre.
- Calcul les éléments relatifs à la quadrature de Gauss.
- Evaluation les fonctions de forme.
- Génération les matrices élémentaires et second membre.
- Assemblage des matrices élémentaires.
- Calcul de la solution.
Permet de générer un maillage sur un domaine 2D définit par les paramètres suivants :
- Frontière sur chaque axe.
- Type d'éléments finis.
- Nombre d'éléments suivant chaque axe.
Exemples
Fonction qui donne les poids et les noeuds de Gauss en 2D suivant le type d'élément et le nombre de points de calcul de l'intégration numérique de Gauss.
Evaluation des fonction chapeaux, leur forme globale a étée déterminé par Maple. Ce script calcule les valeur des coefficients des fonctions de forme.
Génération des matrices élémentaires de rigidité et second membre en utilisant les fonctions chapeaux et les évaluant sur chaque élément du maillage.
Assemblage des matrices élémentaires en matrices globales sur la totalité du domaine.
Regroupe les étapes et initie le calcul de la solution évaluée aux noeuds du maillage du domaine.
Affiche la solution interpolée sur un maillage secondaire pour éviter les problèmes d'affichage triangulaire au lieu d'un affichage rectangulaire.
Script principal du programme, contient toutes les entrées et l'enchainement des étapes ainsi que l'affichage final et l'évaluation de l'erreur.
- Domaine :
(-2 2)x(-1 1) - Type d'éléments :
nq1 - Points de Gauss :
4 - Nombre d'éléments :
8x8
-
Approchée
-
Exacte
| Nombre d'éléments total | Erreur norme L2 | Erreur relative |
|---|---|---|
| 64 | 6.4076500339013243e-01 | 3.7666761480123384e+01 |
| 256 | 3.1571845434976814e-01 | 1.5289501206979500e+02 |
| 1024 | 1.5727399187148874e-01 | 6.1385635955881662e+02 |
| 4096 | 7.8563778713928500e-02 | 2.4577139001903406e+03 |