This is a personal project of mine, a sort of proof of concept of what I know about FEM, from the ground up, with no external libraries except for SFML for visualization (and my own libraries). I tend to move most of the general methods (e.g. numerical integration, linear solvers) into my utilities library, so myBasicFemSolver makes extensive use of it.

The focus is FEM for solid mechanics problems, but I am leaving room for other physics and even other discretization methods or domains of computational mechanics (e.g. rigid body, free surface flow), as well as things like meshing and boundary condition definition. Currently, linear FEM for static and dynamic (using Newmark-beta method) 2D solid mechanics works very well. Here are some random pictures and gifs:

Youngs = Density = Horizontal length = 1.0; Poissons = 0.4:

Material properties ~= Aluminium; Horizontal length = 1.0:

Static: