Как правильно написать код конечного элемента

Question

Я хочу разработать свое понимание метода конечных элементов (структурной механики) и программирования на языке Python. По этой причине я хочу создать какое-то простое приложение. По моему мнению, приложение будет своего рода препроцессором и решателем. Хотя процесс слияния будет выполняться в Gmsh и после обработки в Paraview или Mayavi.

Основная проблема - я не начинаю ... Я могу создать простые скрипты python, которые вычисляют фермы или проблемы с лучом, однако матрицы записываются вручную и решаются с помощью SciPy. Моя цель - автоматизировать это и связать с упомянутым программным обеспечением.

Я признаю теорию ООП, сделал пару проектов студентов на C++, но это все, и я хочу двигаться вперед. Спасибо за любой совет!

Answer 1

Если вам не нужно, написание собственного кода конечного элемента является контрпродуктивным.

Попробуйте использовать существующую библиотеку. Я использовал как Deal.II, так и libmesh - оба являются C++ и open-source. (Отказ от ответственности: у меня тоже очень незначительный вклад в обоих). Но есть много вариантов.

Даже если вы в первую очередь хотите лучше понять предмет (в отличие от фактического решения реальных проблем с ним), то, начиная с одной из существующих фреймворков, вы получите кусок инфраструктуры, в которой вы будете нуждаться (Матрицы, решатели, обработка сетки, уточнение сетки). В противном случае вам нужно будет найти и обработать все эти биты самостоятельно.