Esse repositório contempla os códigos utilizados na disciplina COC752 da COPPE/UFRJ que aborda o Método dos Elementos Finitos. A disponibilização do código aqui presente tem o objetivo de facilitar futuras consultas nessa área, assim como auxiliar futuros estudantes da disciplina.
A disciplina de Métodos dos Elementos Finitos 1 possui um trabalho que deve ser desenvolvido ao longo do curso. Um exemplo é fornecido, mas o trabalho se divide em algumas partes:
O professor sugere que algumas malhas sejam desenvolvidas para a realização do trabalho. Essas malhas podem variar de 16x16 até 512x512 elementos. Os elementos incluem triângulos e quadriláteros. Para entender o formato de entrada do programa, verifique aqui.
Essa parte do trabalho não precisa ser realizada necessariamente em Fortran, existem precedentes em Python e usando o Excel. O importante é entender a lógica de como a malha é gerada para os quadriláteros e depois como é feita a divisão do quadrilátero em 2 elementos triangulares.
Um exemplo de implementação em Fortran está disponível na pasta gerador_malhas. Um exemplo em Python apenas para quadriláteros feito pelo ex-aluno Fabio Cesar Canesin está disponível aqui.
Caso você nem saiba por onde começar, um tutorial sobre como configurar um ambiente de desenvolvimento está disponível aqui.
O programa fornecido não suporta o visualizador Paraview. A segunda parte consiste em adaptar o programa para que sua saída seja convertida no formato compatível com o visualizador. Uma alternativa é fazer um conversor entre os formatos.
Para entender melhor o formato de saída do programa do trabalho, verifique aqui. Para entender como é estruturado o formato VTK, veja a documentação geral disponível aqui.
O programa fornecido não suporta o elemento quadrilátero bilinear. A terceira parte consiste na inclusão desse elemento.