matfem - un package de méthode d'éléments finis MATlab
Un package matlab hautement vectorisé orienté objet pour résoudre des modèles d'éléments finis. Ne prend en charge que les problèmes 2D. Développé en 2016b. Je n'ai testé dans aucune autre version de Matlab.
Ce package est une réimplémentation (orientée objet et vectorisée) de ma solution de devoirs pour une classe d'éléments finis. Ce code est développé pour m'aider à comprendre FEM et est
- bon pour : l'analyse numérique, essayer des idées rapidement et échouer rapidement dans un petit modèle,
- pas bon pour : le code de production, ou toute analyse pour laquelle l'algorithme FEM lui-même n'est pas intéressant.
Avantages
- Plus facile à coder que Deal.II
- Niveau inférieur à FEniCS
- Plus rapide que de coder soi-même
- Moins cher que la boîte à outils Matlab PDE
Limites
- 2D uniquement
- Au plus élément de Lagrange de degré 2
- Même si c'est hautement vectorisé, ce n'est pas très rapide
- Ne fonctionne pas en parallèle
Exemples (en construction ?)
Tous les fichiers .xlm
du dossier demos
sont des exemples de codes. Le bonjour tout le monde a_poisson.xlm
est un bon point de départ.
Codes en direct Matlab :
- Équation de Poisson
- Équation de Poisson avec limites de Neumann
Script Matlab (pas encore documenté) :
- Équation de Poisson avec limites de Robin
- Équation de Poisson avec coefficient matriciel
- Équation d'élasticité
- Équation de Stokes à l'état d'équilibre
- Équation de Stokes (en fonction du temps)
- Modèle Stokes couplé à double porosité
Feuille de route