matfem - um pacote de método de elementos finitos MATlab
Um pacote matlab altamente vetorizado orientado a objetos para resolver modelos de elementos finitos. Suporta apenas problemas 2D. Desenvolvido em 2016b. Não testei em nenhuma outra versão do Matlab.
Este pacote é uma reimplementação (orientada a objetos e vetorizada) da minha solução de lição de casa para uma classe de elementos finitos. Este código foi desenvolvido para me ajudar a entender o FEM e é
- bom para : análise numérica, testar ideias rapidamente e falhar rapidamente em modelos pequenos,
- não é bom para : código de produção ou qualquer análise de que o algoritmo FEM em si não seja de interesse.
Benefícios
- Mais fácil de codificar do que Deal.II
- Nível inferior ao FEniCS
- Mais rápido do que codificar você mesmo
- Mais barato que a caixa de ferramentas Matlab PDE
Limitações
- Somente 2D
- No máximo elemento Lagrange de grau 2
- Mesmo sendo altamente vetorizado, não é muito rápido
- Não funciona em paralelo
Exemplos (em construção?)
Todos os arquivos .xlm na pasta demos são códigos de exemplo. O hello world a_poisson.xlm é um bom lugar para começar.
Códigos ativos do Matlab:
- Equação de Poisson
- Equação de Poisson com limites de Neumann
Script Matlab (ainda não documentado):
- Equação de Poisson com limites de Robin
- Equação de Poisson com coeficiente com valor de matriz
- Equação de Elasticidade
- Equação de Stokes em estado estacionário
- Equação de Stokes (dependente do tempo)
- Modelo acoplado de dupla porosidade-Stokes
Roteiro
